Nonlinear stability of expanding stars in the mass-critical Euler-Poisson system
Hadzic, Mahir
The gravitational Euler-Poisson system is a fundamental astrophysics model of a Newtonian star. We first give a brief overview of the existing results on the free-boundary compressible Euler-Poisson system. We then study the question of nonlinear stability of homogeneous expanding star-solutions discovered by Goldreich and Weber in 1980's in the mass-critical gravitational Euler-Poisson system. We show that these solutions are nonlinearly stable with respect to small perturbations. We thus construct a new class of global-in-time solutions, which are not homogeneous and therefore not encompassed by the existing works.The problem is mass-critical with respect to an invariant rescaling and the nonlinear analysis is carried out in suitably chosen similarity coordinates. We present some interesting open questions at the end. This is a joint work with Juhi Jang
 I would first like to thank the organizers for this wonderful conference. It's the first time I'm here. So before I say anything about math, everything I will talk about is a joint work with Juhi Jang who is based at the University of Southern California and this year at the IAS. So the good thing about this long title is that it basically tells you the statement of the main theorem. There are certain special solutions to the system of PDEs I will look at. Euler Poisson system, it comes from astrophysics, it's a model of a star. And we will be studying the nonlinear stability of that class of solutions. So let me give you an outline. Well first, I intend to spend quite a bit of time introducing the system and talking about general properties, perhaps at the expense of giving you some technical details in the second half of the talk. But I can promise that there will be two different stability results pertaining to two different types of special solutions that we will look at. And one of them, the first one is easier. This is the one where I will probably give you a pretty complete idea how the proof works and hopefully there will be time to talk about the second type of expansion. So before any of that, let me just tell you that the gravitational Euler Poisson system is the model you would find in a typical astrophysics book describing a star. So it is a Newtonian model. The idea is that the star is a gas or a fluid kept together by gravitational forces. So let me tell you what this talk is not about. It's not about relativistic models of stars, but nevertheless part of the motivation for studying this model indeed comes from relativity. And I will mention that hopefully at the end of the talk. So as I said, the star is idealized as a compact fluid body. And the word compact is crucial because it suggests that there is a boundary between the fluid and the vacuum. And if you want to model such a problem and if you want to track the behavior of the boundary, you have to allow for that boundary to move. So this causes, this gives us a moving star vacuum interface. So what is the model? Well, first I'll tell you the unknowns. The first three are unsurprising and you're used to them. The fluid density rho and the pressure P, the velocity field U, and the gravitational field phi that you need to couple to the Euler equations. The fourth unknown or the further unknown is the support of the fluid because we are trying to treat a free boundary problem. Okay, and the equations, again, you have seen some fluid models this time around. Don't be scared. We will soon forget about this slide. The first equation is the continuity equation. The second one expresses the conservation of momentum. And on the right hand side, you will see a term caused by the gravitational field phi. And the gravitational field phi is self-consistent. It's generated by the mass density here. So Laplace phi is 4 pi rho. This is the Poisson equation and these are the equations of compressible fluid dynamics. On top of that, we have to supply this problem with some boundary conditions and remember this boundary is moving. So the kinematic boundary condition is the natural statement that the boundary is moved by the fluid particles. So the normal velocity of the boundary is, in fact, the normal component of the velocity field at the free boundary. And the pressure vanishes at the boundary. Okay? Now, if you count the unknowns, you will notice that there is one superfluous unknown here, namely the pressure. There is no evolution equation for pressure and you need to prescribe a further equation to or a relation to actually close this system. There is one equation for rho, three for u, one for phi, but there is nothing for pressure. And this is typically done by prescribing the so-called equation of state. Now equations of state is something that in the physics literature, you would look at a particular model of a star and you would hope to have a good understanding of what the equation of state is for that particular star. What is commonly done, people study the so-called polytropic equation of state where p of rho is rho to a certain power gamma. And this gamma is the so-called adiabatic index and it's allowed to vary between one and two. So the way I sort of would like you to think about this model, I would like to think of this nonlinear problem as parameterized by gamma. Okay? And as we vary gamma, the qualities, the qualitative properties of the solutions of this system will change. Okay. So let me comment a little bit on the literature. The vacuum interface itself is known to cause severe analytical difficulties in improving local, even local well-positiveness. So the basic result in the theory of compressible fluids with vacuum interfaces are basically these two results, they came independently in 2010 by Couton and Scholar and Jang and Masmoudi. And so what they have to contend with is a certain degeneracy. I will explain it as we go through the talk that is caused by the presence of the vacuum. But one of the key insights is the use of the Lagrangian coordinates. So if you use Lagrangian coordinates and if you rephrase the compressible Euler equations surrounded by vacuum, you will discover that there is a wave equation, a quasi-linear wave equation lurking in the background. Okay? So that's somehow the one of the big revelations. It is a degenerate wave equation, degenerate in the sense that I will specify later, but the word degenerate has to do with the fact that some of the coefficients of the wave equation will vanish at the vacuum boundary. And another key assumption, so another key assumption in these works is that of a physical vacuum boundary condition. So here I'm stating what it is. It's a statement that the fluid enthalpy, which is this quantity, rho to the power gamma minus 1, has to have its normal derivative at the vacuum boundary, has to have a sign. This is perhaps not a well-motivated condition. So physically it expresses the fact that particles have to accelerate in the normal direction at the boundary. But I will soon actually motivate this, sorry. I will soon motivate this condition by finding particular steady-state solutions of the Euler-Pollisson system that satisfy this. And indeed, if you open an astrophysics textbook, the steady-state star models that they look at all satisfy this type of condition, the compactly supported ones. Okay, so there's a little bit of a discrepancy here. Here I'm citing works that pertain to a compressible Euler equation. You could ask me, what about Euler-Pollisson? Because I'm adding a gravitational field. From the point of view of well-pollinist theory, this is a lower-order perturbation. So if you understand how to prove well-pollinist here, you will understand how to prove well-pollinist for Euler-Pollisson system. But nevertheless, there is a reference that I'm missing. Where is the chalk? Oh, fantastic. Thank you. So prior to this work, there's one result which deals only with a priori estimates with the Lienblad, Couton-Lienblad and Scholar. So this is one missing. So an a priori, Bounds, Couton-Lienblad, Scholar. And for the Euler-Pollisson system, this came much later, entirely based on these ideas. There's a work of Lu-Shin and Zhang. Okay, this is for you. Not for the vacuum interface. So for the liquid case, when the density is not vanishing at the boundary. So this is a significantly less degenerate model. Okay, so instead of, I will now try to give you, to do three things for you. I will give you an example of this famous class of steady states known as Lien-Emden-Stars. Then I will, I will, I'm expecting a sigh of relief. I will produce an invariant scale, rescaling for this problem and explain a certain notion of criticality that enters and that will motivate the main results. So most famous class of special solutions are the steady states. And the way you produce them, you set rho to be only a function of x and u to be zero. You plug it in and the whole equation reduces to this statement here, gradient of rho star gamma plus rho star gradient phi star equals zero. So if you use the fact that this satisfies the Poisson equation and if you call this object here, this is the fluid enthalpy, if you call it w, then in spherical symmetry, and this is a very classical stuff, you will see this in astrophysics textbooks, you reduce the problem to solving this ordinary differential equation. And the key question is, for what powers of gamma can I find a compactly supported solution with finite mass? This is your star. So the answer is, if gamma lies in this magical range, then you can find profiles that are indeed compactly supported, say if it supports zero, one, and they look exactly like that. Once gamma hits six-fifths and below, no such solutions exist. It's just the ODE analysis. And moreover, these special solutions in fact satisfy the physical vacuum conditions. So w prime at one is strictly less than zero. So the density profile is not smooth across the vacuum interface. And this is one of the critical points that works of Couton, Scholar, Jang and Masmoudi had to address, the fact that the non-smoothness is a generic part of the theory. So you have to find a way, what does it even mean to differentiate the equation in a proper sense. You can ask however, and this is what physicists did, are these solutions stable? So now the numerology gets even more interesting. If gamma lies in this range here, the lane end and stars are linearly unstable. And if gamma is bigger, equal than four-thirds, it's in fact linearly stable. So it's just some basic linear analysis. Of course it begs a question about two questions. Are these results, can they be upgraded to a nonlinear statement? So the first result can, it is not terribly surprising. Once you have a growing mode, which is indeed the main source of instability in this regime, and once you have a good well-posedness theory, then you can prove nonlinear instability. And this was shown by Jang in 2014. Nonlinear stability is a much more contentious point. So there is a conditional nonlinear stability result in a fairly weak topology generated by basically just the energy by Ryan from 2003, which says if the solution exists, it will be nonlinearly stable. However, there is absolutely no reason why the solution should exist, why small perturbations of these steady states should. And it's in the same regime, four-thirds to two. Pardon? Four? The result of Ryan, is it in that regime? Right, so the result of Ryan, the topology is so weak that it doesn't really see the free boundary. In particular, this type of stability result, if you think about it, will actually allow you to have support splitting up and things like that, the topology will not see it. So you really need something else. And in fact, I will talk about it at the very end. I will highlight this as one of the important open problems in this development. Okay, so let me try to perhaps explain these exponents that show up in a different way. And this hopefully will help you appreciate them. There is an invariant rescaling of this problem. To each gamma, you can associate this invariant rescaling where you see where you just work it out and you realize that if rho and u are solutions, then rho tilde and u tilde are also solutions. And the way you rescale the time depends on gamma. So really the value of gamma has a lot to say about the structure of the solutions of the equation. So similar to the pressure and the potential, you can work it out from these how to rescale. Okay, the second ingredient that I'm sure you will appreciate, what are the conserved quantities? Well, the mass, this is easy, the integral of the density, and the total energy. Let's just briefly talk about the total energy. This bit here is the kinetic energy. This object here, remember, rho is Laplace phi. So if you integrate by parts, this becomes negative gradient phi squared. So this is just the gravitational energy, which gives a negative contribution. And this thing here, rho to the power gamma, is some sort of thermal energy of the star. Okay, and the question that is quite natural in this context is what is the behavior under the rescaling? And then you discover that the mass remains invariant exactly when gamma equals 4 thirds and energy remains invariant exactly when gamma equals 6 fifths. And these are the two numbers that are popped up here. So this is to suggest that these numbers, this numerology has a special meaning. All of the results for the remainder of the article of this presentation will pertain to the mass critical case gamma equals 4 thirds. I said article because the actual article is supposed to appear tomorrow on the pre-print server on the archive. And there are several spots where we said and erased all the results for the remainder of this article pertain to this case. Anyway, but okay, this is the key point to remember. So gamma equals 4 thirds is the mass critical case. I'll tell you why physicists like this particular value of gamma. You see, it's a rescaling that keeps the mass preserved. So it's a very appealing to think of gravitational collapse or expansion as some process that sort of cascades through these scales, lambda, but it preserves the mass. It's very natural. So the beautiful thing about this is that in fact there is a result of Goldreich and Weber who were physicists in 1980 where they found explicit examples of collapsing and expanding solutions of examples of stars in this mass critical case. And now let me do some computation. So if gamma equals 4 thirds and if my self-similar rescaling is this, so you can work out what this coefficient is here. So 2 minus 4 thirds is 2 thirds. One over that is 3 halves. So this becomes P over lambda to the power 3 halves. So to have a self-similar collapse or a self-similar rescaling would correspond to having a collapsing solution whose radius scales like T minus T to the power 2 thirds or an expanding solutions that expands at the same rate. This would be the corresponding coefficient if you wanted to call this a self-similar collapse. However, if you follow through the paper of Goldreich and Weber, already implicitly there and then explicitly in the works of Makino and Lin, you will find a second family of expanding and collapsing solutions that expand and collapse at the linear rate. And this was part of the mystery for us to understand, first of all, the reasons why this is so. Yes? Sorry. Which one corresponds to the explicit solution? This one. This is not explicit. This is just, you see, I'm putting this sign here. So this is informal. In a few slides, I will give you the exact formula. There is a correction term. For this guy, there is no correction term. This is perfect. This is corrected. By your log? No. No, no, no. It's just something linear. You will see. You will see. No logs. Well, we'll come to that at the end of the talk as well. So there are two families of such a family? There are two families. No, no, no. It's not very self-similar. I mean this. But there are two families of collapsing and expanding solutions. So for whatever reasons, in these works, you will never read this off. They will never talk about the rates. They analyze this and they say there exist these solutions. But if you really look at the rates, this is what will happen. I will give a more precise statement later on. So this is an informal statement. Informal statement. Don't ask me for any norms and so on at this moment. This will also come later. The basic result is that these self-similar expanding profiles, so I'm not talking about collapse, the ones that behave like this, like 2 to the power 2 thirds, are co-dimension one stable. And the linearly expanding guys will be just stable. In this case, we can actually talk about asymptotic stability in this co-dimension one set of perturbations. Whereas here, we do not characterize the asymptotic attractors. There are just some objects that do not necessarily belong to some sort of nearby linearly expanding profile. OK. This is the gist of the results. I will get into the technical statement slightly later. In particular, these are, to my knowledge at least, first non-trivial examples of global solutions for this free boundary problem. And what is very important for us, we can actually say that the support of the star, of the nearby perturbed star actually grows at approximately the same rate as the underlying background solution. OK. Now let me take you back to what ideas go into the well-posedness theory. And this will bring us then to the proofs of these results. As I mentioned, fundamental idea in this work is to use the Lagrangian coordinates. OK. So you've seen them. Eta dot, OK. Eta is simply the flow of particle trajectories. And because I will now make an assumption of spherical symmetry, it is actually not fundamental to these results. But I want to make it because it will display in a very concise way the fundamental structure that we need to worry about. OK. So assume everything is spherically symmetric. Now the, maybe it's slightly pompous to say structural miracle, but it is a beautiful structure of these equations. If you use Lagrangian coordinates in spherical symmetry, you discover that chi, which is now my Lagrangian map, satisfies a second order equation, chi-t-t plus f-w of chi equals zero, where this f-w is some nonlinear operator which is of second order. OK. And this is an effective quasi-linear wave equation on a compact domain. OK. Now there's this thing here. I don't want to display f-w because it's not beautiful. But I will tell you later how it linearizes. And you will be convinced that it is a second order, it's a second order self-adjoint operator in a certain space. Is this an artifact of the fact that this would be irrotational? Correct. Absolutely. Yeah. So I can, OK. So now that you raised this question. So if it were not spherically symmetric or irrotational. If it were pure, if it were the full system. The flow map eta would satisfy an effective wave map type equation. The divergence of the flow map, which roughly corresponds, which captures the irrotational part, would satisfy any equation like this. So divergence of the jacobo, divergence of the eta t of the velocity vector field. As the curl would satisfy, because, OK, so let me explain this. Because remember in the Lagrangian coordinates rho dt u, where u is now pullback with respect to the t is my Lagrangian time and u is my Lagrangian velocity. Plus the pullback of the pressure would be something. And your curl of this guy vanishes. So this would give you, if you wish, a transport equation for curl and the sort of estimates of curl of u for free. So the idea is that when you deal with the full problem, you estimate the divergence in the curl in two separate ways. So you have a, you can decouple. Absolutely. Yeah. So this is in some sense at the heart of the approach of Couton, Scholar, Jeng and Nassmoudi. The framework, the functional framework is slightly different, but this is the main idea. OK, so there is a wave equation. That's my point. And it's posed on a compact domain, because you pulled it back onto the, to fix the free boundary, you had to pull it back onto something fixed. So it's a wave equation on a compact domain. So let me now re-derive, really, those solutions of Goldreich and Weber and other people. I will call them homogeneous solutions, because I will make the ansatz that chi is just lambda of t. It doesn't depend on r. Now this will tell you, this ansatz will reduce the previous equation to this thing here, lambda square, lambda double dot, plus something which doesn't depend on t equals zero. So you can separate variables, and if you do that, you discover that there is an ODE satisfied by lambda. And there is another ODE in r, this time satisfied by w. You are interested in finding compactly supported solutions, because that's what qualifies as a star with a vacuum interface, and you can do that. There is a magical value of delta star, such that for any delta bigger than delta star, you can find compactly supported solutions to this equation. And this is just an ODE. You can really, for a range of values of delta and initial conditions, you can really classify and see what are the rates of collapse or expansion and so on. So I will do that for you in a second. As I said, you can solve this ODE, prescribe the initial conditions. There are three parameters in the problem, the delta that I mentioned, lambda 1, initial velocity, and lambda 0, the initial radius of the star. There is a conserved energy associated with this ODE. It's trivial. It's this quantity here, lambda dot squared plus 2 delta over lambda. In the Eulerian description, so you can always go back and recover the original solution in the Eulerian description, the density and the velocity will look like this. So the key point, the density is like lambda to the minus 3. So if lambda is shrinking, density is blowing up. If lambda is expanding, the density, the fluid wants to tear itself apart. It's decaying. And remember, as I mentioned, self-similarity would correspond to these numbers here, 2 thirds. So it turns out that self-similar solutions exist, but only when the energy is zero. If energy is not zero, they're not going to obey the 2 third law. And instead of solving the ODE, I also could not resist not putting this on the slides. This is a kind of a bifurcation diagram that explains what's going on. So this is a delta axis. This is the lambda 1 axis. And I'm fixing initial radius to be exactly 1. So what you see here, the light blue region, this is the range of parameters that will give you always a linearly expanding profile. A profile that behaves asymptotically like some constant times t. These guys here, the thick blue line, this is exactly the self-similar expansion. And there is an explicit formula that I will give you for them on the next slide. So you have to insist that the energy is zero. So they're obviously unstable, because if you can go either way from the surface, the question is, if you look at perturbations of energy zero, are they going to remain stable? And then you have the collapsing. So as you can see, most of the collapsing solutions in this range here exhibit self-similar collapse. And then there is this lower dimension line where the collapse happens at a linear rate. So this is the structure of this space in this two-dimensional parameter space. And the basic results are this light blue region is stable, and this thick blue region is a co-dimension 1 stable. So for this line, you have an explicit formula, which is just a statement that E equals zero. Do you have a question? You mean here? Yes. Well, I can tell you a very weak statement. I can prove that these are non-linearly stable in this self-similar frame, but this is not a good statement. You really need to, I don't know whether, so for instance, I don't know if you perturb an element here, whether it will collapse. Okay, so this is an open problem. I will state it at the end as well. Is a linear rate of expansion the physical one, the one that has to be fixed? So Goldreich and Weber completely discarded this for whatever reason. And they were interested in collapse and expansion at self-similar rate. If you look at other physics literature, they actually mostly don't want to deal with the free boundary. They are aware of the difficulties that come with that. So in their literature, they think of Goldreich and Weber as the only example where you have a sharp boundary. In most cases, they are worried about constructing some type of solution that exhibits either collapse or supernova or the expansion, where you have infinite support, but some sort of decay at infinity. This would be for them good enough. And this also leads to good problems because most of it is also not rigorous. So part two, how to study the stability of the expanding guys. This is the formula that Frank, I believe, asked me about. This is the explicit formula for self-similar expansion, whereas you don't have explicit formulas for the linear expansion. They are sort of approximate, but they are the leading word or just this. So I will from now on focus only on the expanding profiles and perhaps pompously call it supernova, but let's just say expansion. And I want you to recall that they satisfy a wave equation of this type, and FW is some second-order operator. Now how to study this object? Well, you adapt your unknown to what you believe to be stable by dividing by it. Basically you define psi to be key, he divided by lambda bar. And then you expect psi to be a small perturbation of one, which happens to be a steady state of this equation in this case. So if you make this ansatz and plug it back in here, you get this equation, second-order PDE. Now how to deal with this? Well, self-similar expansion really behaves well with respect to the scaling invariances of the problem. So the natural thing to do is to rescale time with respect to the corresponding self-similar rescaling. As you will see, the linear expansion doesn't sort of, the equations don't digest it very well. So it will cause certain difficulties that will not be present in the self-similar case. Okay, so let's move on to the, are there any questions about this? So is it absolutely clear what's going on? Okay. So let me tell you the basic structure that will allow you to understand the stability. So let's now focus purely on the self-similar profiles. So remember, let me write it down. Self-similar profile is something that expands like t to the 2 thirds. So this is the radius of the star. Okay, XE is a steady state of the equation on this slide. So I will write XE as 1 plus a perturbation. XE equal 1 is a steady state. You do that, you demand that the initial energy is zero, that the energy of perturbation is zero. So this is a conserved quantity. And of course the key idea, as I mentioned, is to pass to the self-similar time. If you do that, so this roughly means that the new time behaves like log t. Okay, because lambda expands like 2 thirds. If you do that, you discover that this coefficient here is constant. I call it B to honor some of the notation that's popular in the literature. This quantity is a constant which is strictly less than zero. What about the non-linearity? So you expand the non-linearity around 1. And you discover, and this is now the point that will perhaps clarify one of the questions, you discover that in the linearization of this operator is precisely this second order elliptic operator here. But now I want you to stare at it for a second. Notice that when r is zero, this term vanishes. And when r is one, remember w delta is the density profile, is a density profile which roughly looks like this. And it also vanishes at r equal 1. So in this sense, this operator is degenerate because it causes, it has these vanishing weights, if you wish, at zero and at one. This is not surprising. This is the type of difficulty you will see in the study of just pure Euler. But it is suggestive. It tells you what sort of functional framework to use to, in fact, deal with this problem. It tells you that you need to use some type of weighted spaces. Now the cool thing about this, if you pass through the self-similar time and you express, you now try to find the equation for phi, this is what you discover. So phi ss, the second derivative of phi, plus a strictly positive term times phi s minus something negative times phi, plus this linearized operator equals the right-hand side, which is nonlinear. Okay, and you immediately see, you realize what should be the mechanism for stability. It is precisely this blue term here which acts as an effective damping. This is a consequence of the fact that we are expanding our profile around an expanding solution. If I did it around a collapsing solution, the sign here would change and this would be sort of the opposite of damping. So it's a wave equation, quasi-linear wave equation on a compact domain, but with a this damping term. So the correct statement is that there is a high-order energy. If it is initially sufficiently small, and if the energy, the physical energy initially is exactly zero, then there is a global unique solution to the above equation and in fact it decays exponentially fast in the s variable. In the t variable, this translates into some sort of algebraic decay. How? The higher the energy is bound to the quesit area, I guess, right? Precisely. This is what you need. This thing eats you the derivative, right? Thank you for the question because I will now, in next slide, I will write down the energy. The L delta is a self-adjoint operator with respect to this weight. Remember this vanishes at one, this vanishes at zero. Spectral gap, okay, one is in the null space of this operator. So if you're orthogonal to the null space, you're good. You have a quercivity of this operator. This per se is obvious. What I do want to pay attention to, this first term is just the L2 norm in the first term in the corresponding weighted space. If you want to control the first derivative, you have to raise the level of degeneracy in the weight W delta. So we interpret it as first term. Correct. Correct. So this is a critical observation. No, it's only in this direction. Absolutely. Excuse me, the B you mentioned is, sorry, forget, is a constant? It's a constant. Okay. So B is a constant due to the, because of the self-similar structure. Because I'm expanding around self-similar solution. So it really captures that. This is the correct information about the self-similarity that you need to address this. Term is obviously the good guy, this term here. This is the bad guy. This will give you exactly one unstable mode. But this is not surprising. I told you that it has to be unstable. Where the self-similar expansion was a co-dimension one phenomenon. So clearly this has to be the case. But this unstable mode, which is caused by this negative term, is, well, remember, we are starting with zero energy data. And so the unstable direction you can check this is transversal to this zero energy surface. So you can really control the instability, which translates in this case in controlling the inner product of phi and the generator of the null space. So this allows you to carry on the estimates. So think of this as, so this is the relationship that you get by linearizing energy equals zero. This is the quantity that you will control just by doing some simple energy identity. This is the quantity that you want to control to get a spectral gap. Now answer to Pierre's question. Indeed, the energy is motivated in this case quite concretely, in the second case more philosophically, in the case of linear expansion, by the works of these people. And now we come to your observation. Indeed, I don't want you to gain any intuition from this. All I want you to see is that each time I take a further normal derivative, I have to raise some index in my weight. And this is how this works. These delta k spaces basically mean you add k, you add a k-fold power of w delta to your weight. And this is the only way you can do it. So more spatial derivatives implies more degenerate weights. And now you will trust me if I tell you how the energy method works. You try to prove something like this. This term here kind of comes from that damping term that I mentioned. Of course, you need to control the energy by this dissipative term here. You can do that because of the spectral gap property. It's a technical thing, but you can do it. And this, in principle, gives you the exponential decay and completes the proof. But perhaps I want to say, going back one slide, this is a mixed space time norm. In reality, you do the estimates only with the time vector field because it commutes with the equation entirely. And then you use elliptic estimate to control the higher space derivatives with the time derivatives. So this is what I just said. Use the time derivatives to build an energy, elliptic estimates. Now what are the technical tools that go into closing the estimates? Remember there is a nonlinear right hand side which actually looks quite dirty when you write it down. So the tools are hardy inequalities which exactly allow you to deal with those weights that degenerate at 0 and 1 and embeddings between weighted sub-allows spaces. If you take enough derivatives, you can control lower norm by the energy, L infinity norms of lower order terms by the energy. So I will not go into any of this. This is technical. I have eight minutes left. I will, OK, so let me then say just a few words about this second problem here and maybe try to hint at least at why is it more difficult in the first one. The difficulty appears not to be conceptual but more technical. But it is there and it took us some time to understand. So the reason why this problem is harder conceptually is because you cannot, it doesn't honor the self-similar structure of the problem. Doesn't honor the scaling in variances. These solutions are there but how do you prove that they are stable? So it will necessitate the use of a, so first of all, the fact that it is not honoring the invariance will disallow us from using time vector fields. As you will see, they will not commute well with the problem because now when I write down the equation for my perturbation phi, it will have time dependent coefficients. OK, so we have to come up with a correct spatial derivative that will capture the degeneracies at the vacuum. And this is precisely this vector field here, which just happens to be a five dimensional aplosion but this is a coincidence. OK, so what we do, we rescale time so that tau, the new time becomes log t again. So this means dividing by this linearly expanding profile. Let me write it down. OK, here. So lambda tilde of t, think of it as something that behaves like a constant times t. And if you write down the equation for the perturbation, now you discover that you have these tau dependent factors in front and lambda tilde is growing exponentially in this new variable. OK, so you have e to the, let's say e to the tau plus e to the tau plus this stuff equals right inside. This is still causing a certain damping effect. It is well known that the expansion produces a sort of a stabilizing effect. It's reminiscent of works, for instance, in a special relativity, in the general relativity when you have cosmological constant, for instance. But, and this is the statement, there is a high order energy e tilde that will allow you that will give you global existence for small data. However, you cannot classify the asymptotic attractor. So all you can say is that the time derivative, the tau derivative of the solution decays, but you cannot say that in the limit you belong to a nearby member of the linearly expanding family. This energy is constructed by using powers of this elliptic operator that I just wrote a second ago. And so to be honest, I can talk about it. I have a few more slides at the end, but let me not go into the technical aspects of this proof. Let me just say that this operator is designed to capture very precisely. So in other words, you cannot just blindly apply dr derivative or dx. If you do that, you will create singularities at zero and they are hard to handle. So you have to do it with a particular combination. And this operator s is a reflection of the structure of that weighted-generate operator L delta. And let's leave it at that. So I will not go into the technical details. Yeah, that's a good question. So it's just something that sits nearby, but it's not one of the... But you have to identify the proof, right? Yeah, yeah. Well, you can prove that the time derivative of that is decaying. So that gives you convergence to something. Yeah, yeah. And it's nearby. Yeah. So the reminiscent of a result I heard recently in a talk by Falker, Falker Schuhe, and sort of in these expanding space times. Who's here, I believe? Yeah, very good. Let me spend the remaining three minutes discussing perhaps questions that are... I find very interesting that are open. So here's this kind of a silly lemma. I call it lemma. It's too ambitious to call it lemma. Do there exist collapsing or expanding compactly supported stars? So this is the key for other values of gamma. They should exist in the supercritical range, which is between 6 fifths and 4 thirds. This is an open question. What you can show, however, is that self-similarity and a spherical symmetry cannot coexist for this type of collapse. Lama being 4 thirds is absolutely critical for the two to be able to combine the two. And this is a very silly thing to show. It's just some conservation law type calculation. So the question is, can you construct exisymmetric self-similar collapse in this mass supercritical range? I think this is a very interesting problem. And now we come to the question that I believe Sergio asked. Are the Lane-Emdon stars that are mentioned before in the subcritical range, where they are shown to be conditionally stable, are they in fact stable? So a possible mechanism, and you will see this discussed sometimes in the physics literature, is some type of shock formation. This is not known, and this is really important. Let me remind you of this picture. The pink elephant in the room is, of course, this region here, the collapsing region. So as I said, it is not difficult to show that it's nonlinearly unstable. The moment you have correct well-posedness theory and you have a growing mode, it's very easy. But so again, this should be a lemma. But it's really unsatisfactory, because it doesn't tell you anything about the real instability. In particular, it could it be that the nearby guys decay, collapse at perhaps a corrected rate? Do they collapse at all? So this is just in spherical symmetry. You don't have to go outside spherical symmetry to ask these questions. And finally, I would like to sort of point out this viewpoint that I believe studying this family of problems is really the right thing, because you really want to have a robust understanding of how equations of state affect the collapse. This is partly motivated by general relativity. If you think of the most famous example of a collapsing star in general relativity, it's the Oppenheimer-Sneider solution, which corresponds to the case where you don't have any pressure. So by a result of Christodulu from 83, I believe, you can show that these are generically unstable solutions. If you go outside the homogeneity class, you generically form naked singularities. So you want to, well, you can devise various conclusions out of that. One conclusion that I like to draw out of that is that studying non-trivial equations of state is an important problem, and this is just the Newtonian version of the problem. But realistic equations of state are, this is just a final comment, is in principle quite a hard problem, even for physicists. So this is my favorite relativistic astrophysics book by Zeldovich and Noviko, where I sort of learned this stuff about. And they discuss one of the sort of high density regimes for the equations of state. And here's what Zeldovich says. So it reminds him of Avocenko's parody. So it is precisely like this. So if you look at the astrophysical literature, there's a huge disparity between what people think is the right equation of state. Yeah. So this is a teaser, I think of it as a teaser trailer for this wonderful textbook because it's a very serious book, but as I browse through it every now and then I find a gem of this type. Thank you. So the question is, how about? So it seems, since the non-year expanded rate is unstable, does it mean that basically the type of solutions which would converge to the stable one, which is a new one? Okay, so if you generically perturb a self-similar expanding rate, if you generically perturb it with perturbations that have positive energy, then they will converge to the linear guy. So it's generically unstable and these guys are stable really only in the co-dimension one sense. So the linear is stable. The linear is stable. The linear takes over, not the self-similar. Yeah. Yeah, yeah, exactly. So do the physicists know this part? I mean, I'm not... It should be of great interest to them, I guess. Well, if you're... If you're very well, I should be linear. Yeah, exactly. So if you read Goldreich and Weber, which was one of the starting points for this project, for some reason they really say, we will not be interested in these linear rates. They drop a constant of integration in a second order ODE, which allows them to go down one level and get an explicit solution, which is this. And they literally say, ah, we kind of don't want to look at that. So but this is the finding. I mean, part of the mess was also to understand the exact structure of this parameter space. Yes. The non-linear unstable result covers the whole regime of gamma written there or...? Actually it actually goes up to here. I did not tell you what happens at 6-fifths. At 6-fifths you still get a steady state, but it is supported all the way up to... It's infinite. It has an infinite support. And it turns out to be... So in this case, the steady state is exactly the steady state of the critical Quintic wave equation. So it's one of these guys. I don't know. Once there's some constant here, and then I believe one over two, something like that. So this is the rate of decay of this guy. So it is the most complicated case to prove the most important? In this work, well, it's quite different because here you are working on a compactly supported domain. So you really have to rely on the well-posed in this theory for vacuum interface problems. Here it's different. Here the weights are at infinity. So it has a flavor of it, but it's quite different. I'm not sure how to... So this result is also due to Jang. So this equality case, all of this is due to Jang. Sorry, in this range. I have to be careful. This is the instability range. What this result shows in a way is that the fourth-third case, which is linearly stable, is in fact widely unstable. So your result is radial from the stability. This is entirely radial. So the paper that will show up is entirely radial, but it holds also for non-radial perturbations. So I alluded earlier to that. Here's the reason why the radial case is simpler, not just for obvious reasons, but also because it reduces the whole problem to just a quasi-linear wave equation. When you work with a non-radial case, so the fluid is not irrotational anymore in general, so you have to account for the vorticity, which does not have this wave-like structure, but it satisfies a certain transport equation that allows you to control it. What makes the vorticity better than in the standard equation where you know that you cannot control the vorticity? Is it because of the? Well, it's a local in-time result. I mean, yeah. No, but is it draw? Yeah, yeah. So all I'm saying, I'm not sure if I understand your question, but all I'm saying is that Ah, the long time. For the long time result, you don't have a long time. Oh, we haven't produced it yet. No, no. But for the irrotational case, it will be straightforward. And for the vorticity in the long time result, I don't want to make the claim yet. I don't think- Because I don't see why it should be any better than in the regular case. Maybe. No, no, no. You have this damping effect, so it will be better. It will be better. But will it affect the vorticity? No. You have, for instance, for Friedman, Lomet, Roberts and Walker space times, you have a stability without irrotational assumption. So when you have a cosmological constant, for instance, and you perturb away from that, you see this. You can control the full flow. Of course, this is not a free boundary problem, because this free boundary problem in GR is more- Okay, but the mechanism is massive. Yeah, yeah. You will kill it with that. All right. So thank you. So we start again at 11.30. So let's go to the next one. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you.
 I would first like to thank the organizers for this wonderful conference.  It's the first time I'm here.  So before I say anything about math, everything I will talk about is a joint work with Juhi  Jang who is based at the University of Southern California and this year at the IAS.  So the good thing about this long title is that it basically tells you the statement  of the main theorem.  There are certain special solutions to the system of PDEs I will look at.  Euler Poisson system, it comes from astrophysics, it's a model of a star.  And we will be studying the nonlinear stability of that class of solutions.  So let me give you an outline.  Well first, I intend to spend quite a bit of time introducing the system and talking  about general properties, perhaps at the expense of giving you some technical details in the  second half of the talk.  But I can promise that there will be two different stability results pertaining to  two different types of special solutions that we will look at.  And one of them, the first one is easier.  This is the one where I will probably give you a pretty complete idea how the proof works  and hopefully there will be time to talk about the second type of expansion.  So before any of that, let me just tell you that the gravitational Euler Poisson system  is the model you would find in a typical astrophysics book describing a star.  So it is a Newtonian model.  The idea is that the star is a gas or a fluid kept together by gravitational forces.  So let me tell you what this talk is not about.  It's not about relativistic models of stars, but nevertheless part of the motivation for  studying this model indeed comes from relativity.  And I will mention that hopefully at the end of the talk.  So as I said, the star is idealized as a compact fluid body.  And the word compact is crucial because it suggests that there is a boundary between  the fluid and the vacuum.  And if you want to model such a problem and if you want to track the behavior of the boundary,  you have to allow for that boundary to move.  So this causes, this gives us a moving star vacuum interface.  So what is the model?  Well, first I'll tell you the unknowns.  The first three are unsurprising and you're used to them.  The fluid density rho and the pressure P, the velocity field U, and the gravitational  field phi that you need to couple to the Euler equations.  The fourth unknown or the further unknown is the support of the fluid because we are  trying to treat a free boundary problem.  Okay, and the equations, again, you have seen some fluid models this time around.  Don't be scared.  We will soon forget about this slide.  The first equation is the continuity equation.  The second one expresses the conservation of momentum.  And on the right hand side, you will see a term caused by the gravitational field phi.  And the gravitational field phi is self-consistent.  It's generated by the mass density here.  So Laplace phi is 4 pi rho.  This is the Poisson equation and these are the equations of compressible fluid dynamics.  On top of that, we have to supply this problem with some boundary conditions and remember  this boundary is moving.  So the kinematic boundary condition is the natural statement that the boundary is moved  by the fluid particles.  So the normal velocity of the boundary is, in fact, the normal component of the velocity  field at the free boundary.  And the pressure vanishes at the boundary.  Okay?  Now, if you count the unknowns, you will notice that there is one superfluous unknown here,  namely the pressure.  There is no evolution equation for pressure and you need to prescribe a further equation  to or a relation to actually close this system.  There is one equation for rho, three for u, one for phi, but there is nothing for pressure.  And this is typically done by prescribing the so-called equation of state.  Now equations of state is something that in the physics literature, you would look at  a particular model of a star and you would hope to have a good understanding of what  the equation of state is for that particular star.  What is commonly done, people study the so-called polytropic equation of state where p of rho  is rho to a certain power gamma.  And this gamma is the so-called adiabatic index and it's allowed to vary between one  and two.  So the way I sort of would like you to think about this model, I would like to think of  this nonlinear problem as parameterized by gamma.  Okay?  And as we vary gamma, the qualities, the qualitative properties of the solutions of  this system will change.  Okay.  So let me comment a little bit on the literature.  The vacuum interface itself is known to cause severe analytical difficulties in improving  local, even local well-positiveness.  So the basic result in the theory of compressible fluids with vacuum interfaces are basically  these two results, they came independently in 2010 by Couton and Scholar and Jang and  Masmoudi.  And so what they have to contend with is a certain degeneracy.  I will explain it as we go through the talk that is caused by the presence of the vacuum.  But one of the key insights is the use of the Lagrangian coordinates.  So if you use Lagrangian coordinates and if you rephrase the compressible Euler equations  surrounded by vacuum, you will discover that there is a wave equation, a quasi-linear wave  equation lurking in the background.  Okay?  So that's somehow the one of the big revelations.  It is a degenerate wave equation, degenerate in the sense that I will specify later, but  the word degenerate has to do with the fact that some of the coefficients of the wave  equation will vanish at the vacuum boundary.  And another key assumption, so another key assumption in these works is that of a physical  vacuum boundary condition.  So here I'm stating what it is.  It's a statement that the fluid enthalpy, which is this quantity, rho to the power gamma  minus 1, has to have its normal derivative at the vacuum boundary, has to have a sign.  This is perhaps not a well-motivated condition.  So physically it expresses the fact that particles have to accelerate in the normal direction  at the boundary.  But I will soon actually motivate this, sorry.  I will soon motivate this condition by finding particular steady-state solutions of the Euler-Pollisson  system that satisfy this.  And indeed, if you open an astrophysics textbook, the steady-state star models that they look  at all satisfy this type of condition, the compactly supported ones.  Okay, so there's a little bit of a discrepancy here.  Here I'm citing works that pertain to a compressible Euler equation.  You could ask me, what about Euler-Pollisson?  Because I'm adding a gravitational field.  From the point of view of well-pollinist theory, this is a lower-order perturbation.  So if you understand how to prove well-pollinist here, you will understand how to prove well-pollinist  for Euler-Pollisson system.  But nevertheless, there is a reference that I'm missing.  Where is the chalk?  Oh, fantastic.  Thank you.  So prior to this work, there's one result which deals only with a priori estimates with  the Lienblad, Couton-Lienblad and Scholar.  So this is one missing.  So an a priori, Bounds, Couton-Lienblad, Scholar.  And for the Euler-Pollisson system, this came much later, entirely based on these ideas.  There's a work of Lu-Shin and Zhang.  Okay, this is for you.  Not for the vacuum interface.  So for the liquid case, when the density is not vanishing at the boundary.  So this is a significantly less degenerate model.  Okay, so instead of, I will now try to give you, to do three things for you.  I will give you an example of this famous class of steady states known as Lien-Emden-Stars.  Then I will, I will, I'm expecting a sigh of relief.  I will produce an invariant scale, rescaling for this problem and explain a certain notion  of criticality that enters and that will motivate the main results.  So most famous class of special solutions are the steady states.  And the way you produce them, you set rho to be only a function of x and u to be zero.  You plug it in and the whole equation reduces to this statement here, gradient of rho star  gamma plus rho star gradient phi star equals zero.  So if you use the fact that this satisfies the Poisson equation and if you call this  object here, this is the fluid enthalpy, if you call it w, then in spherical symmetry,  and this is a very classical stuff, you will see this in astrophysics textbooks, you reduce  the problem to solving this ordinary differential equation.  And the key question is, for what powers of gamma can I find a compactly supported solution  with finite mass?  This is your star.  So the answer is, if gamma lies in this magical range, then you can find profiles that are  indeed compactly supported, say if it supports zero, one, and they look exactly like that.  Once gamma hits six-fifths and below, no such solutions exist.  It's just the ODE analysis.  And moreover, these special solutions in fact satisfy the physical vacuum conditions.  So w prime at one is strictly less than zero.  So the density profile is not smooth across the vacuum interface.  And this is one of the critical points that works of Couton, Scholar, Jang and Masmoudi  had to address, the fact that the non-smoothness is a generic part of the theory.  So you have to find a way, what does it even mean to differentiate the equation in a proper  sense.  You can ask however, and this is what physicists did, are these solutions stable?  So now the numerology gets even more interesting.  If gamma lies in this range here, the lane end and stars are linearly unstable.  And if gamma is bigger, equal than four-thirds, it's in fact linearly stable.  So it's just some basic linear analysis.  Of course it begs a question about two questions.  Are these results, can they be upgraded to a nonlinear statement?  So the first result can, it is not terribly surprising.  Once you have a growing mode, which is indeed the main source of instability in this regime,  and once you have a good well-posedness theory, then you can prove nonlinear instability.  And this was shown by Jang in 2014.  Nonlinear stability is a much more contentious point.  So there is a conditional nonlinear stability result in a fairly weak topology generated  by basically just the energy by Ryan from 2003, which says if the solution exists, it  will be nonlinearly stable.  However, there is absolutely no reason why the solution should exist, why small perturbations  of these steady states should.  And it's in the same regime, four-thirds to two.  Pardon?  Four?  The result of Ryan, is it in that regime?  Right, so the result of Ryan, the topology is so weak that it doesn't really see the  free boundary.  In particular, this type of stability result, if you think about it, will actually allow  you to have support splitting up and things like that, the topology will not see it.  So you really need something else.  And in fact, I will talk about it at the very end.  I will highlight this as one of the important open problems in this development.  Okay, so let me try to perhaps explain these exponents that show up in a different way.  And this hopefully will help you appreciate them.  There is an invariant rescaling of this problem.  To each gamma, you can associate this invariant rescaling where you see where you just work  it out and you realize that if rho and u are solutions, then rho tilde and u tilde are  also solutions.  And the way you rescale the time depends on gamma.  So really the value of gamma has a lot to say about the structure of the solutions of  the equation.  So similar to the pressure and the potential, you can work it out from these how to rescale.  Okay, the second ingredient that I'm sure you will appreciate, what are the conserved  quantities?  Well, the mass, this is easy, the integral of the density, and the total energy.  Let's just briefly talk about the total energy.  This bit here is the kinetic energy.  This object here, remember, rho is Laplace phi.  So if you integrate by parts, this becomes negative gradient phi squared.  So this is just the gravitational energy, which gives a negative contribution.  And this thing here, rho to the power gamma, is some sort of thermal energy of the star.  Okay, and the question that is quite natural in this context is what is the behavior under  the rescaling?  And then you discover that the mass remains invariant exactly when gamma equals 4 thirds  and energy remains invariant exactly when gamma equals 6 fifths.  And these are the two numbers that are popped up here.  So this is to suggest that these numbers, this numerology has a special meaning.  All of the results for the remainder of the article of this presentation will pertain  to the mass critical case gamma equals 4 thirds.  I said article because the actual article is supposed to appear tomorrow on the pre-print  server on the archive.  And there are several spots where we said and erased all the results for the remainder  of this article pertain to this case.  Anyway, but okay, this is the key point to remember.  So gamma equals 4 thirds is the mass critical case.  I'll tell you why physicists like this particular value of gamma.  You see, it's a rescaling that keeps the mass preserved.  So it's a very appealing to think of gravitational collapse or expansion as some process that  sort of cascades through these scales, lambda, but it preserves the mass.  It's very natural.  So the beautiful thing about this is that in fact there is a result of Goldreich and  Weber who were physicists in 1980 where they found explicit examples of collapsing and  expanding solutions of examples of stars in this mass critical case.  And now let me do some computation.  So if gamma equals 4 thirds and if my self-similar rescaling is this, so you can work out what  this coefficient is here.  So 2 minus 4 thirds is 2 thirds.  One over that is 3 halves.  So this becomes P over lambda to the power 3 halves.  So to have a self-similar collapse or a self-similar rescaling would correspond to having a collapsing  solution whose radius scales like T minus T to the power 2 thirds or an expanding solutions  that expands at the same rate.  This would be the corresponding coefficient if you wanted to call this a self-similar collapse.  However, if you follow through the paper of Goldreich and Weber, already implicitly there  and then explicitly in the works of Makino and Lin, you will find a second family of  expanding and collapsing solutions that expand and collapse at the linear rate.  And this was part of the mystery for us to understand, first of all, the reasons why  this is so.  Yes?  Sorry.  Which one corresponds to the explicit solution?  This one.  This is not explicit.  This is just, you see, I'm putting this sign here.  So this is informal.  In a few slides, I will give you the exact formula.  There is a correction term.  For this guy, there is no correction term.  This is perfect.  This is corrected.  By your log?  No.  No, no, no.  It's just something linear.  You will see.  You will see.  No logs.  Well, we'll come to that at the end of the talk as well.  So there are two families of such a family?  There are two families.  No, no, no.  It's not very self-similar.  I mean this.  But there are two families of collapsing and expanding solutions.  So for whatever reasons, in these works, you will never read this off.  They will never talk about the rates.  They analyze this and they say there exist these solutions.  But if you really look at the rates, this is what will happen.  I will give a more precise statement later on.  So this is an informal statement.  Informal statement.  Don't ask me for any norms and so on at this moment.  This will also come later.  The basic result is that these self-similar expanding profiles, so I'm not talking about  collapse, the ones that behave like this, like 2 to the power 2 thirds, are co-dimension  one stable.  And the linearly expanding guys will be just stable.  In this case, we can actually talk about asymptotic stability in this co-dimension  one set of perturbations.  Whereas here, we do not characterize the asymptotic attractors.  There are just some objects that do not necessarily belong to some sort of nearby linearly expanding  profile.  OK.  This is the gist of the results.  I will get into the technical statement slightly later.  In particular, these are, to my knowledge at least, first non-trivial examples of global  solutions for this free boundary problem.  And what is very important for us, we can actually say that the support of the star,  of the nearby perturbed star actually grows at approximately the same rate as the underlying  background solution.  OK.  Now let me take you back to what ideas go into the well-posedness theory.  And this will bring us then to the proofs of these results.  As I mentioned, fundamental idea in this work is to use the Lagrangian coordinates.  OK.  So you've seen them.  Eta dot, OK.  Eta is simply the flow of particle trajectories.  And because I will now make an assumption of spherical symmetry, it is actually not  fundamental to these results.  But I want to make it because it will display in a very concise way the fundamental structure  that we need to worry about.  OK.  So assume everything is spherically symmetric.  Now the, maybe it's slightly pompous to say structural miracle, but it is a beautiful  structure of these equations.  If you use Lagrangian coordinates in spherical symmetry, you discover that chi, which is  now my Lagrangian map, satisfies a second order equation, chi-t-t plus f-w of chi equals  zero, where this f-w is some nonlinear operator which is of second order.  OK.  And this is an effective quasi-linear wave equation on a compact domain.  OK.  Now there's this thing here.  I don't want to display f-w because it's not beautiful.  But I will tell you later how it linearizes.  And you will be convinced that it is a second order, it's a second order self-adjoint operator  in a certain space.  Is this an artifact of the fact that this would be irrotational?  Correct.  Absolutely.  Yeah.  So I can, OK.  So now that you raised this question.  So if it were not spherically symmetric or irrotational.  If it were pure, if it were the full system.  The flow map eta would satisfy an effective wave map type equation.  The divergence of the flow map, which roughly corresponds, which captures the irrotational  part, would satisfy any equation like this.  So divergence of the jacobo, divergence of the eta t of the velocity vector field.  As the curl would satisfy, because, OK, so let me explain this.  Because remember in the Lagrangian coordinates rho dt u, where u is now pullback with respect  to the t is my Lagrangian time and u is my Lagrangian velocity.  Plus the pullback of the pressure would be something.  And your curl of this guy vanishes.  So this would give you, if you wish, a transport equation for curl and the sort of estimates  of curl of u for free.  So the idea is that when you deal with the full problem, you estimate the divergence  in the curl in two separate ways.  So you have a, you can decouple.  Absolutely.  Yeah.  So this is in some sense at the heart of the approach of Couton, Scholar, Jeng and Nassmoudi.  The framework, the functional framework is slightly different, but this is the main idea.  OK, so there is a wave equation.  That's my point.  And it's posed on a compact domain, because you pulled it back onto the, to fix the free  boundary, you had to pull it back onto something fixed.  So it's a wave equation on a compact domain.  So let me now re-derive, really, those solutions of Goldreich and Weber and other people.  I will call them homogeneous solutions, because I will make the ansatz that chi is just lambda  of t.  It doesn't depend on r.  Now this will tell you, this ansatz will reduce the previous equation to this thing  here, lambda square, lambda double dot, plus something which doesn't depend on t equals  zero.  So you can separate variables, and if you do that, you discover that there is an ODE  satisfied by lambda.  And there is another ODE in r, this time satisfied by w.  You are interested in finding compactly supported solutions, because that's what qualifies as  a star with a vacuum interface, and you can do that.  There is a magical value of delta star, such that for any delta bigger than delta star,  you can find compactly supported solutions to this equation.  And this is just an ODE.  You can really, for a range of values of delta and initial conditions, you can really classify  and see what are the rates of collapse or expansion and so on.  So I will do that for you in a second.  As I said, you can solve this ODE, prescribe the initial conditions.  There are three parameters in the problem, the delta that I mentioned, lambda 1, initial  velocity, and lambda 0, the initial radius of the star.  There is a conserved energy associated with this ODE.  It's trivial.  It's this quantity here, lambda dot squared plus 2 delta over lambda.  In the Eulerian description, so you can always go back and recover the original solution  in the Eulerian description, the density and the velocity will look like this.  So the key point, the density is like lambda to the minus 3.  So if lambda is shrinking, density is blowing up.  If lambda is expanding, the density, the fluid wants to tear itself apart.  It's decaying.  And remember, as I mentioned, self-similarity would correspond to these numbers here, 2 thirds.  So it turns out that self-similar solutions exist, but only when the energy is zero.  If energy is not zero, they're not going to obey the 2 third law.  And instead of solving the ODE, I also could not resist not putting this on the slides.  This is a kind of a bifurcation diagram that explains what's going on.  So this is a delta axis.  This is the lambda 1 axis.  And I'm fixing initial radius to be exactly 1.  So what you see here, the light blue region, this is the range of parameters that will  give you always a linearly expanding profile.  A profile that behaves asymptotically like some constant times t.  These guys here, the thick blue line, this is exactly the self-similar expansion.  And there is an explicit formula that I will give you for them on the next slide.  So you have to insist that the energy is zero.  So they're obviously unstable, because if you can go either way from the surface, the  question is, if you look at perturbations of energy zero, are they going to remain stable?  And then you have the collapsing.  So as you can see, most of the collapsing solutions in this range here exhibit self-similar  collapse.  And then there is this lower dimension line where the collapse happens at a linear rate.  So this is the structure of this space in this two-dimensional parameter space.  And the basic results are this light blue region is stable, and this thick blue region  is a co-dimension 1 stable.  So for this line, you have an explicit formula, which is just a statement that E equals zero.  Do you have a question?  You mean here?  Yes.  Well, I can tell you a very weak statement.  I can prove that these are non-linearly stable in this self-similar frame, but this is not  a good statement.  You really need to, I don't know whether, so for instance, I don't know if you perturb  an element here, whether it will collapse.  Okay, so this is an open problem.  I will state it at the end as well.  Is a linear rate of expansion the physical one, the one that has to be fixed?  So Goldreich and Weber completely discarded this for whatever reason.  And they were interested in collapse and expansion at self-similar rate.  If you look at other physics literature, they actually mostly don't want to deal with the  free boundary.  They are aware of the difficulties that come with that.  So in their literature, they think of Goldreich and Weber as the only example where you have  a sharp boundary.  In most cases, they are worried about constructing some type of solution that exhibits either  collapse or supernova or the expansion, where you have infinite support, but some sort of  decay at infinity.  This would be for them good enough.  And this also leads to good problems because most of it is also not rigorous.  So part two, how to study the stability of the expanding guys.  This is the formula that Frank, I believe, asked me about.  This is the explicit formula for self-similar expansion, whereas you don't have explicit  formulas for the linear expansion.  They are sort of approximate, but they are the leading word or just this.  So I will from now on focus only on the expanding profiles and perhaps pompously call it supernova,  but let's just say expansion.  And I want you to recall that they satisfy a wave equation of this type, and FW is some  second-order operator.  Now how to study this object?  Well, you adapt your unknown to what you believe to be stable by dividing by it.  Basically you define psi to be key, he divided by lambda bar.  And then you expect psi to be a small perturbation of one, which happens to be a steady state  of this equation in this case.  So if you make this ansatz and plug it back in here, you get this equation, second-order  PDE.  Now how to deal with this?  Well, self-similar expansion really behaves well with respect to the scaling invariances  of the problem.  So the natural thing to do is to rescale time with respect to the corresponding self-similar  rescaling.  As you will see, the linear expansion doesn't sort of, the equations don't digest it very  well.  So it will cause certain difficulties that will not be present in the self-similar case.  Okay, so let's move on to the, are there any questions about this?  So is it absolutely clear what's going on?  Okay.  So let me tell you the basic structure that will allow you to understand the stability.  So let's now focus purely on the self-similar profiles.  So remember, let me write it down.  Self-similar profile is something that expands like t to the 2 thirds.  So this is the radius of the star.  Okay, XE is a steady state of the equation on this slide.  So I will write XE as 1 plus a perturbation.  XE equal 1 is a steady state.  You do that, you demand that the initial energy is zero, that the energy of perturbation is  zero.  So this is a conserved quantity.  And of course the key idea, as I mentioned, is to pass to the self-similar time.  If you do that, so this roughly means that the new time behaves like log t.  Okay, because lambda expands like 2 thirds.  If you do that, you discover that this coefficient here is constant.  I call it B to honor some of the notation that's popular in the literature.  This quantity is a constant which is strictly less than zero.  What about the non-linearity?  So you expand the non-linearity around 1.  And you discover, and this is now the point that will perhaps clarify one of the questions,  you discover that in the linearization of this operator is precisely this second order  elliptic operator here.  But now I want you to stare at it for a second.  Notice that when r is zero, this term vanishes.  And when r is one, remember w delta is the density profile, is a density profile which  roughly looks like this.  And it also vanishes at r equal 1.  So in this sense, this operator is degenerate because it causes, it has these vanishing weights,  if you wish, at zero and at one.  This is not surprising.  This is the type of difficulty you will see in the study of just pure Euler.  But it is suggestive.  It tells you what sort of functional framework to use to, in fact, deal with this problem.  It tells you that you need to use some type of weighted spaces.  Now the cool thing about this, if you pass through the self-similar time and you express,  you now try to find the equation for phi, this is what you discover.  So phi ss, the second derivative of phi, plus a strictly positive term times phi s minus  something negative times phi, plus this linearized operator equals the right-hand side, which  is nonlinear.  Okay, and you immediately see, you realize what should be the mechanism for stability.  It is precisely this blue term here which acts as an effective damping.  This is a consequence of the fact that we are expanding our profile around an expanding  solution.  If I did it around a collapsing solution, the sign here would change and this would be sort  of the opposite of damping.  So it's a wave equation, quasi-linear wave equation on a compact domain, but with a  this damping term.  So the correct statement is that there is a high-order energy.  If it is initially sufficiently small, and if the energy, the physical energy initially  is exactly zero, then there is a global unique solution to the above equation and in fact  it decays exponentially fast in the s variable.  In the t variable, this translates into some sort of algebraic decay.  How?  The higher the energy is bound to the quesit area, I guess, right?  Precisely.  This is what you need.  This thing eats you the derivative, right?  Thank you for the question because I will now, in next slide, I will write down the  energy.  The L delta is a self-adjoint operator with respect to this weight.  Remember this vanishes at one, this vanishes at zero.  Spectral gap, okay, one is in the null space of this operator.  So if you're orthogonal to the null space, you're good.  You have a quercivity of this operator.  This per se is obvious.  What I do want to pay attention to, this first term is just the L2 norm in the first  term in the corresponding weighted space.  If you want to control the first derivative, you have to raise the level of degeneracy  in the weight W delta.  So we interpret it as first term.  Correct.  Correct.  So this is a critical observation.  No, it's only in this direction.  Absolutely.  Excuse me, the B you mentioned is, sorry, forget, is a constant?  It's a constant.  Okay.  So B is a constant due to the, because of the self-similar structure.  Because I'm expanding around self-similar solution.  So it really captures that.  This is the correct information about the self-similarity that you need to address this.  Term is obviously the good guy, this term here.  This is the bad guy.  This will give you exactly one unstable mode.  But this is not surprising.  I told you that it has to be unstable.  Where the self-similar expansion was a co-dimension one phenomenon.  So clearly this has to be the case.  But this unstable mode, which is caused by this negative term, is, well, remember, we  are starting with zero energy data.  And so the unstable direction you can check this is transversal to this zero energy surface.  So you can really control the instability, which translates in this case in controlling  the inner product of phi and the generator of the null space.  So this allows you to carry on the estimates.  So think of this as, so this is the relationship that you get by linearizing energy equals  zero.  This is the quantity that you will control just by doing some simple energy identity.  This is the quantity that you want to control to get a spectral gap.  Now answer to Pierre's question.  Indeed, the energy is motivated in this case quite concretely, in the second case more  philosophically, in the case of linear expansion, by the works of these people.  And now we come to your observation.  Indeed, I don't want you to gain any intuition from this.  All I want you to see is that each time I take a further normal derivative, I have to  raise some index in my weight.  And this is how this works.  These delta k spaces basically mean you add k, you add a k-fold power of w delta to your  weight.  And this is the only way you can do it.  So more spatial derivatives implies more degenerate weights.  And now you will trust me if I tell you how the energy method works.  You try to prove something like this.  This term here kind of comes from that damping term that I mentioned.  Of course, you need to control the energy by this dissipative term here.  You can do that because of the spectral gap property.  It's a technical thing, but you can do it.  And this, in principle, gives you the exponential decay and completes the proof.  But perhaps I want to say, going back one slide, this is a mixed space time norm.  In reality, you do the estimates only with the time vector field because it commutes  with the equation entirely.  And then you use elliptic estimate to control the higher space derivatives with the time  derivatives.  So this is what I just said.  Use the time derivatives to build an energy, elliptic estimates.  Now what are the technical tools that go into closing the estimates?  Remember there is a nonlinear right hand side which actually looks quite dirty when you  write it down.  So the tools are hardy inequalities which exactly allow you to deal with those weights  that degenerate at 0 and 1 and embeddings between weighted sub-allows spaces.  If you take enough derivatives, you can control lower norm by the energy, L infinity norms  of lower order terms by the energy.  So I will not go into any of this.  This is technical.  I have eight minutes left.  I will, OK, so let me then say just a few words about this second problem here and maybe  try to hint at least at why is it more difficult in the first one.  The difficulty appears not to be conceptual but more technical.  But it is there and it took us some time to understand.  So the reason why this problem is harder conceptually is because you cannot, it doesn't honor the  self-similar structure of the problem.  Doesn't honor the scaling in variances.  These solutions are there but how do you prove that they are stable?  So it will necessitate the use of a, so first of all, the fact that it is not honoring the  invariance will disallow us from using time vector fields.  As you will see, they will not commute well with the problem because now when I write down  the equation for my perturbation phi, it will have time dependent coefficients.  OK, so we have to come up with a correct spatial derivative that will capture the degeneracies  at the vacuum.  And this is precisely this vector field here, which just happens to be a five dimensional  aplosion but this is a coincidence.  OK, so what we do, we rescale time so that tau, the new time becomes log t again.  So this means dividing by this linearly expanding profile.  Let me write it down.  OK, here.  So lambda tilde of t, think of it as something that behaves like a constant times t.  And if you write down the equation for the perturbation, now you discover that you have  these tau dependent factors in front and lambda tilde is growing exponentially in this new  variable.  OK, so you have e to the, let's say e to the tau plus e to the tau plus this stuff equals  right inside.  This is still causing a certain damping effect.  It is well known that the expansion produces a sort of a stabilizing effect.  It's reminiscent of works, for instance, in a special relativity, in the general relativity  when you have cosmological constant, for instance.  But, and this is the statement, there is a high order energy e tilde that will allow  you that will give you global existence for small data.  However, you cannot classify the asymptotic attractor.  So all you can say is that the time derivative, the tau derivative of the solution decays,  but you cannot say that in the limit you belong to a nearby member of the linearly expanding  family.  This energy is constructed by using powers of this elliptic operator that I just wrote  a second ago.  And so to be honest, I can talk about it.  I have a few more slides at the end, but let me not go into the technical aspects of this  proof.  Let me just say that this operator is designed to capture very precisely.  So in other words, you cannot just blindly apply dr derivative or dx.  If you do that, you will create singularities at zero and they are hard to handle.  So you have to do it with a particular combination.  And this operator s is a reflection of the structure of that weighted-generate operator  L delta.  And let's leave it at that.  So I will not go into the technical details.  Yeah, that's a good question.  So it's just something that sits nearby, but it's not one of the...  But you have to identify the proof, right?  Yeah, yeah.  Well, you can prove that the time derivative of that is decaying.  So that gives you convergence to something.  Yeah, yeah.  And it's nearby.  Yeah.  So the reminiscent of a result I heard recently in a talk by Falker, Falker Schuhe, and sort  of in these expanding space times.  Who's here, I believe?  Yeah, very good.  Let me spend the remaining three minutes discussing perhaps questions that are...  I find very interesting that are open.  So here's this kind of a silly lemma.  I call it lemma.  It's too ambitious to call it lemma.  Do there exist collapsing or expanding compactly supported stars?  So this is the key for other values of gamma.  They should exist in the supercritical range, which is between 6 fifths and 4 thirds.  This is an open question.  What you can show, however, is that self-similarity and a spherical symmetry cannot coexist for  this type of collapse.  Lama being 4 thirds is absolutely critical for the two to be able to combine the two.  And this is a very silly thing to show.  It's just some conservation law type calculation.  So the question is, can you construct exisymmetric self-similar collapse in this mass supercritical  range?  I think this is a very interesting problem.  And now we come to the question that I believe Sergio asked.  Are the Lane-Emdon stars that are mentioned before in the subcritical range, where they  are shown to be conditionally stable, are they in fact stable?  So a possible mechanism, and you will see this discussed sometimes in the physics literature,  is some type of shock formation.  This is not known, and this is really important.  Let me remind you of this picture.  The pink elephant in the room is, of course, this region here, the collapsing region.  So as I said, it is not difficult to show that it's nonlinearly unstable.  The moment you have correct well-posedness theory and you have a growing mode, it's very  easy.  But so again, this should be a lemma.  But it's really unsatisfactory, because it doesn't tell you anything about the real instability.  In particular, it could it be that the nearby guys decay, collapse at perhaps a corrected  rate?  Do they collapse at all?  So this is just in spherical symmetry.  You don't have to go outside spherical symmetry to ask these questions.  And finally, I would like to sort of point out this viewpoint that I believe studying  this family of problems is really the right thing, because you really want to have a robust  understanding of how equations of state affect the collapse.  This is partly motivated by general relativity.  If you think of the most famous example of a collapsing star in general relativity, it's  the Oppenheimer-Sneider solution, which corresponds to the case where you don't have any pressure.  So by a result of Christodulu from 83, I believe, you can show that these are generically unstable  solutions.  If you go outside the homogeneity class, you generically form naked singularities.  So you want to, well, you can devise various conclusions out of that.  One conclusion that I like to draw out of that is that studying non-trivial equations  of state is an important problem, and this is just the Newtonian version of the problem.  But realistic equations of state are, this is just a final comment, is in principle quite  a hard problem, even for physicists.  So this is my favorite relativistic astrophysics book by Zeldovich and Noviko, where I sort  of learned this stuff about.  And they discuss one of the sort of high density regimes for the equations of state.  And here's what Zeldovich says.  So it reminds him of Avocenko's parody.  So it is precisely like this.  So if you look at the astrophysical literature, there's a huge disparity between what people  think is the right equation of state.  Yeah.  So this is a teaser, I think of it as a teaser trailer for this wonderful textbook because  it's a very serious book, but as I browse through it every now and then I find a gem  of this type.  Thank you.  So the question is, how about?  So it seems, since the non-year expanded rate is unstable, does it mean that basically  the type of solutions which would converge to the stable one, which is a new one?  Okay, so if you generically perturb a self-similar expanding rate, if you generically perturb  it with perturbations that have positive energy, then they will converge to the linear guy.  So it's generically unstable and these guys are stable really only in the co-dimension  one sense.  So the linear is stable.  The linear is stable.  The linear takes over, not the self-similar.  Yeah.  Yeah, yeah, exactly.  So do the physicists know this part?  I mean, I'm not...  It should be of great interest to them, I guess.  Well, if you're...  If you're very well, I should be linear.  Yeah, exactly.  So if you read Goldreich and Weber, which was one of the starting points for this project,  for some reason they really say, we will not be interested in these linear rates.  They drop a constant of integration in a second order ODE, which allows them to go down one  level and get an explicit solution, which is this.  And they literally say, ah, we kind of don't want to look at that.  So but this is the finding.  I mean, part of the mess was also to understand the exact structure of this parameter space.  Yes.  The non-linear unstable result covers the whole regime of gamma written there or...?  Actually it actually goes up to here.  I did not tell you what happens at 6-fifths.  At 6-fifths you still get a steady state, but it is supported all the way up to...  It's infinite.  It has an infinite support.  And it turns out to be...  So in this case, the steady state is exactly the steady state of the critical Quintic wave  equation.  So it's one of these guys.  I don't know.  Once there's some constant here, and then I believe one over two, something like that.  So this is the rate of decay of this guy.  So it is the most complicated case to prove the most important?  In this work, well, it's quite different because here you are working on a compactly  supported domain.  So you really have to rely on the well-posed in this theory for vacuum interface problems.  Here it's different.  Here the weights are at infinity.  So it has a flavor of it, but it's quite different.  I'm not sure how to...  So this result is also due to Jang.  So this equality case, all of this is due to Jang.  Sorry, in this range.  I have to be careful.  This is the instability range.  What this result shows in a way is that the fourth-third case, which is linearly stable,  is in fact widely unstable.  So your result is radial from the stability.  This is entirely radial.  So the paper that will show up is entirely radial, but it holds also for non-radial perturbations.  So I alluded earlier to that.  Here's the reason why the radial case is simpler, not just for obvious reasons, but also because  it reduces the whole problem to just a quasi-linear wave equation.  When you work with a non-radial case, so the fluid is not irrotational anymore in general,  so you have to account for the vorticity, which does not have this wave-like structure,  but it satisfies a certain transport equation that allows you to control it.  What makes the vorticity better than in the standard equation where you know that you  cannot control the vorticity?  Is it because of the?  Well, it's a local in-time result.  I mean, yeah.  No, but is it draw?  Yeah, yeah.  So all I'm saying, I'm not sure if I understand your question, but all I'm saying is that  Ah, the long time.  For the long time result, you don't have a long time.  Oh, we haven't produced it yet.  No, no.  But for the irrotational case, it will be straightforward.  And for the vorticity in the long time result, I don't want to make the claim yet.  I don't think- Because I don't see why it should be any better than in the regular case.  Maybe.  No, no, no.  You have this damping effect, so it will be better.  It will be better.  But will it affect the vorticity?  No.  You have, for instance, for Friedman, Lomet, Roberts and Walker space times, you have  a stability without irrotational assumption.  So when you have a cosmological constant, for instance, and you perturb away from that,  you see this.  You can control the full flow.  Of course, this is not a free boundary problem, because this free boundary problem in GR is  more- Okay, but the mechanism is massive.  Yeah, yeah.  You will kill it with that.  All right.  So thank you.  So we start again at 11.30.  So let's go to the next one.  Thank you.  Thank you.  Thank you.  Thank you.  Thank you.  Thank you.  Thank you.  Thank you.
Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_0_1223.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_1224_1559.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_1560_1679.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_1680_2075.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_2076_2183.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_2184_2627.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_2628_3263.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_3264_3539.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_3540_5087.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_5088_6323.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_6324_9743.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_9744_9899.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_9900_10607.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_10608_10763.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_10764_12623.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_12624_12671.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_12672_13187.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_13188_13247.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_13248_16631.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_16632_16679.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_16680_18491.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_18492_18863.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_18864_19655.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_19656_20591.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_20592_21143.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_21144_21191.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_21192_21479.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_21480_26123.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_26124_26471.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_26472_26951.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_26952_27875.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_27876_28919.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_28920_29051.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_29052_29507.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_29508_30107.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_30108_30503.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_30504_33563.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_33564_35471.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_35472_35819.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_35820_37259.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_37260_37907.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_37908_42551.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_42552_44855.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_44856_44903.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_44904_45479.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_45480_48887.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_48888_49127.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_49128_49235.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_49236_49847.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_49848_54611.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_54612_54683.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_54684_55019.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_55020_55847.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_55848_57635.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_57636_57719.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_57720_62891.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_62892_62963.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_62964_63371.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_63372_64103.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_64104_65327.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_65328_68447.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_68448_68555.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_68556_71111.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_71112_72743.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_72744_73727.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_73728_74351.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_74352_74471.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_74472_77375.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_77376_78335.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_78336_78467.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_78468_80687.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_80688_80747.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_80748_81287.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_81288_81599.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_81600_81923.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_81924_82019.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_82020_82067.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_82068_82847.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_82848_83663.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_83664_83735.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_83736_84287.png Mahir_HADZIC_-_Nonlinear_stability_of_expanding_stars_in_the_mass-critical_Euler-Poisson_system_84288_85506.png
I would first like to thank the organizers for this wonderful conference. It's the first time I'm here. So before I say anything about math, everything I will talk about is a joint work with Juhi Jang who is based at the University of Southern California and this year at the IAS. So the good thing about this long title is that it basically tells you the statement of the main theorem. There are certain special solutions to the system of PDEs I will look at. Euler Poisson system, it comes from astrophysics, it's a model of a star. And we will be studying the nonlinear stability of that class of solutions. So let me give you an outline. Well first, I intend to spend quite a bit of time introducing the system and talking about general properties, perhaps at the expense of giving you some technical details in the second half of the talk. But I can promise that there will be two different stability results pertaining to two different types of special solutions that we will look at. And one of them, the first one is easier. This is the one where I will probably give you a pretty complete idea how the proof works and hopefully there will be time to talk about the second type of expansion. So before any of that, let me just tell you that the gravitational Euler Poisson system is the model you would find in a typical astrophysics book describing a star. So it is a Newtonian model. The idea is that the star is a gas or a fluid kept together by gravitational forces. So let me tell you what this talk is not about. It's not about relativistic models of stars, but nevertheless part of the motivation for studying this model indeed comes from relativity. And I will mention that hopefully at the end of the talk. So as I said, the star is idealized as a compact fluid body. And the word compact is crucial because it suggests that there is a boundary between the fluid and the vacuum. And if you want to model such a problem and if you want to track the behavior of the boundary, you have to allow for that boundary to move. So this causes, this gives us a moving star vacuum interface. So what is the model? Well, first I'll tell you the unknowns. The first three are unsurprising and you're used to them. The fluid density rho and the pressure P, the velocity field U, and the gravitational field phi that you need to couple to the Euler equations. The fourth unknown or the further unknown is the support of the fluid because we are trying to treat a free boundary problem. Okay, and the equations, again, you have seen some fluid models this time around. Don't be scared. We will soon forget about this slide. The first equation is the continuity equation. The second one expresses the conservation of momentum. And on the right hand side, you will see a term caused by the gravitational field phi. And the gravitational field phi is self-consistent. It's generated by the mass density here. So Laplace phi is 4 pi rho. This is the Poisson equation and these are the equations of compressible fluid dynamics. On top of that, we have to supply this problem with some boundary conditions and remember this boundary is moving. So the kinematic boundary condition is the natural statement that the boundary is moved by the fluid particles. So the normal velocity of the boundary is, in fact, the normal component of the velocity field at the free boundary. And the pressure vanishes at the boundary. Okay? Now, if you count the unknowns, you will notice that there is one superfluous unknown here, namely the pressure. There is no evolution equation for pressure and you need to prescribe a further equation to or a relation to actually close this system. There is one equation for rho, three for u, one for phi, but there is nothing for pressure. And this is typically done by prescribing the so-called equation of state. Now equations of state is something that in the physics literature, you would look at a particular model of a star and you would hope to have a good understanding of what the equation of state is for that particular star. What is commonly done, people study the so-called polytropic equation of state where p of rho is rho to a certain power gamma. And this gamma is the so-called adiabatic index and it's allowed to vary between one and two. So the way I sort of would like you to think about this model, I would like to think of this nonlinear problem as parameterized by gamma. Okay? And as we vary gamma, the qualities, the qualitative properties of the solutions of this system will change. Okay. So let me comment a little bit on the literature. The vacuum interface itself is known to cause severe analytical difficulties in improving local, even local well-positiveness. So the basic result in the theory of compressible fluids with vacuum interfaces are basically these two results, they came independently in 2010 by Couton and Scholar and Jang and Masmoudi. And so what they have to contend with is a certain degeneracy. I will explain it as we go through the talk that is caused by the presence of the vacuum. But one of the key insights is the use of the Lagrangian coordinates. So if you use Lagrangian coordinates and if you rephrase the compressible Euler equations surrounded by vacuum, you will discover that there is a wave equation, a quasi-linear wave equation lurking in the background. Okay? So that's somehow the one of the big revelations. It is a degenerate wave equation, degenerate in the sense that I will specify later, but the word degenerate has to do with the fact that some of the coefficients of the wave equation will vanish at the vacuum boundary. And another key assumption, so another key assumption in these works is that of a physical vacuum boundary condition. So here I'm stating what it is. It's a statement that the fluid enthalpy, which is this quantity, rho to the power gamma minus 1, has to have its normal derivative at the vacuum boundary, has to have a sign. This is perhaps not a well-motivated condition. So physically it expresses the fact that particles have to accelerate in the normal direction at the boundary. But I will soon actually motivate this, sorry. I will soon motivate this condition by finding particular steady-state solutions of the Euler-Pollisson system that satisfy this. And indeed, if you open an astrophysics textbook, the steady-state star models that they look at all satisfy this type of condition, the compactly supported ones. Okay, so there's a little bit of a discrepancy here. Here I'm citing works that pertain to a compressible Euler equation. You could ask me, what about Euler-Pollisson? Because I'm adding a gravitational field. From the point of view of well-pollinist theory, this is a lower-order perturbation. So if you understand how to prove well-pollinist here, you will understand how to prove well-pollinist for Euler-Pollisson system. But nevertheless, there is a reference that I'm missing. Where is the chalk? Oh, fantastic. Thank you. So prior to this work, there's one result which deals only with a priori estimates with the Lienblad, Couton-Lienblad and Scholar. So this is one missing. So an a priori, Bounds, Couton-Lienblad, Scholar. And for the Euler-Pollisson system, this came much later, entirely based on these ideas. There's a work of Lu-Shin and Zhang. Okay, this is for you. Not for the vacuum interface. So for the liquid case, when the density is not vanishing at the boundary. So this is a significantly less degenerate model. Okay, so instead of, I will now try to give you, to do three things for you. I will give you an example of this famous class of steady states known as Lien-Emden-Stars. Then I will, I will, I'm expecting a sigh of relief. I will produce an invariant scale, rescaling for this problem and explain a certain notion of criticality that enters and that will motivate the main results. So most famous class of special solutions are the steady states. And the way you produce them, you set rho to be only a function of x and u to be zero. You plug it in and the whole equation reduces to this statement here, gradient of rho star gamma plus rho star gradient phi star equals zero. So if you use the fact that this satisfies the Poisson equation and if you call this object here, this is the fluid enthalpy, if you call it w, then in spherical symmetry, and this is a very classical stuff, you will see this in astrophysics textbooks, you reduce the problem to solving this ordinary differential equation. And the key question is, for what powers of gamma can I find a compactly supported solution with finite mass? This is your star. So the answer is, if gamma lies in this magical range, then you can find profiles that are indeed compactly supported, say if it supports zero, one, and they look exactly like that. Once gamma hits six-fifths and below, no such solutions exist. It's just the ODE analysis. And moreover, these special solutions in fact satisfy the physical vacuum conditions. So w prime at one is strictly less than zero. So the density profile is not smooth across the vacuum interface. And this is one of the critical points that works of Couton, Scholar, Jang and Masmoudi had to address, the fact that the non-smoothness is a generic part of the theory. So you have to find a way, what does it even mean to differentiate the equation in a proper sense. You can ask however, and this is what physicists did, are these solutions stable? So now the numerology gets even more interesting. If gamma lies in this range here, the lane end and stars are linearly unstable. And if gamma is bigger, equal than four-thirds, it's in fact linearly stable. So it's just some basic linear analysis. Of course it begs a question about two questions. Are these results, can they be upgraded to a nonlinear statement? So the first result can, it is not terribly surprising. Once you have a growing mode, which is indeed the main source of instability in this regime, and once you have a good well-posedness theory, then you can prove nonlinear instability. And this was shown by Jang in 2014. Nonlinear stability is a much more contentious point. So there is a conditional nonlinear stability result in a fairly weak topology generated by basically just the energy by Ryan from 2003, which says if the solution exists, it will be nonlinearly stable. However, there is absolutely no reason why the solution should exist, why small perturbations of these steady states should. And it's in the same regime, four-thirds to two. Pardon? Four? The result of Ryan, is it in that regime? Right, so the result of Ryan, the topology is so weak that it doesn't really see the free boundary. In particular, this type of stability result, if you think about it, will actually allow you to have support splitting up and things like that, the topology will not see it. So you really need something else. And in fact, I will talk about it at the very end. I will highlight this as one of the important open problems in this development. Okay, so let me try to perhaps explain these exponents that show up in a different way. And this hopefully will help you appreciate them. There is an invariant rescaling of this problem. To each gamma, you can associate this invariant rescaling where you see where you just work it out and you realize that if rho and u are solutions, then rho tilde and u tilde are also solutions. And the way you rescale the time depends on gamma. So really the value of gamma has a lot to say about the structure of the solutions of the equation. So similar to the pressure and the potential, you can work it out from these how to rescale. Okay, the second ingredient that I'm sure you will appreciate, what are the conserved quantities? Well, the mass, this is easy, the integral of the density, and the total energy. Let's just briefly talk about the total energy. This bit here is the kinetic energy. This object here, remember, rho is Laplace phi. So if you integrate by parts, this becomes negative gradient phi squared. So this is just the gravitational energy, which gives a negative contribution. And this thing here, rho to the power gamma, is some sort of thermal energy of the star. Okay, and the question that is quite natural in this context is what is the behavior under the rescaling? And then you discover that the mass remains invariant exactly when gamma equals 4 thirds and energy remains invariant exactly when gamma equals 6 fifths. And these are the two numbers that are popped up here. So this is to suggest that these numbers, this numerology has a special meaning. All of the results for the remainder of the article of this presentation will pertain to the mass critical case gamma equals 4 thirds. I said article because the actual article is supposed to appear tomorrow on the pre-print server on the archive. And there are several spots where we said and erased all the results for the remainder of this article pertain to this case. Anyway, but okay, this is the key point to remember. So gamma equals 4 thirds is the mass critical case. I'll tell you why physicists like this particular value of gamma. You see, it's a rescaling that keeps the mass preserved. So it's a very appealing to think of gravitational collapse or expansion as some process that sort of cascades through these scales, lambda, but it preserves the mass. It's very natural. So the beautiful thing about this is that in fact there is a result of Goldreich and Weber who were physicists in 1980 where they found explicit examples of collapsing and expanding solutions of examples of stars in this mass critical case. And now let me do some computation. So if gamma equals 4 thirds and if my self-similar rescaling is this, so you can work out what this coefficient is here. So 2 minus 4 thirds is 2 thirds. One over that is 3 halves. So this becomes P over lambda to the power 3 halves. So to have a self-similar collapse or a self-similar rescaling would correspond to having a collapsing solution whose radius scales like T minus T to the power 2 thirds or an expanding solutions that expands at the same rate. This would be the corresponding coefficient if you wanted to call this a self-similar collapse. However, if you follow through the paper of Goldreich and Weber, already implicitly there and then explicitly in the works of Makino and Lin, you will find a second family of expanding and collapsing solutions that expand and collapse at the linear rate. And this was part of the mystery for us to understand, first of all, the reasons why this is so. Yes? Sorry. Which one corresponds to the explicit solution? This one. This is not explicit. This is just, you see, I'm putting this sign here. So this is informal. In a few slides, I will give you the exact formula. There is a correction term. For this guy, there is no correction term. This is perfect. This is corrected. By your log? No. No, no, no. It's just something linear. You will see. You will see. No logs. Well, we'll come to that at the end of the talk as well. So there are two families of such a family? There are two families. No, no, no. It's not very self-similar. I mean this. But there are two families of collapsing and expanding solutions. So for whatever reasons, in these works, you will never read this off. They will never talk about the rates. They analyze this and they say there exist these solutions. But if you really look at the rates, this is what will happen. I will give a more precise statement later on. So this is an informal statement. Informal statement. Don't ask me for any norms and so on at this moment. This will also come later. The basic result is that these self-similar expanding profiles, so I'm not talking about collapse, the ones that behave like this, like 2 to the power 2 thirds, are co-dimension one stable. And the linearly expanding guys will be just stable. In this case, we can actually talk about asymptotic stability in this co-dimension one set of perturbations. Whereas here, we do not characterize the asymptotic attractors. There are just some objects that do not necessarily belong to some sort of nearby linearly expanding profile. OK. This is the gist of the results. I will get into the technical statement slightly later. In particular, these are, to my knowledge at least, first non-trivial examples of global solutions for this free boundary problem. And what is very important for us, we can actually say that the support of the star, of the nearby perturbed star actually grows at approximately the same rate as the underlying background solution. OK. Now let me take you back to what ideas go into the well-posedness theory. And this will bring us then to the proofs of these results. As I mentioned, fundamental idea in this work is to use the Lagrangian coordinates. OK. So you've seen them. Eta dot, OK. Eta is simply the flow of particle trajectories. And because I will now make an assumption of spherical symmetry, it is actually not fundamental to these results. But I want to make it because it will display in a very concise way the fundamental structure that we need to worry about. OK. So assume everything is spherically symmetric. Now the, maybe it's slightly pompous to say structural miracle, but it is a beautiful structure of these equations. If you use Lagrangian coordinates in spherical symmetry, you discover that chi, which is now my Lagrangian map, satisfies a second order equation, chi-t-t plus f-w of chi equals zero, where this f-w is some nonlinear operator which is of second order. OK. And this is an effective quasi-linear wave equation on a compact domain. OK. Now there's this thing here. I don't want to display f-w because it's not beautiful. But I will tell you later how it linearizes. And you will be convinced that it is a second order, it's a second order self-adjoint operator in a certain space. Is this an artifact of the fact that this would be irrotational? Correct. Absolutely. Yeah. So I can, OK. So now that you raised this question. So if it were not spherically symmetric or irrotational. If it were pure, if it were the full system. The flow map eta would satisfy an effective wave map type equation. The divergence of the flow map, which roughly corresponds, which captures the irrotational part, would satisfy any equation like this. So divergence of the jacobo, divergence of the eta t of the velocity vector field. As the curl would satisfy, because, OK, so let me explain this. Because remember in the Lagrangian coordinates rho dt u, where u is now pullback with respect to the t is my Lagrangian time and u is my Lagrangian velocity. Plus the pullback of the pressure would be something. And your curl of this guy vanishes. So this would give you, if you wish, a transport equation for curl and the sort of estimates of curl of u for free. So the idea is that when you deal with the full problem, you estimate the divergence in the curl in two separate ways. So you have a, you can decouple. Absolutely. Yeah. So this is in some sense at the heart of the approach of Couton, Scholar, Jeng and Nassmoudi. The framework, the functional framework is slightly different, but this is the main idea. OK, so there is a wave equation. That's my point. And it's posed on a compact domain, because you pulled it back onto the, to fix the free boundary, you had to pull it back onto something fixed. So it's a wave equation on a compact domain. So let me now re-derive, really, those solutions of Goldreich and Weber and other people. I will call them homogeneous solutions, because I will make the ansatz that chi is just lambda of t. It doesn't depend on r. Now this will tell you, this ansatz will reduce the previous equation to this thing here, lambda square, lambda double dot, plus something which doesn't depend on t equals zero. So you can separate variables, and if you do that, you discover that there is an ODE satisfied by lambda. And there is another ODE in r, this time satisfied by w. You are interested in finding compactly supported solutions, because that's what qualifies as a star with a vacuum interface, and you can do that. There is a magical value of delta star, such that for any delta bigger than delta star, you can find compactly supported solutions to this equation. And this is just an ODE. You can really, for a range of values of delta and initial conditions, you can really classify and see what are the rates of collapse or expansion and so on. So I will do that for you in a second. As I said, you can solve this ODE, prescribe the initial conditions. There are three parameters in the problem, the delta that I mentioned, lambda 1, initial velocity, and lambda 0, the initial radius of the star. There is a conserved energy associated with this ODE. It's trivial. It's this quantity here, lambda dot squared plus 2 delta over lambda. In the Eulerian description, so you can always go back and recover the original solution in the Eulerian description, the density and the velocity will look like this. So the key point, the density is like lambda to the minus 3. So if lambda is shrinking, density is blowing up. If lambda is expanding, the density, the fluid wants to tear itself apart. It's decaying. And remember, as I mentioned, self-similarity would correspond to these numbers here, 2 thirds. So it turns out that self-similar solutions exist, but only when the energy is zero. If energy is not zero, they're not going to obey the 2 third law. And instead of solving the ODE, I also could not resist not putting this on the slides. This is a kind of a bifurcation diagram that explains what's going on. So this is a delta axis. This is the lambda 1 axis. And I'm fixing initial radius to be exactly 1. So what you see here, the light blue region, this is the range of parameters that will give you always a linearly expanding profile. A profile that behaves asymptotically like some constant times t. These guys here, the thick blue line, this is exactly the self-similar expansion. And there is an explicit formula that I will give you for them on the next slide. So you have to insist that the energy is zero. So they're obviously unstable, because if you can go either way from the surface, the question is, if you look at perturbations of energy zero, are they going to remain stable? And then you have the collapsing. So as you can see, most of the collapsing solutions in this range here exhibit self-similar collapse. And then there is this lower dimension line where the collapse happens at a linear rate. So this is the structure of this space in this two-dimensional parameter space. And the basic results are this light blue region is stable, and this thick blue region is a co-dimension 1 stable. So for this line, you have an explicit formula, which is just a statement that E equals zero. Do you have a question? You mean here? Yes. Well, I can tell you a very weak statement. I can prove that these are non-linearly stable in this self-similar frame, but this is not a good statement. You really need to, I don't know whether, so for instance, I don't know if you perturb an element here, whether it will collapse. Okay, so this is an open problem. I will state it at the end as well. Is a linear rate of expansion the physical one, the one that has to be fixed? So Goldreich and Weber completely discarded this for whatever reason. And they were interested in collapse and expansion at self-similar rate. If you look at other physics literature, they actually mostly don't want to deal with the free boundary. They are aware of the difficulties that come with that. So in their literature, they think of Goldreich and Weber as the only example where you have a sharp boundary. In most cases, they are worried about constructing some type of solution that exhibits either collapse or supernova or the expansion, where you have infinite support, but some sort of decay at infinity. This would be for them good enough. And this also leads to good problems because most of it is also not rigorous. So part two, how to study the stability of the expanding guys. This is the formula that Frank, I believe, asked me about. This is the explicit formula for self-similar expansion, whereas you don't have explicit formulas for the linear expansion. They are sort of approximate, but they are the leading word or just this. So I will from now on focus only on the expanding profiles and perhaps pompously call it supernova, but let's just say expansion. And I want you to recall that they satisfy a wave equation of this type, and FW is some second-order operator. Now how to study this object? Well, you adapt your unknown to what you believe to be stable by dividing by it. Basically you define psi to be key, he divided by lambda bar. And then you expect psi to be a small perturbation of one, which happens to be a steady state of this equation in this case. So if you make this ansatz and plug it back in here, you get this equation, second-order PDE. Now how to deal with this? Well, self-similar expansion really behaves well with respect to the scaling invariances of the problem. So the natural thing to do is to rescale time with respect to the corresponding self-similar rescaling. As you will see, the linear expansion doesn't sort of, the equations don't digest it very well. So it will cause certain difficulties that will not be present in the self-similar case. Okay, so let's move on to the, are there any questions about this? So is it absolutely clear what's going on? Okay. So let me tell you the basic structure that will allow you to understand the stability. So let's now focus purely on the self-similar profiles. So remember, let me write it down. Self-similar profile is something that expands like t to the 2 thirds. So this is the radius of the star. Okay, XE is a steady state of the equation on this slide. So I will write XE as 1 plus a perturbation. XE equal 1 is a steady state. You do that, you demand that the initial energy is zero, that the energy of perturbation is zero. So this is a conserved quantity. And of course the key idea, as I mentioned, is to pass to the self-similar time. If you do that, so this roughly means that the new time behaves like log t. Okay, because lambda expands like 2 thirds. If you do that, you discover that this coefficient here is constant. I call it B to honor some of the notation that's popular in the literature. This quantity is a constant which is strictly less than zero. What about the non-linearity? So you expand the non-linearity around 1. And you discover, and this is now the point that will perhaps clarify one of the questions, you discover that in the linearization of this operator is precisely this second order elliptic operator here. But now I want you to stare at it for a second. Notice that when r is zero, this term vanishes. And when r is one, remember w delta is the density profile, is a density profile which roughly looks like this. And it also vanishes at r equal 1. So in this sense, this operator is degenerate because it causes, it has these vanishing weights, if you wish, at zero and at one. This is not surprising. This is the type of difficulty you will see in the study of just pure Euler. But it is suggestive. It tells you what sort of functional framework to use to, in fact, deal with this problem. It tells you that you need to use some type of weighted spaces. Now the cool thing about this, if you pass through the self-similar time and you express, you now try to find the equation for phi, this is what you discover. So phi ss, the second derivative of phi, plus a strictly positive term times phi s minus something negative times phi, plus this linearized operator equals the right-hand side, which is nonlinear. Okay, and you immediately see, you realize what should be the mechanism for stability. It is precisely this blue term here which acts as an effective damping. This is a consequence of the fact that we are expanding our profile around an expanding solution. If I did it around a collapsing solution, the sign here would change and this would be sort of the opposite of damping. So it's a wave equation, quasi-linear wave equation on a compact domain, but with a this damping term. So the correct statement is that there is a high-order energy. If it is initially sufficiently small, and if the energy, the physical energy initially is exactly zero, then there is a global unique solution to the above equation and in fact it decays exponentially fast in the s variable. In the t variable, this translates into some sort of algebraic decay. How? The higher the energy is bound to the quesit area, I guess, right? Precisely. This is what you need. This thing eats you the derivative, right? Thank you for the question because I will now, in next slide, I will write down the energy. The L delta is a self-adjoint operator with respect to this weight. Remember this vanishes at one, this vanishes at zero. Spectral gap, okay, one is in the null space of this operator. So if you're orthogonal to the null space, you're good. You have a quercivity of this operator. This per se is obvious. What I do want to pay attention to, this first term is just the L2 norm in the first term in the corresponding weighted space. If you want to control the first derivative, you have to raise the level of degeneracy in the weight W delta. So we interpret it as first term. Correct. Correct. So this is a critical observation. No, it's only in this direction. Absolutely. Excuse me, the B you mentioned is, sorry, forget, is a constant? It's a constant. Okay. So B is a constant due to the, because of the self-similar structure. Because I'm expanding around self-similar solution. So it really captures that. This is the correct information about the self-similarity that you need to address this. Term is obviously the good guy, this term here. This is the bad guy. This will give you exactly one unstable mode. But this is not surprising. I told you that it has to be unstable. Where the self-similar expansion was a co-dimension one phenomenon. So clearly this has to be the case. But this unstable mode, which is caused by this negative term, is, well, remember, we are starting with zero energy data. And so the unstable direction you can check this is transversal to this zero energy surface. So you can really control the instability, which translates in this case in controlling the inner product of phi and the generator of the null space. So this allows you to carry on the estimates. So think of this as, so this is the relationship that you get by linearizing energy equals zero. This is the quantity that you will control just by doing some simple energy identity. This is the quantity that you want to control to get a spectral gap. Now answer to Pierre's question. Indeed, the energy is motivated in this case quite concretely, in the second case more philosophically, in the case of linear expansion, by the works of these people. And now we come to your observation. Indeed, I don't want you to gain any intuition from this. All I want you to see is that each time I take a further normal derivative, I have to raise some index in my weight. And this is how this works. These delta k spaces basically mean you add k, you add a k-fold power of w delta to your weight. And this is the only way you can do it. So more spatial derivatives implies more degenerate weights. And now you will trust me if I tell you how the energy method works. You try to prove something like this. This term here kind of comes from that damping term that I mentioned. Of course, you need to control the energy by this dissipative term here. You can do that because of the spectral gap property. It's a technical thing, but you can do it. And this, in principle, gives you the exponential decay and completes the proof. But perhaps I want to say, going back one slide, this is a mixed space time norm. In reality, you do the estimates only with the time vector field because it commutes with the equation entirely. And then you use elliptic estimate to control the higher space derivatives with the time derivatives. So this is what I just said. Use the time derivatives to build an energy, elliptic estimates. Now what are the technical tools that go into closing the estimates? Remember there is a nonlinear right hand side which actually looks quite dirty when you write it down. So the tools are hardy inequalities which exactly allow you to deal with those weights that degenerate at 0 and 1 and embeddings between weighted sub-allows spaces. If you take enough derivatives, you can control lower norm by the energy, L infinity norms of lower order terms by the energy. So I will not go into any of this. This is technical. I have eight minutes left. I will, OK, so let me then say just a few words about this second problem here and maybe try to hint at least at why is it more difficult in the first one. The difficulty appears not to be conceptual but more technical. But it is there and it took us some time to understand. So the reason why this problem is harder conceptually is because you cannot, it doesn't honor the self-similar structure of the problem. Doesn't honor the scaling in variances. These solutions are there but how do you prove that they are stable? So it will necessitate the use of a, so first of all, the fact that it is not honoring the invariance will disallow us from using time vector fields. As you will see, they will not commute well with the problem because now when I write down the equation for my perturbation phi, it will have time dependent coefficients. OK, so we have to come up with a correct spatial derivative that will capture the degeneracies at the vacuum. And this is precisely this vector field here, which just happens to be a five dimensional aplosion but this is a coincidence. OK, so what we do, we rescale time so that tau, the new time becomes log t again. So this means dividing by this linearly expanding profile. Let me write it down. OK, here. So lambda tilde of t, think of it as something that behaves like a constant times t. And if you write down the equation for the perturbation, now you discover that you have these tau dependent factors in front and lambda tilde is growing exponentially in this new variable. OK, so you have e to the, let's say e to the tau plus e to the tau plus this stuff equals right inside. This is still causing a certain damping effect. It is well known that the expansion produces a sort of a stabilizing effect. It's reminiscent of works, for instance, in a special relativity, in the general relativity when you have cosmological constant, for instance. But, and this is the statement, there is a high order energy e tilde that will allow you that will give you global existence for small data. However, you cannot classify the asymptotic attractor. So all you can say is that the time derivative, the tau derivative of the solution decays, but you cannot say that in the limit you belong to a nearby member of the linearly expanding family. This energy is constructed by using powers of this elliptic operator that I just wrote a second ago. And so to be honest, I can talk about it. I have a few more slides at the end, but let me not go into the technical aspects of this proof. Let me just say that this operator is designed to capture very precisely. So in other words, you cannot just blindly apply dr derivative or dx. If you do that, you will create singularities at zero and they are hard to handle. So you have to do it with a particular combination. And this operator s is a reflection of the structure of that weighted-generate operator L delta. And let's leave it at that. So I will not go into the technical details. Yeah, that's a good question. So it's just something that sits nearby, but it's not one of the... But you have to identify the proof, right? Yeah, yeah. Well, you can prove that the time derivative of that is decaying. So that gives you convergence to something. Yeah, yeah. And it's nearby. Yeah. So the reminiscent of a result I heard recently in a talk by Falker, Falker Schuhe, and sort of in these expanding space times. Who's here, I believe? Yeah, very good. Let me spend the remaining three minutes discussing perhaps questions that are... I find very interesting that are open. So here's this kind of a silly lemma. I call it lemma. It's too ambitious to call it lemma. Do there exist collapsing or expanding compactly supported stars? So this is the key for other values of gamma. They should exist in the supercritical range, which is between 6 fifths and 4 thirds. This is an open question. What you can show, however, is that self-similarity and a spherical symmetry cannot coexist for this type of collapse. Lama being 4 thirds is absolutely critical for the two to be able to combine the two. And this is a very silly thing to show. It's just some conservation law type calculation. So the question is, can you construct exisymmetric self-similar collapse in this mass supercritical range? I think this is a very interesting problem. And now we come to the question that I believe Sergio asked. Are the Lane-Emdon stars that are mentioned before in the subcritical range, where they are shown to be conditionally stable, are they in fact stable? So a possible mechanism, and you will see this discussed sometimes in the physics literature, is some type of shock formation. This is not known, and this is really important. Let me remind you of this picture. The pink elephant in the room is, of course, this region here, the collapsing region. So as I said, it is not difficult to show that it's nonlinearly unstable. The moment you have correct well-posedness theory and you have a growing mode, it's very easy. But so again, this should be a lemma. But it's really unsatisfactory, because it doesn't tell you anything about the real instability. In particular, it could it be that the nearby guys decay, collapse at perhaps a corrected rate? Do they collapse at all? So this is just in spherical symmetry. You don't have to go outside spherical symmetry to ask these questions. And finally, I would like to sort of point out this viewpoint that I believe studying this family of problems is really the right thing, because you really want to have a robust understanding of how equations of state affect the collapse. This is partly motivated by general relativity. If you think of the most famous example of a collapsing star in general relativity, it's the Oppenheimer-Sneider solution, which corresponds to the case where you don't have any pressure. So by a result of Christodulu from 83, I believe, you can show that these are generically unstable solutions. If you go outside the homogeneity class, you generically form naked singularities. So you want to, well, you can devise various conclusions out of that. One conclusion that I like to draw out of that is that studying non-trivial equations of state is an important problem, and this is just the Newtonian version of the problem. But realistic equations of state are, this is just a final comment, is in principle quite a hard problem, even for physicists. So this is my favorite relativistic astrophysics book by Zeldovich and Noviko, where I sort of learned this stuff about. And they discuss one of the sort of high density regimes for the equations of state. And here's what Zeldovich says. So it reminds him of Avocenko's parody. So it is precisely like this. So if you look at the astrophysical literature, there's a huge disparity between what people think is the right equation of state. Yeah. So this is a teaser, I think of it as a teaser trailer for this wonderful textbook because it's a very serious book, but as I browse through it every now and then I find a gem of this type. Thank you. So the question is, how about? So it seems, since the non-year expanded rate is unstable, does it mean that basically the type of solutions which would converge to the stable one, which is a new one? Okay, so if you generically perturb a self-similar expanding rate, if you generically perturb it with perturbations that have positive energy, then they will converge to the linear guy. So it's generically unstable and these guys are stable really only in the co-dimension one sense. So the linear is stable. The linear is stable. The linear takes over, not the self-similar. Yeah. Yeah, yeah, exactly. So do the physicists know this part? I mean, I'm not... It should be of great interest to them, I guess. Well, if you're... If you're very well, I should be linear. Yeah, exactly. So if you read Goldreich and Weber, which was one of the starting points for this project, for some reason they really say, we will not be interested in these linear rates. They drop a constant of integration in a second order ODE, which allows them to go down one level and get an explicit solution, which is this. And they literally say, ah, we kind of don't want to look at that. So but this is the finding. I mean, part of the mess was also to understand the exact structure of this parameter space. Yes. The non-linear unstable result covers the whole regime of gamma written there or...? Actually it actually goes up to here. I did not tell you what happens at 6-fifths. At 6-fifths you still get a steady state, but it is supported all the way up to... It's infinite. It has an infinite support. And it turns out to be... So in this case, the steady state is exactly the steady state of the critical Quintic wave equation. So it's one of these guys. I don't know. Once there's some constant here, and then I believe one over two, something like that. So this is the rate of decay of this guy. So it is the most complicated case to prove the most important? In this work, well, it's quite different because here you are working on a compactly supported domain. So you really have to rely on the well-posed in this theory for vacuum interface problems. Here it's different. Here the weights are at infinity. So it has a flavor of it, but it's quite different. I'm not sure how to... So this result is also due to Jang. So this equality case, all of this is due to Jang. Sorry, in this range. I have to be careful. This is the instability range. What this result shows in a way is that the fourth-third case, which is linearly stable, is in fact widely unstable. So your result is radial from the stability. This is entirely radial. So the paper that will show up is entirely radial, but it holds also for non-radial perturbations. So I alluded earlier to that. Here's the reason why the radial case is simpler, not just for obvious reasons, but also because it reduces the whole problem to just a quasi-linear wave equation. When you work with a non-radial case, so the fluid is not irrotational anymore in general, so you have to account for the vorticity, which does not have this wave-like structure, but it satisfies a certain transport equation that allows you to control it. What makes the vorticity better than in the standard equation where you know that you cannot control the vorticity? Is it because of the? Well, it's a local in-time result. I mean, yeah. No, but is it draw? Yeah, yeah. So all I'm saying, I'm not sure if I understand your question, but all I'm saying is that Ah, the long time. For the long time result, you don't have a long time. Oh, we haven't produced it yet. No, no. But for the irrotational case, it will be straightforward. And for the vorticity in the long time result, I don't want to make the claim yet. I don't think- Because I don't see why it should be any better than in the regular case. Maybe. No, no, no. You have this damping effect, so it will be better. It will be better. But will it affect the vorticity? No. You have, for instance, for Friedman, Lomet, Roberts and Walker space times, you have a stability without irrotational assumption. So when you have a cosmological constant, for instance, and you perturb away from that, you see this. You can control the full flow. Of course, this is not a free boundary problem, because this free boundary problem in GR is more- Okay, but the mechanism is massive. Yeah, yeah. You will kill it with that. All right. So thank you. So we start again at 11.30. So let's go to the next one. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you. Thank you.
everything I will talk about is a joint work with Juhi Jang who is based at the University of system, it comes from astrophysics, it's a model of a star. And we will be studying the nonlinear stability of that class of solutions. So let me compressible Euler equations surrounded by should exist, why small perturbations of these There is an invariant rescaling of this problem. pressure and the potential, you can work it out then to the proofs of these results. As I resist not putting this on the slides. This is a will give you always a linearly expanding non-linearly stable in this self-similar frame, open problem. I will state it at the end as well. So part two, how to study the stability of the self-adjoint operator with respect to this difficult in the first one. The difficulty the question is, can you construct exisymmetric you can show that these are generically unstable Goldreich and Weber, which was one of the