Researcher profile

Matthew Russo

Matthew Russo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
7works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2015arXiv

Radially Magnetized Protoplanetary Disk: Vertical Profile

This paper studies the response of a thin accretion disk to an external radial magnetic field. Our focus is on protoplanetary disks (PPDs), which are exposed during their later evolution to an intense, magnetized wind from the central star. A radial magnetic field is mixed into a thin surface layer, is wound up by the disk shear, and is pushed downward by a combination of turbulent mixing and ambipolar and Ohmic drift. The toroidal field reaches much greater strengths than the seed vertical field that is usually invoked in PPD models, even becoming superthermal. Linear stability analysis indicates that the disk experiences the magnetorotational instability (MRI) at a higher magnetization than a vertically magnetized disk when both the effects of ambipolar and Hall drift are taken into account. Steady vertical profiles of density and magnetic field are obtained at several radii between 0.06 and 1 AU in response to a wind magnetic field $B_r \sim (10^{-4}$-$10^{-2})(r/{\rm AU})^{-2}$ G. Careful attention is given to the radial and vertical ionization structure resulting from irradiation by stellar X-rays. The disk is more strongly magnetized closer to the star, where it can support a higher rate of mass transfer. As a result, the inner $\sim 1$ AU of a PPD is found to evolve toward lower surface density. Mass transfer rates around $10^{-8}\,M_\odot$ yr$^{-1}$ are obtained under conservative assumptions about the MRI-generated stress. The evolution of the disk, and the implications for planet migration, are investigated in the accompanying paper.

preprint2014arXiv

Integrable Generalized KdV and MKdV Equations with Spatiotemporally Varying Coefficients

A technique based on extended Lax Pairs is first considered to derive variable-coefficient generalizations of various Lax-integrable NLPDE hierarchies recently introduced in the literature. As illustrative examples, we consider generalizations of KdV equations and three variants of generalized MKdV equations. It is demonstrated that the techniques yield Lax- or S-integrable NLPDEs with both time- AND space-dependent coefficients which are thus more general than almost all cases considered earlier via other methods such as the Painleve Test, Bell Polynomials, and various similarity methods. However, this technique, although operationally effective, has the significant disadvantage that, for any integrable system with spatiotemporally varying coefficients, one must guess a generalization of the structure of the known Lax Pair for the corresponding system with constant coefficients. Motivated by the somewhat arbitrary nature of the above procedure, we embark in this paper on an attempt to systematize the derivation of Lax-integrable sytems with variable coefficients. Hence we attempt to apply the Estabrook-Wahlquist (EW) prolongation technique, a relatively self-consistent procedure requiring little prior information. However, this immediately requires that the technique be significantly generalized or broadened in several different ways, including solving matrix partial differential equations instead of algebraic ones. The new and extended EW technique whch results is illustrated by algorithmically deriving generalized Lax-integrable versions of the generalized fifth-order KdV, and MKdV equations.

preprint2014arXiv

Integrable Spatiotemporally Varying NLS, PT-Symmetric NLS, and DNLS Equations: Generalized Lax Pairs and Lie Algebras

This paper develops two approaches to Lax-integrbale systems with spatiotemporally varying coefficients. A technique based on extended Lax Pairs is first considered to derive variable-coefficient generalizations of various Lax-integrable NLPDE hierarchies recently introduced in the literature. As illustrative examples, we consider generalizations of the NLS and DNLS equations, as well as a PT-symmetric version of the NLS equation. It is demonstrated that the techniques yield Lax- or S-integrable NLPDEs with both time- AND space-dependent coefficients which are thus more general than almost all cases considered earlier via other methods such as the Painleve Test, Bell Polynomials, and various similarity methods. However, this technique, although operationally effective, has the significant disadvantage that, for any integrable system with spatiotemporally varying coefficients, one must 'guess' a generaliza- tion of the structure of the known Lax Pair for the corresponding system with constant coefficients. Motivated by the somewhat arbitrary nature of the above procedure, we therefore next attempt to systematize the derivation of Lax-integrable sytems with variable coefficients. We attempt to apply the Estabrook- Wahlquist (EW) prolongation technique, a relatively self-consistent procedure requiring little prior infomation. However, this immediately requires that the technique be significantly generalized or broadened in several different ways. The new and extended EW technique which results is illustrated by algorithmically deriving generalized Lax-integrable versions of NLS, PT-symmetric NLS, and DNLS equations.

preprint2014arXiv

The Extended Estabrook-Wahlquist Method

Variable Coefficient Korteweg de Vries (vcKdV), Modified Korteweg de Vries (vcMKdV), and nonlinear Schrodinger (NLS) equations have a long history dating from their derivation in various applications. A technique based on extended Lax Pairs has been devised recently to derive variable-coefficient generalizations of various Lax-integrable NLPDE hierarchies. The resulting Lax- or S-integrable NLPDEs have both time- AND space-dependent coefficients, and are thus more general than almost all cases considered earlier via other methods such as the Painlevé Test, Bell Polynomials, and various similarity methods. However, this technique, although operationally effective, has the significant disadvantage that, for any integrable system with spatiotemporally varying coefficients, one must 'guess' a generalization of the structure of the known Lax Pair for the corresponding system with constant coefficients. Motivated by the somewhat arbitrary nature of the above procedure, we embark in this paper on an attempt to systematize the derivation of Lax-integrable systems with variable coefficients. An ideal approach would be a method which does not require knowledge of the Lax pair to an associated constant coefficient system, and also involves little to no guesswork. Hence we attempt to apply the Estabrook-Wahlquist (EW) prolongation technique, a relatively self-consistent procedure requiring little prior information. The new and extended EW technique which results is illustrated by algorithmically deriving generalized Lax-integrable versions of the NLS, generalized fifth-order KdV, MKdV, and derivative nonlinear Schrodinger (DNLS) equations.

preprint2013arXiv

Hot Electromagnetic Outflows I: Acceleration and Spectra

The theory of cold, relativistic, magnetohydrodynamic outflows is generalized by the inclusion of an intense radiation source. In some contexts, such the breakout of a gamma-ray burst jet from a star, the outflow is heated to a high temperature at a large optical depth. Eventually it becomes transparent and is pushed to a higher Lorentz factor by a combination of the Lorentz force and radiation pressure. We obtain its profile, both inside and outside the fast magnetosonic critical point, when the poloidal magnetic field is radial and monopolar. Most of the energy flux is carried by the radiation field and the toroidal magnetic field that is wound up close to the rapidly rotating engine. Although the entrained matter carries little energy, it couples the radiation field to the magnetic field. Then the fast critical point is pushed inward from infinity and, above a critical radiation intensity, the outflow is accelerated mainly by radiation pressure. We identify a distinct observational signature of this hybrid outflow: a hardening of the radiation spectrum above the peak of the seed photon distribution, driven by bulk Compton scattering. The non-thermal spectrum -- obtained by a Monte Carlo method -- is most extended when the Lorentz force dominates the acceleration, and the seed photon beam is wider than the Lorentz cone of the MHD fluid. This effect is a generic feature of hot, magnetized outflows interacting with slower relativistic material. It may explain why some GRB spectra appear to peak at photon energies above the original Amati et al. scaling. A companion paper addresses the case of jet breakout, where diverging magnetic flux surfaces yield strong MHD acceleration over a wider range of Lorentz factor.

preprint2013arXiv

Hot Electromagnetic Outflows II: Jet Breakout

We consider the interaction between radiation, matter and a magnetic field in a compact, relativistic jet. The entrained matter accelerates outward as the jet breaks out of a star or other confining medium. In some circumstances, such as gamma-ray bursts (GRBs), the magnetization of the jet is greatly reduced by an advected radiation field while the jet is optically thick to scattering. Where magnetic flux surfaces diverge rapidly, a strong outward Lorentz force develops and radiation and matter begin to decouple. The increase in magnetization is coupled to a rapid growth in Lorentz factor. We take two approaches to this problem. The first examines the flow outside the fast magnetosonic critical surface, and calculates the flow speed and the angular distribution of the radiation field over a range of scattering depths. The second considers the flow structure on both sides of the critical surface in the optically thin regime, using a relaxation method. In both approaches, we find how the terminal Lorentz factor, and radial profile of the outflow, depend on the radiation intensity and optical depth at breakout. The effect of bulk Compton scattering on the radiation spectrum is calculated by a Monte Carlo method, while neglecting the effects of internal dissipation. The peak of the scattered spectrum sits near the seed peak if radiation pressure dominates the acceleration, but is pushed to a higher frequency if the Lorentz force dominates, and especially if the seed photon cone is broadened by interaction with a slower component of the outflow.

preprint2011arXiv

Light Loop Echoes and Blinking Black Holes

Radiation emitted near a black hole reaches the observer by multiple paths; and when this radiation varies in time, the time-delays between the various paths generate a &#34;blinking&#34; effect in the observed light curve L(t) or its auto-correlation function xi(T)= <L(t)L(t-T)>. For the particularly important &#34;face-on&#34; configuration (in which the hole is viewed roughly along its spin axis, while the emission comes roughly from its equatorial plane -- e.g. from the inner edge of its accretion disk, or from the violent flash of a nearby/infalling star) we calculate the blinking in detail by computing the time delay Delta t_{j}(r,a) and magnification mu_{j}(r,a) of the jth path (j=1,2,3,...), relative to the primary path (j=0), as a function of the emission radius r and black hole spin 0<a/M<1. The particular geometry and symmetry of the nearly-face-on configuration enhances and &#34;protects&#34; the blinking signal, making it more detectable and more independent of certain astrophysical and observational details. The effect can be surprisingly strong: e.g. for radiation from the innermost stable circular orbit (&#34;ISCO&#34;) of a black hole of critical spin (a_{crit}/M = 0.853), the j=1,2,3 fluxes are, respectively, 27%, 2% and 0.1% of the j=0 flux.