-
Stability of hypermassive neutron stars with realistic rotation and entropy profiles
Authors:
Nishad Muhammed,
Matthew D. Duez,
Pavan Chawhan,
Noora Ghadiri,
Luisa T. Buchman,
Francois Foucart,
Patrick Chi-Kit Cheong,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Mark A. Scheel
Abstract:
Binary neutron star mergers produce massive, hot, rapidly differentially rotating neutron star remnants; electromagnetic and gravitational wave signals associated with the subsequent evolution depend on the stability of these remnants. Stability of relativistic stars has previously been studied for uniform rotation and for a class of differential rotation with monotonic angular velocity profiles.…
▽ More
Binary neutron star mergers produce massive, hot, rapidly differentially rotating neutron star remnants; electromagnetic and gravitational wave signals associated with the subsequent evolution depend on the stability of these remnants. Stability of relativistic stars has previously been studied for uniform rotation and for a class of differential rotation with monotonic angular velocity profiles. Stability of those equilibria to axisymmetric perturbations was found to respect a turning point criterion: along a constant angular momentum sequence, the onset of unstable stars is found at maximum density less than but close to the density of maximum mass. In this paper, we test this turning point criterion for non-monotonic angular velocity profiles and non-isentropic entropy profiles, both chosen to more realistically model post-merger equilibria. Stability is assessed by evolving perturbed equilibria in 2D using the Spectral Einstein Code. We present tests of the code's new capability for axisymmetric metric evolution. We confirm the turning point theorem and determine the region of our rotation law parameter space that provides highest maximum mass for a given angular momentum.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
High angular momentum hot differentially rotating equilibrium star evolutions in conformally flat spacetime
Authors:
Patrick Chi-Kit Cheong,
Nishad Muhammed,
Pavan Chawhan,
Matthew D. Duez,
Francois Foucart,
Lawrence E. Kidder,
Harald P. Pfeiffer,
Mark A. Scheel
Abstract:
The conformal flatness approximation to the Einstein equations has been successfully used in many astrophysical applications such as initial data constructions and dynamical simulations. Although it has been shown that full general relativistic strongly differentially rotating equilibrium models deviate by at most a few percent from their conformally flat counterparts, whether those conformally fl…
▽ More
The conformal flatness approximation to the Einstein equations has been successfully used in many astrophysical applications such as initial data constructions and dynamical simulations. Although it has been shown that full general relativistic strongly differentially rotating equilibrium models deviate by at most a few percent from their conformally flat counterparts, whether those conformally flat solutions remain stable has not been fully addressed. To further understand the limitations of the conformal flatness approximation, in this work, we construct spatially-conformally-flat hot hypermassive neutron stars with post-merger-like rotation laws, and perform conformally flat evolutions and analysis over dynamical timescales. We find that enforcing conformally-flat spacetime could change the equilibrium of quasi-toroidal models with high angular momentum for $J \gtrsim 9 ~G M_{\odot}^2 / c$ compared to fully general relativistic cases. In contrast, all the quasi-spherical models considered in this work remain stable even with high angular momentum $J=9~G M_{\odot}^2 / c$. Our investigation suggests that the quasi-spherical models are suitable initial data for long-lived hypermassive neutron star modeling in conformally flat spacetime.
△ Less
Submitted 2 May, 2024; v1 submitted 28 February, 2024;
originally announced February 2024.
-
Universal relations for fundamental modes of rotating neutron stars with differential rotations
Authors:
Anson Ka Long Yip,
Patrick Chi-Kit Cheong,
Tjonnie Guang Feng Li
Abstract:
Violent astrophysical events, including core-collapse supernovae and binary neutron star mergers, can result in rotating neutron stars with diverse degrees of differential rotation. Oscillation modes of these neutron stars could be excited and emit strong gravitational waves. Detecting these modes may provide information about neutron stars, including their structures and dynamics. Hence, dynamica…
▽ More
Violent astrophysical events, including core-collapse supernovae and binary neutron star mergers, can result in rotating neutron stars with diverse degrees of differential rotation. Oscillation modes of these neutron stars could be excited and emit strong gravitational waves. Detecting these modes may provide information about neutron stars, including their structures and dynamics. Hence, dynamical simulations were employed to construct relations for quantifying the oscillation mode frequency in previous studies. Specifically, linear relations for the frequencies of fundamental $l=0$ quasi-radial mode $f_{F}$ and fundamental $l=2$ quadrupolar mode $f_{^2f}$ were constructed by simulations with the Cowling approximation. Nevertheless, these relations can overestimate $f_{F}$ and underestimate $f_{^2f}$ up to $\sim 30\%$. Furthermore, it has yet to be fully studied how the degree of differential rotation affects $f_{F}$ and $f_{^2f}$. Here, for the first time, we consider both various degrees of differential rotation $\tilde{A}$ and dynamical spacetime to construct linear relations for quantifying $f_{F}$ and $f_{^2f}$. Through 2D axisymmetric simulations, we first show that both $f_{F}$ and $f_{^2f}$ scale almost linearly with the stellar compactness $M/R$ for different values of $\tilde{A}$. We also observe the quasi-linear relations for both $f_{F}$ and $f_{^2f}$ with the kinetic-to-binding energy ratio $T/|W|$ for different $\tilde{A}$ values. Finally, we constructed linear fits that can quantify $f_{F}$ and $f_{^2f}$ by $T/|W|$. Consequently, this work presented universal relations for the fundamental modes of rotating neutron stars with differential rotations in dynamical spacetime.
△ Less
Submitted 25 January, 2024;
originally announced January 2024.
-
General-relativistic radiation transport scheme in Gmunu II: Implementation of novel microphysical library for neutrino radiation -- Weakhub
Authors:
Harry Ho-Yin Ng,
Patrick Chi-Kit Cheong,
Alan Tsz-Lok Lam,
Tjonnie Guang Feng Li
Abstract:
We introduce \texttt{Weakhub}, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino-matter, neutrino-photon, and neutrino-neutrino interactions, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of $β$-processes, incorporating v…
▽ More
We introduce \texttt{Weakhub}, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino-matter, neutrino-photon, and neutrino-neutrino interactions, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of $β$-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino-antineutrino annihilation, and nuclear de-excitation are included. We also present the detailed derivations of weak interactions and the coupling of them to the two-moment based general-relativistic multi-group radiation transport in the \texttt{G}eneral-relativistic \texttt{mu}ltigrid \texttt{nu}merical (\texttt{Gmunu}) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernova and binary neutron star postmerger remnant, and predict the effects of improved opacities in comparison to conventional ones for a binary neutron star postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library \texttt{NuLib} in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than $\sim 10\%$ in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modelling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.
△ Less
Submitted 1 December, 2023; v1 submitted 7 September, 2023;
originally announced September 2023.
-
Gravitational wave signatures from the phase-transition-induced collapse of a magnetized neutron star
Authors:
Anson Ka Long Yip,
Patrick Chi-Kit Cheong,
Tjonnie Guang Feng Li
Abstract:
Strong magnetic fields make neutron stars potential sources of detectable electromagnetic and gravitational-wave signals. Hence, inferring these magnetic fields is critical to understand the emissions of neutron stars. However, due to the lack of direct observational evidence, the interior magnetic field configuration remains ambiguous. Here, for the first time, we show that the internal magnetic…
▽ More
Strong magnetic fields make neutron stars potential sources of detectable electromagnetic and gravitational-wave signals. Hence, inferring these magnetic fields is critical to understand the emissions of neutron stars. However, due to the lack of direct observational evidence, the interior magnetic field configuration remains ambiguous. Here, for the first time, we show that the internal magnetic field strength along with the composition of a neutron star can be directly constrained by detecting the gravitational waves from the phase-transition-induced collapse of a magnetized neutron star. By dynamically simulating this collapsing event, we first find that the dominant peaks in the gravitational waveform are the fundamental $l=0$ quasi-radial $F$ mode and the fundamental $l=2$ quadrupolar $^2f$ mode. We next show that the maximum gravitational wave amplitude $|h|_\mathrm{max}$ increases with the maximum magnetic field strength of the interior toroidal field $\mathcal{B}_\mathrm{max}$ until the maximum rest-mass density at bounce $ρ_\mathrm{max,b}$ decreases due to the increasing $\mathcal{B}_\mathrm{max}$. We then demonstrated that the magnetic suppression of fundamental modes found in our previous work remains valid for the hybrid stars formed after the phase-transition-induced collapses. We finally show that measuring the frequency ratio between the two fundamental modes $f_{^2f}/f_{F}$ allows one to infer $\mathcal{B}_\mathrm{max}$ and the baryonic mass fraction of matter in the mixed phase $M_\mathrm{mp} / M_{0}$ of the resulting hybrid star. Consequently, taking $\mathcal{B}_\mathrm{max}$ and $M_\mathrm{mp} / M_{0}$ as examples, this work has demonstrated that much information inside neutron stars could be extracted similarly through measuring the oscillation modes of the stars.
△ Less
Submitted 24 May, 2023;
originally announced May 2023.
-
General-relativistic simulations of the formation of a magnetized hybrid star
Authors:
Anson Ka Long Yip,
Patrick Chi-Kit Cheong,
Tjonnie Guang Feng Li
Abstract:
Strongly magnetized neutron stars are popular candidates for producing detectable electromagnetic and gravitational-wave signals. A rapid density increase in a neutron star core could also trigger the phase transition from hadrons to deconfined quarks and form a hybrid star. This formation process could release a considerable amount of energy in the form of gravitational waves and neutrinos. Hence…
▽ More
Strongly magnetized neutron stars are popular candidates for producing detectable electromagnetic and gravitational-wave signals. A rapid density increase in a neutron star core could also trigger the phase transition from hadrons to deconfined quarks and form a hybrid star. This formation process could release a considerable amount of energy in the form of gravitational waves and neutrinos. Hence, the formation of a magnetized hybrid star is an interesting scenario for detecting all these signals. These detections may provide essential probes for the magnetic field and composition of such stars. Thus far, a dynamical study of the formation of a magnetized hybrid star has yet to be realized. Here, we investigate the formation dynamics and the properties of a magnetized hybrid star through dynamical simulations. We find that the maximum values of rest-mass density and magnetic field strength increase slightly and these two quantities are coupled in phase during the formation. We then demonstrate that all microscopic and macroscopic quantities of the resulting hybrid star vary dramatically when the maximum magnetic field strength goes beyond a threshold of $\sim 5 \times 10^{17}$ G but they are insensitive to the magnetic field below this threshold. Specifically, the magnetic deformation makes the rest-mass density drop significantly, suppressing the matter fraction in the mixed phase. Therefore, this work provides a solid support for the magnetic effects on a hybrid star, so it is possible to link observational signals from the star to its magnetic field configuration.
△ Less
Submitted 29 March, 2023;
originally announced March 2023.
-
Oscillations of Highly Magnetized Non-rotating Neutron Stars
Authors:
Man Yin Leung,
Anson Ka Long Yip,
Patrick Chi-Kit Cheong,
Tjonnie Guang Feng Li
Abstract:
Highly magnetized neutron stars are promising candidates to explain some of the most peculiar astronomical phenomena, for instance, fast radio bursts, gamma-ray bursts, and superluminous supernovae. Pulsations of these highly magnetized neutron stars are also speculated to produce detectable gravitational waves. In addition, pulsations are important probes of the structure and equation of state of…
▽ More
Highly magnetized neutron stars are promising candidates to explain some of the most peculiar astronomical phenomena, for instance, fast radio bursts, gamma-ray bursts, and superluminous supernovae. Pulsations of these highly magnetized neutron stars are also speculated to produce detectable gravitational waves. In addition, pulsations are important probes of the structure and equation of state of the neutron stars. The major challenge in studying the pulsations of highly magnetized neutron stars is the demanding numerical cost of consistently solving the nonlinear Einstein and Maxwell equations under minimum assumptions. With the recent breakthroughs in numerical solvers, we investigate pulsation modes of non-rotating neutron stars which harbour strong purely toroidal magnetic fields of $10^{15-17}$ G through two-dimensional axisymmetric general-relativistic magnetohydrodynamics simulations. We show that stellar oscillations are insensitive to magnetization effects until the magnetic to binding energy ratio goes beyond 10%, where the pulsation mode frequencies are strongly suppressed. We further show that this is the direct consequence of the decrease in stellar compactness when the extreme magnetic fields introduce strong deformations of the neutron stars.
△ Less
Submitted 9 March, 2023;
originally announced March 2023.
-
General-relativistic radiation transport scheme in $\texttt{Gmunu}$ I: Implementation of two-moment based multi-frequency radiative transfer and code tests
Authors:
Patrick Chi-Kit Cheong,
Harry Ho-Yin Ng,
Alan Tsz-Lok Lam,
Tjonnie Guang Feng Li
Abstract:
We present the implementation of two-moment based general-relativistic multi-group radiation transport module in the $\texttt{G}$eneral-relativistic $\texttt{mu}$ltigrid $\texttt{nu}$merical ($\texttt{Gmunu}$) code. On top of solving the general-relativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth-…
▽ More
We present the implementation of two-moment based general-relativistic multi-group radiation transport module in the $\texttt{G}$eneral-relativistic $\texttt{mu}$ltigrid $\texttt{nu}$merical ($\texttt{Gmunu}$) code. On top of solving the general-relativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth- and first-order moments of the radiations. Analytic closure relation is used to obtain the higher order moments and close the system. The finite-volume discretisation has been adopted for the radiation moments. The advection in spatial and frequency spaces are handled explicitly. In addition, the radiation-matter interaction terms, which are very stiff in the optically thick region, are solved implicitly. Implicit-explicit Runge-Kutta schemes are adopted for time integration. We test the implementation with a number of numerical benchmarks from frequency-integrated to frequency dependent cases. Furthermore, we also illustrate the astrophysical applications in hot neutron star and core-collapse supernovae modellings, and compare with other neutrino transport codes.
△ Less
Submitted 8 August, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS
Authors:
Travis J. Thieme,
Shih-Ping Lai,
Sheng-Jun Lin,
Pou-Ieng Cheong,
Chin-Fei Lee,
Hsi-Wei Yen,
Zhi-Yun Li,
Ka Ho Lam,
Bo Zhao
Abstract:
Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flo…
▽ More
Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized due to their low intensity and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow-like structures in C$^{18}$O that follow the edges of the outflows. We make various types of position-velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five sub-structures in position-position-velocity space. Four out of the five sub-structures fit well ($>$95%) with our simple infall model, with specific angular momenta between $2.7-6.9\times10^{-4}\,$km$\,$s$^{-1}\,$pc and mass-infall rates of $0.5-1.1\times10^{-6}\,M_{\odot}\,$yr$^{-1}$. Better characterization of the physical structure in the supposed "outflow-cavities" is important to disentangle the true outflow cavities and accretion flows.
△ Less
Submitted 7 November, 2021;
originally announced November 2021.
-
An extension of Gmunu: General-relativistic resistive magnetohydrodynamics based on staggered-meshed constrained transport with elliptic cleaning
Authors:
Patrick Chi-Kit Cheong,
David Yat Tung Pong,
Anson Ka Long Yip,
Tjonnie Guang Feng Li
Abstract:
We present the implementation of general-relativistic resistive magnetohydrodynamics solvers and three divergence-free handling approaches adopted in the General-relativistic multigrid numerical (Gmunu) code.
In particular, implicit-explicit Runge-Kutta schemes are used to deal with the stiff terms in the evolution equations for small resistivity.
Three divergence-free handling methods are (i)…
▽ More
We present the implementation of general-relativistic resistive magnetohydrodynamics solvers and three divergence-free handling approaches adopted in the General-relativistic multigrid numerical (Gmunu) code.
In particular, implicit-explicit Runge-Kutta schemes are used to deal with the stiff terms in the evolution equations for small resistivity.
Three divergence-free handling methods are (i) hyperbolic divergence cleaning through a generalised Lagrange multiplier (GLM); (ii) staggered-meshed constrained transport (CT) schemes and (iii) elliptic cleaning though multigrid (MG) solver which is applicable in both cell-centred and face-centred (stagger grid) magnetic field.
The implementation has been test with a number of numerical benchmarks from special-relativistic to general-relativistic cases.
We demonstrate that our code can robustly recover a very wide range of resistivity.
We also illustrate the applications in modelling magnetised neutron stars, and compare how different divergence-free handling affects the evolution of the stars.
Furthermore, we show that the preservation of the divergence-free condition of magnetic field when staggered-meshed constrained transport schemes can be significantly improved by applying elliptic cleaning.
△ Less
Submitted 20 January, 2022; v1 submitted 7 October, 2021;
originally announced October 2021.
-
Gravitational-wave asteroseismology with f-modes from neutron star binaries at the merger phase
Authors:
Harry Ho-Yin Ng,
Patrick Chi-Kit Cheong,
Lap-Ming Lin,
Tjonnie Guang Feng Li
Abstract:
Gravitational-wave signals from coalescing binary neutron stars can yield important information about the properties of nuclear-matter equation of state from the early part of the signal through tidal effects to the properties and oscillation frequencies of the merger product. We investigate a direct link between the properties of isolated neutron stars and their merger, by comparing the frequency…
▽ More
Gravitational-wave signals from coalescing binary neutron stars can yield important information about the properties of nuclear-matter equation of state from the early part of the signal through tidal effects to the properties and oscillation frequencies of the merger product. We investigate a direct link between the properties of isolated neutron stars and their merger, by comparing the frequency of the fundamental oscillation mode ($f$-mode) of neutron stars with the gravitational-wave frequency associated with the merger of two neutron stars. We calculate the quadrupolar ($l=2$) $f$-mode oscillation ($f_{2f}$) of non-rotating and rotating neutron stars using a nonlinear hydrodynamics code in the conformally-flat approximation and obtain the gravitational-wave frequency associated with the peak amplitude ($f_{\rm max}$) of binary-neutron stars from a set of publicly available simulations. We find that $f_{\rm max}$ and $f_{2f}$ differ by about 1\%, on average, across forty-five equal-mass systems with different total mass and equations of state. Assuming that the gravitational-wave frequency is still approximately equal to twice the orbital frequency $Ω$ near the merger, the result indicates that the condition for tidal resonance $|m|Ω= f_{2f}$ is satisfied to high accuracy near the merger, where $m=2$ is the azimuthal quantum number. Moreover, the well established universal relation between $f_{\rm max}$ and the tidal deformability of equal-mass binary systems can now be explained by a similar relation between $f_{2f}$ and the tidal deformability of isolated neutron stars, which has been demonstrated to be associated with the nearly incompressible properties of neutron stars. Our findings suggest that it is possible to relate the gravitational-wave signal at the merger of a binary neutron star system directly to the fundamental oscillation modes and the mass ratio.
△ Less
Submitted 15 December, 2020;
originally announced December 2020.
-
Gmunu: Paralleled, grid-adaptive, general-relativistic magnetohydrodynamics in curvilinear geometries in dynamical spacetimes
Authors:
Patrick Chi-Kit Cheong,
Alan Tsz-Lok Lam,
Harry Ho-Yin Ng,
Tjonnie Guang Feng Li
Abstract:
We present an update of the General-relativistic multigrid numerical (Gmunu) code, a parallelized, multi-dimensional curvilinear, general relativistic magnetohydrodynamics code with an efficient non-linear cell-centred multigrid (CCMG) elliptic solver, which is fully coupled with an efficient block-based adaptive mesh refinement modules. Currently, Gmunu is able to solve the elliptic metric equati…
▽ More
We present an update of the General-relativistic multigrid numerical (Gmunu) code, a parallelized, multi-dimensional curvilinear, general relativistic magnetohydrodynamics code with an efficient non-linear cell-centred multigrid (CCMG) elliptic solver, which is fully coupled with an efficient block-based adaptive mesh refinement modules. Currently, Gmunu is able to solve the elliptic metric equations in the conformally flat condition (CFC) approximation with the multigrid approach and the equations of ideal general-relativistic magnetohydrodynamics by means of high-resolution shock-capturing finite volume method with reference-metric formularise multi-dimensionally in cartesian, cylindrical or spherical geometries. To guarantee the absence of magnetic monopoles during the evolution, we have developed an elliptical divergence cleaning method by using multigrid solver. In this paper, we present the methodology, full evolution equations and implementation details of our code Gmunu and its properties and performance in some benchmarking and challenging relativistic magnetohydrodynamics problems.
△ Less
Submitted 10 September, 2021; v1 submitted 14 December, 2020;
originally announced December 2020.
-
Gmunu: Toward multigrid based Einstein field equations solver for general-relativistic hydrodynamics simulations
Authors:
Patrick Chi-Kit Cheong,
Lap-Ming Lin,
Tjonnie Guang-Feng Li
Abstract:
We present a new open-source axisymmetric general relativistic hydrodynamics code Gmunu (General-relativistic multigrid numerical solver) which uses a multigrid method to solve the elliptic metric equations in the conformally flat condition (CFC) approximation on a spherical grid. Most of the existing relativistic hydrodynamics codes are based on formulations which rely on a free-evolution approac…
▽ More
We present a new open-source axisymmetric general relativistic hydrodynamics code Gmunu (General-relativistic multigrid numerical solver) which uses a multigrid method to solve the elliptic metric equations in the conformally flat condition (CFC) approximation on a spherical grid. Most of the existing relativistic hydrodynamics codes are based on formulations which rely on a free-evolution approach of numerical relativity, where the metric variables are determined by hyperbolic equations without enforcing the constraint equations in the evolution. On the other hand, although a fully constrained-evolution formulation is theoretical more appealing and should lead to more stable and accurate simulations, such an approach is not widely used because solving the elliptic-type constraint equations during the evolution is in general more computationally expensive than hyperbolic free-evolution schemes. Multigrid methods solve differential equations with a hierarchy of discretizations and its computational cost is generally lower than other methods such as direct methods, relaxation methods, successive over-relaxation. With multigrid acceleration, one can solve the metric equations on a comparable time scale as solving the hydrodynamics equations. This would potentially make a fully constrained-evolution formulation more affordable in numerical relativity simulations. As a first step to assess the performance and robustness of multigrid methods in relativistic simulations, we develop a hydrodynamics code that makes use of standard finite-volume methods coupled with a multigrid metric solver to solve the Einstein equations in the CFC approximation. In this paper, we present the methodology and implementation of our code Gmunu and its properties and performance in some benchmarking relativistic hydrodynamics problems.
△ Less
Submitted 23 September, 2020; v1 submitted 16 January, 2020;
originally announced January 2020.
-
Ultrabroadband Density of States of Amorphous In-Ga-Zn-O
Authors:
Kyle T. Vogt,
Christopher E. Malmberg,
Jacob C. Buchanan,
George W. Mattson,
G. Mirek Brandt,
Dylan B. Fast,
Paul H. -Y. Cheong,
John F. Wager,
Matt W. Graham
Abstract:
The sub-gap density of states of amorphous indium gallium zinc oxide ($a$-IGZO) is obtained using the ultrabroadband photoconduction (UBPC) response of thin-film transistors (TFTs). Density functional theory simulations classify the origin of the measured sub-gap density of states peaks as a series of donor-like oxygen vacancy states and acceptor-like Zn vacancy states. Donor peaks are found both…
▽ More
The sub-gap density of states of amorphous indium gallium zinc oxide ($a$-IGZO) is obtained using the ultrabroadband photoconduction (UBPC) response of thin-film transistors (TFTs). Density functional theory simulations classify the origin of the measured sub-gap density of states peaks as a series of donor-like oxygen vacancy states and acceptor-like Zn vacancy states. Donor peaks are found both near the conduction band and deep in the sub-gap, with peak densities of $10^{17}-10^{18}$ cm$^{-3}$eV$^{-1}$. Two deep acceptor-like metal vacancy peaks with peak densities in the range of $10^{18}$ cm$^{-3}$eV$^{-1}$ and lie adjacent to the valance band Urbach tail region at 2.0 to 2.5 eV below the conduction band edge. By applying detailed charge balance, we show increasing the density of metal vacancy deep-acceptors strongly shifts the $a$-IGZO TFT threshold voltage to more positive values. Photoionization (h$ν$ > 2.0 eV) of metal vacancy acceptors is one cause of transfer curve hysteresis in $a$-IGZO TFTs owing to longer recombination lifetimes as they get captured into acceptor-like vacancies.
△ Less
Submitted 3 July, 2020; v1 submitted 14 October, 2019;
originally announced October 2019.
-
Determining the physical conditions of extremely young Class 0 circumbinary disk around VLA1623A
Authors:
Cheng-Han Hsieh,
Shih-Ping Lai,
Pou-Ieng Cheong,
Chia-Lin Ko,
Zhi-Yun Li,
Nadia M. Murillo
Abstract:
We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbina…
▽ More
We present detailed analysis of high-resolution C18O (2-1), SO (88-77), CO (3-2) and DCO+ (3-2) data obtained by the Atacama Large Millimeter/sub-millimeter Array (ALMA) towards a Class 0 Keplerian circumbinary disk around VLA1623A, which represents one of the most complete analysis towards a Class 0 source. From the dendrogram analysis, we identified several accretion flows feeding the circumbinary disk in a highly anisotropic manner. Stream-like SO emission around the circumbinary disk reveals the complicated shocks caused by the interactions between the disk, accretion flows and outflows. A wall-like structure is discovered south of VLA1623B. The discovery of two outflow cavity walls at the same position traveling at different velocities suggests the two outflows from both VLA1623A and VLA1623B overlays on top of each other in the plane of sky. Our detailed flat and flared disk modeling shows that Cycle 2 C18O J = 2-1 data is inconsistent with the combined binary mass of 0.2 Msun as suggested by early Cycle 0 studies. The combined binary mass for VLA1623A should be modified to 0.3 ~ 0.5 Msun.
△ Less
Submitted 27 February, 2020; v1 submitted 29 March, 2019;
originally announced April 2019.
-
Numerical studies on core collapse supernova in self-interacting massive scalar-tensor gravity
Authors:
Patrick Chi-Kit Cheong,
Tjonnie Guang Feng Li
Abstract:
We investigate stellar core collapse in scalar-tensor theory with a massive self-interacting scalar field. In these theories, strong long-lived inverse chirp signals could be induced during the stellar core collapse, which provides us with several potential smoking-gun signatures that could be found using current ground-based detectors. We show that the existence of self-interaction in the potenti…
▽ More
We investigate stellar core collapse in scalar-tensor theory with a massive self-interacting scalar field. In these theories, strong long-lived inverse chirp signals could be induced during the stellar core collapse, which provides us with several potential smoking-gun signatures that could be found using current ground-based detectors. We show that the existence of self-interaction in the potential of the scalar field can significantly suppress spontaneous scalarization and the amplitude of the monopole gravitational-wave radiation. Moreover, this suppression due to self-interaction is frequency dependent and may be discernible in LIGO/Virgo's sensitive band. Therefore, self-interaction should be considered when constraining scalar-tensor coupling parameters with gravitational-wave detections. Alternatively, if such monopole gravitational-wave signals are detected, then one may be able to infer the presence of self-interaction.
△ Less
Submitted 29 August, 2019; v1 submitted 12 December, 2018;
originally announced December 2018.
-
Commensurate magnetic structures ofRMn2O5 (R=Y, Ho, Bi) determined by Single-crystal Neutron Diffraction
Authors:
C. Vecchini L. C. Chapon,
P. J. Brown,
T. Chatterji,
S. Park S-W. Cheong,
P. G. Radaelli
Abstract:
Precise magnetic structures of RMn2O5, with R= Y, Ho, Bi in the commensurate/ferroelectric regime, have been determined by single-crystal neutron diffraction. For each system, the integrated intensities of a large number of independent magnetic Bragg reflections have been measured, allowing unconstrained least-squares refinement of the structures. The analysis confirms the previously reported ma…
▽ More
Precise magnetic structures of RMn2O5, with R= Y, Ho, Bi in the commensurate/ferroelectric regime, have been determined by single-crystal neutron diffraction. For each system, the integrated intensities of a large number of independent magnetic Bragg reflections have been measured, allowing unconstrained least-squares refinement of the structures. The analysis confirms the previously reported magnetic configuration in the ab-plane, in particular the existence of zig-zag antiferromagnetic chains. For the Y and Ho compounds additional weak magnetic components parallel to the c-axis were detected which are modulated in phase quadrature with the a-b components. This component is extremely small in the BiMn2O5 sample, therefore supporting symmetric exchange as the principal mechanism inducing ferroelectricity. For HoMn2O5, a magnetic ordering of the Ho moments was observed, which is consistent with a super-exchange interaction through the oxygens. For all three compounds, the point symmetry in the magnetically ordered state is m2m, allowing the polar b-axis found experimentally.
△ Less
Submitted 10 October, 2007;
originally announced October 2007.