We report results from global ideal MHD simulations that study thin accretion disks (with thermal scale height H/R = 0.1 and 0.05) threaded by net vertical magnetic fields. Our computations span ...three orders of magnitude in radius, extend all the way to the pole, and are evolved for more than 1000 innermost orbits. We find that (1) inward accretion occurs mostly in the upper magnetically dominated regions of the disk at z ∼ R, similar to predictions from some previous analytical work and the "coronal accretion" flows found in GRMHD simulations. (2) A quasi-static global field geometry is established in which flux transport by inflows at the surface is balanced by turbulent diffusion. The resulting field is strongly pinched inwards at the surface. A steady-state advection-diffusion model, with a turbulent magnetic Prandtl number of order unity, reproduces this geometry well. (3) Weak unsteady disk winds are launched beyond the disk corona with the Alfvén radius RA/R0 ∼ 3. Although the surface inflow is filamentary and the wind is episodic, we show that the time-averaged properties are well-described by steady-wind theory. Even with strong fields, β0 = 103 at the midplane initially, only 5% of the angular momentum transport is driven by the wind, and the wind mass flux from the inner decade of the radius is only ∼0.4% of the mass accretion rate. (4) Within the disk, most of the accretion is driven by the Rφ stress from the MRI and global magnetic fields. Our simulations have many applications to astrophysical accretion systems.
The design and implementation of a new framework for adaptive mesh refinement calculations are described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible ...and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both nonrelativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.
•A fourth-order accurate method for the numerical solution of the equations of ideal magnetohydrodynamics (MHD) is proposed.•The piecewise parabolic method (PPM) in the upwind constrained transport ...framework satisfies the divergence-free constraint.•Modern PPM limiters provide MHD shock-capturing abilities while maintaining full fourth-order accuracy for smooth profiles.•Fourth-order accuracy yields significant improvement in solution errors relative to traditional second-order schemes.
We present a fourth-order accurate finite volume method for the solution of ideal magnetohydrodynamics (MHD). The numerical method combines high-order quadrature rules in the solution of semi-discrete formulations of hyperbolic conservation laws with the upwind constrained transport (UCT) framework to ensure that the divergence-free constraint of the magnetic field is satisfied. A novel implementation of UCT that uses the piecewise parabolic method (PPM) for the reconstruction of magnetic fields at cell corners in 2D is introduced. The resulting scheme can be expressed as the extension of the second-order accurate constrained transport (CT) Godunov-type scheme that is currently used in the Athena astrophysics code. After validating the base algorithm on a series of hydrodynamics test problems, we present the results of multidimensional MHD test problems which demonstrate formal fourth-order convergence for smooth problems, robustness for discontinuous problems, and improved accuracy relative to the second-order scheme.
ABSTRACT We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of ...advanced, less diffusive Riemann solvers, in particular HLLC and HLLD. We also employ a staggered-mesh constrained transport algorithm suited for curvilinear coordinate systems in order to maintain the divergence-free constraint of the magnetic field. Our code is designed to work with arbitrary stationary spacetimes in one, two, or three dimensions, and we demonstrate its reliability through a number of tests. We also report on its promising performance and scalability.
We use global three-dimensional radiation magnetohydrodynamical simulations to study accretion disks onto a black hole with accretion rates varying from to . We initialize the disks with a weakly ...magnetized torus centered at either 50 or 80 gravitational radii, leading to self-consistent turbulence generated by the magnetorotational instability (MRI). The inner regions of all disks have radiation pressure ∼104-106 times the gas pressure. Nonaxisymmetric density waves that steepen into spiral shocks form as gas flows toward the black hole. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when the net vertical magnetic flux is sufficiently large. Outflows are formed with a speed of ∼0.1-0.4c. When the accretion rate is smaller than , outflows are launched from ∼10 gravitational radii, and the radiative efficiency is ∼5%-7%. For an accretion rate reaching , most of the funnel region near the rotation axis becomes optically thick, and the outflow is launched from beyond 50 gravitational radii. The radiative efficiency is reduced to 1%. We always find that the kinetic energy luminosity associated with the outflow is at most ∼15%-30% of the radiative luminosity. The mass flux in the outflow is ∼15%-50% of the net mass accretion rates. We discuss the implications of our simulation results on the observational properties of these disks.
The global evolution of protoplanetary disks (PPDs) has recently been shown to be largely controlled by the amount of poloidal magnetic flux threading the disk. The amount of magnetic flux must also ...coevolve with the disk, as a result of magnetic flux transport, a process that is poorly understood. In weakly ionized gas as in PPDs, magnetic flux is largely frozen in the electron fluid, except when resistivity is large. When the disk is largely laminar, we show that the relative drift between the electrons and ions (the Hall drift), and the ions and neutral fluids (ambipolar drift) can play a dominant role on the transport of magnetic flux. Using two-dimensional simulations that incorporate the Hall effect and ambipolar diffusion (AD) with prescribed diffusivities, we show that when large-scale poloidal field is aligned with disk rotation, the Hall effect rapidly drags magnetic flux inward at the midplane region, while it slowly pushes flux outward above/below the midplane. This leads to a highly radially elongated field configuration as a global manifestation of the Hall-shear instability. This field configuration further promotes rapid outward flux transport by AD at the midplane, leading to instability saturation. In quasi-steady state, magnetic flux is transported outward at approximately the same rate at all heights, and the rate is comparable to the Hall-free case. For anti-aligned field polarity, the Hall effect consistently transports magnetic flux outward, leading to a largely vertical field configuration in the midplane region. The field lines in the upper layer first bend radially inward and then outward to launch a disk wind. Overall, the net rate of outward flux transport is about twice as fast as that of the aligned case. In addition, the rate of flux transport increases with increasing disk magnetization. The absolute rate of transport is sensitive to disk microphysics, which remains to be explored in future studies.
ABSTRACT Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral ...patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three-dimensional (3D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the nonlinear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3D hydrodynamical models than images based on two-dimensional models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.
Ketamine: A tale of two enantiomers Jelen, Luke A; Young, Allan H; Stone, James M
Journal of Psychopharmacology,
02/2021, Letnik:
35, Številka:
2
Book Review, Journal Article
Recenzirano
Odprti dostop
The discovery of the rapid antidepressant effects of the dissociative anaesthetic ketamine, an uncompetitive N-Methyl-D-Aspartate receptor antagonist, is arguably the most important breakthrough in ...depression research in the last 50 years. Ketamine remains an off-label treatment for treatment-resistant depression with factors that limit widespread use including its dissociative effects and abuse potential. Ketamine is a racemic mixture, composed of equal amounts of (S)-ketamine and (R)-ketamine. An (S)-ketamine nasal spray has been developed and approved for use in treatment-resistant depression in the United States and Europe; however, some concerns regarding efficacy and side effects remain. Although (R)-ketamine is a less potent N-Methyl-D-Aspartate receptor antagonist than (S)-ketamine, increasing preclinical evidence suggests (R)-ketamine may have more potent and longer lasting antidepressant effects than (S)-ketamine, alongside fewer side effects. Furthermore, a recent pilot trial of (R)-ketamine has demonstrated rapid-acting and sustained antidepressant effects in individuals with treatment-resistant depression. Research is ongoing to determine the specific cellular and molecular mechanisms underlying the antidepressant actions of ketamine and its component enantiomers in an effort to develop future rapid-acting antidepressants that lack undesirable effects. Here, we briefly review findings regarding the antidepressant effects of ketamine and its enantiomers before considering underlying mechanisms including N-Methyl-D-Aspartate receptor antagonism, γ-aminobutyric acid-ergic interneuron inhibition, α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic receptor activation, brain-derived neurotrophic factor and tropomyosin kinase B signalling, mammalian target of rapamycin complex 1 and extracellular signal-regulated kinase signalling, inhibition of glycogen synthase kinase-3 and inhibition of lateral habenula bursting, alongside potential roles of the monoaminergic and opioid receptor systems.
Gaseous disks have been proposed as a mechanism for facilitating mergers of binary black holes. We explore circumbinary disk systems to determine the evolution of the central binary. To do so, we ...perform 3D, hydrodynamic, locally isothermal simulations of circumbinary disks on a Cartesian grid. We focus on binaries of equal mass ratios on fixed circular orbits. To investigate the orbital evolution of the binary, we examine the various torques exerted on the system. For the case where the disk plane and binary orbital plane are aligned, we find that the total torque is positive so that the semimajor axis of the binary increases. For the misaligned case, we run simulations with the binary orbital plane and disk midplane misaligned by 45° and find the same results-the binary grows. The timescale for the circumbinary disk to realign to the plane of the binary is consistent with the global viscous timescale of the disk.
Abstract
We use 2D (axisymmetric) and 3D hydrodynamic simulations to study Bondi–Hoyle–Lyttleton accretion with and without transverse upstream gradients. We mainly focus on the regime of high ...(upstream) Mach number, weak upstream gradients, and small accretor size, which is relevant to neutron star accretion in wind-fed supergiant X-ray binaries (SgXBs). We present a systematic exploration of the flow in this regime. When there are no upstream gradients, the flow is always stable regardless of accretor size or Mach number. For finite upstream gradients, there are three main types of behaviour: stable flow (small upstream gradient), turbulent unstable flow without a disc (intermediate upstream gradient), and turbulent flow with a disc-like structure (relatively large upstream gradient). When the accretion flow is turbulent, the accretion rate decreases non-convergently as the accretor size decreases. The flow is more prone to instability and the disc is less likely to form than previously expected; the parameters of most observed SgXBs place them in the regime of a turbulent, disc-less accretion flow. Among the SgXBs with relatively well-determined parameters, we find OAO 1657−415 to be the only one that is likely to host a persistent disc (or disc-like structure); this finding is consistent with observations.