Engineering topics

Chen, Che-Yu

Publications and source records attributed to Chen, Che-Yu.

11 records

DOE OSTI2024

Modeling the Far-infrared Polarization Spectrum of a High-mass Star-forming Cloud

The polarization spectrum, or wavelength dependence of the polarization fraction, of interstellar dust emission provides important insights into the grain alignment mechanism of interstellar dust grains. We investigate the far-infrared polarization spectrum of a realistic simulated high-mass star-forming cloud under various models of grain alignment and emission. We find that neither a homogeneous grain alignment model nor a grain alignment model that includes collisional dealignment is able to produce the falling spectrum seen in observations. On the other hand, we find that a grain alignment model with grain alignment efficiency dependent on local temperature is capable of producing a falling spectrum that is in qualitative agreement with observations of OMC-1. For the model most in agreement with OMC-1, we find no correlation between the temperature and the slope of the polarization spectrum. However, we do find a positive correlation between the column density and the slope of the polarization spectrum. We suggest this latter correlation to be the result of wavelength-dependent polarization by absorption.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Dynamics in Star-forming Cores (DiSCo): project overview and the first look towards the B1 and NGC 1333 regions in Perseus

The internal velocity structure within dense gaseous cores plays a crucial role in providing the initial conditions for star formation in molecular clouds. However, the kinematic properties of dense gas at core scales (~0.01–0.1 pc) has not been extensively characterized because of instrument limitations until the unique capabilities of GBT-Argus became available. The ongoing GBT-Argus Large Program, Dynamics in Star-forming Cores (DiSCo) thus aims to investigate the origin and distribution of angular momentum of star-forming cores. DiSCo will survey all starless cores and Class 0 protostellar cores in the Perseus molecular complex down to ~0.01 pc scales with <0.05 km s –1 velocity resolution using the dense gas tracer N 2 H + . Here, we present the first data sets from DiSCo towards the B1 and NGC 1333 regions in Perseus. Our results suggest that a dense core’s internal velocity structure has little correlation with other core-scale properties, indicating these gas motions may be originated externally from cloud-scale turbulence. These first data sets also reaffirm the ability of GBT-Argus for studying dense core velocity structure and provided an empirical basis for future studies that address the angular momentum problem with a statistically broad sample.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Modelling CN Zeeman effect observations of the envelopes of a low-mass protostellar disc and a massive protostar

We use the POLARIS radiative transfer code to produce simulated circular polarization Zeeman emission maps of the cyanide (CN) J = 1–0 molecular line transition for two types of protostellar envelope magnetohydrodynamic simulations. Our first model is a low-mass disc envelope system (box length L = 200 au), and our second model is the envelope of a massive protostar (L = 104 au) with a protostellar wind and a CN-enhanced outflow shell. We compute the velocity-integrated Stokes I and V, as well as the implied V/I polarization percentage, for each detector pixel location in our simulated emission maps. Our results show that both types of protostellar environments are in principle accessible with current circular polarization instruments, with each containing swaths of envelope area that yield percentage polarizations that exceed the 1.8 percent nominal sensitivity limit for circular polarization experiments with the Atacama Large Millimeter/submillimeter Array. In both systems, high-polarization (≳1.8 percent) pixels tend to lie at an intermediate distance away from the central star and where the line-centre opacity of the CN emission is moderately optically thin (τLC ~ 0.1–1). Furthermore, our computed V/I values scale roughly with the density-weighted mean line-of-sight magnetic field strength, indicating that Zeeman observations can effectively diagnose the strength of envelope-scale magnetic fields. We also find that pixels with large V/I are preferentially co-located where the absolute value of the velocity-integrated V is also large, suggesting that locations with favourable percentage polarization are also favourable in terms of raw signal.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Relative alignment between magnetic fields and molecular gas structure in molecular clouds

Here, we compare the structure of synthetic dust polarization with synthetic molecular line emission from radiative transfer calculations using a three-dimensional, turbulent collapsing-cloud magnetohydrodynamics simulation. The histogram of relative orientation (HRO) technique and the projected Rayleigh statistic (PRS) are considered. In our trans-Alfvénic (more strongly magnetized) simulation, there is a transition to perpendicular alignment at densities above ~4 × 10 3 cm –3 . This transition is recovered in most of our synthetic observations of optically thin molecular tracers; however, for 12 CO it does not occur and the PRS remains in parallel alignment across the whole observer space. We calculate the physical depth of the optical depth τ = 1 surface and find that for 12 CO it is largely located in front of the cloud mid-plane, suggesting that 12 CO is too optically thick and instead mainly probes low-volume density gas. In our super-Alfvénic simulation, the magnetic field becomes significantly more tangled, and all observed tracers tend towards no preference for perpendicular or parallel alignment. An observable difference in alignment between optically thin and optically thick tracers may indicate the presence of a dynamically important magnetic field, though there is some degeneracy with viewing angle. We convolve our data with a Gaussian beam and compare it with HRO results of the Vela C molecular cloud. We find good agreement between these results and our sub-Alfvénic simulations when viewed with the magnetic field in the plane of the sky (especially when sensitivity limitations are considered), though the observations are also consistent with an intermediately inclined magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Revealing the Drag Instability in One-fluid Nonideal Magnetohydrodynamic Simulations of a 1D Isothermal C-shock

C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in nonideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the post-shock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Generating rotating spacetime in Ricci-based gravity: naked singularity as a black hole mimicker

Motivated by the lack of rotating solutions sourced by matter in General Relativity as well as in modified gravity theories, in this work we extend a recently discovered exact rotating solution of the minimal Einstein-scalar theory to its counterpart in Eddington-inspired Born-Infeld gravity coupled to a Born-Infeld scalar field. This is accomplished with the implementation of a well-developed mapping between solutions of Ricci-Based Palatini theories of gravity and General Relativity. The new solution is parametrized by the scalar charge and the Born-Infeld coupling constant apart from the mass and spin of the compact object. Compared to the spacetime prior to the mapping, we find that the high-energy modifications at the Born-Infeld scale are able to suppress but not remove the curvature divergence of the original naked null singularity. Depending on the sign of the Born-Infeld coupling constant, these modifications may even give rise to an additional timelike singularity exterior to the null one. In spite of that, both of the naked singularities before and after the mapping are capable of casting shadows, and as a consequence of the mapping relation, their shadows turn out to be identical as seen by a distant observer on the equatorial plane. Even though the scalar field induces a peculiar oblateness to the appearance of the shadow with its left and right endpoints held fixed, the closedness condition for the shadow contour sets a small upper bound on the absolute value of the scalar charge, which leads to observational features of the shadow closely resembling those of a Kerr black hole.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Lessons from black hole quasinormal modes in modified gravity

Quasinormal modes (QNMs) of perturbed black holes have recently gained much interest because of their tight relations with the gravitational wave signals emitted during the post-merger phase of a binary black hole coalescence. One of the intriguing features of these modes is that they respect the no-hair theorem, and hence, they can be used to test black hole spacetimes and the underlying gravitational theory. In this paper, we exhibit three different aspects of how black hole QNMs could be altered in theories beyond Einstein’s general relativity (GR). These aspects are (i) the direct alterations of QNM spectra as compared with those in GR, (ii) the violation of the geometric correspondence between the high-frequency QNMs and the photon geodesics around the black hole, and (iii) the breaking of the isospectrality between the axial and polar gravitational perturbations. Several examples will be provided in each individual case. The prospects and possible challenges associated with future observations will be also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Black hole perturbations and quasinormal modes in hybrid metric-Palatini gravity

The rapid advancement of gravitational wave astronomy in recent years has paved the way for the burgeoning development of black hole spectroscopy, which enhances the possibility of testing black holes by their quasinormal modes (QNMs). In this paper, the axial gravitational perturbations and the QNM frequencies of black holes in the hybrid metric-Palatini gravity (HMPG) are investigated. Additionally, the HMPG theory is characterized by a dynamical scalar degree of freedom and is able to explain the late-time accelerating expansion of the universe without introducing any ad hoc screening mechanism to preserve the dynamics at the Solar System scale. We obtain the master equation governing the axial gravitational perturbations of the HMPG black holes and calculate the QNM frequencies. Moreover, in the scrutiny of the black holes and their QNMs, we take into account the constraints on the model parameters based on the post-Newtonian analysis, and show how the QNM frequencies of the HMPG black holes would be altered in the observationally consistent range of parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Asymptotic non-flatness of an effective black hole model based on loop quantum gravity

For any non-rotating effective quantum (uncharged) black hole model to be viable, its asymptotic structure of spacetime should reduce to that of a Schwarzschild black hole. After examining the asymptotic structure of quantum black holes proposed in Ashtekar et al., (2018), in this study we find that the solution is actually not asymptotically flat. Departure from asymptotically flat spacetime is a direct consequence of the quantum parameters included in the model. Moreover, since the property of asymptotic flatness needs to be explicitly invoked, for instance, in relating the black hole and the white hole masses, it thereby raises a serious concern about the consistency of this model in view of the classical limit.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Annihilation-to-nothing: a quantum gravitational boundary condition for the Schwarzschild black hole

The interior of a static Schwarzschild metric can be written in terms of two functions, similar to some models of anisotropic cosmology. With a suitable choice of canonical variables, we solve the Wheeler-DeWitt equation (WDW) inside the horizon of a Schwarzschild black hole. By imposing classicality near the horizon, and requiring boundedness of the wave function, we get a rather generic solution of the WDW equation, whose steepest-descent solution, i.e., the ridge of the wave function, coincides nicely with the classical trajectory. However, there is an ambiguity in defining the arrow of time which leads to two possible interpretations—(i) if there is only one arrow of time, one can infer that the steepest-descent of the wave function follows the classical trajectory throughout: coming from the event horizon and going all the way down to the singularity, while (ii) if there are two different arrows of time in two separate regimes, it can be inferred that the steepest-descent of the wave function comes inwards from the event horizon in one region while it moves outwards from the singularity in the other region, and there exists an annihilation process of these two parts of the wave function inside the horizon. Additionally, adopting the second interpretation could shed light on the information loss paradox: as time goes on, probabilities for histories that include black holes and singularities decay to zero and eventually only trivial geometries dominate.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Eikonal black hole ringings in generalized energy-momentum squared gravity

In the scope of black hole spectroscopy, several attempts have been made in the past decades to test black holes or gravitational theories via black hole quasinormal modes. In the eikonal approximation, the quasinormal modes are generically associated with the photon ring of the black hole. This correspondence is valid for most cases in general relativity, but may not be true in other theories of gravity. In this paper, we consider the generalized energy-momentum squared gravity in which matter fields are nonminimally coupled to geometry. We investigate the axial perturbations of the charged black holes in this model, without assuming any explicit expression of the action functional. After obtaining the modified Klein-Gordon equation and the modified Maxwell equations, we perturb the gravitational equations and the modified Maxwell equations to derive the coupled master equations of the axial perturbations. In the presence of the nonminimal coupling between matter and geometry, the correspondence between the eikonal quasinormal modes and the photon ring is not satisfied in general. Also, the two coupled fields of the axial perturbations are found to propagate independently, and they do not share the same quasinormal frequencies in the eikonal limit.

79 ASTRONOMY AND ASTROPHYSICS↗