Search

We have performed r-process calculations in neutron star mergers (NSM) and jets of magnetohydrodynamically driven (MHD) supernovae. In these very neutron-rich environments the fission model of heavy nuclei has an impact on the shape of the final abundance distribution and the second r-process peak in particular. We have studied the effect of different fission fragment mass distribution models in calculations of low-Yₑ ejecta, ranging from a simple parametrization to extensive statistical treatments (ABLA07). The r-process path ends when it reaches an area in the nuclear chart where fission dominates over further neutron captures. The position of this point is determined by the fission barriers and the neutron separation energies of the nuclei involved. As these values both depend on the choice of the nuclear mass model, so does the r-process path. Here we present calculations using the FRDM (Finite Range Droplet Model) and the ETFSI (Extended Thomas Fermi with Strutinsky Integral) mass model with the related TF and ETFSI fission barrier predictions. Utilizing sophisticated fission fragment distribution leads to a highly improved abundance distribution.

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

We have performed r-process calculations in neutron star mergers (NSM) and jets of magnetohydrodynamically driven (MHD) supernovae. In these very neutron-rich environments the fission model of heavy nuclei has an impact on the shape of the final abundance distribution and the second r-process peak in particular. We have studied the effect of different fission fragment mass distribution models in calculations of low-Yₑ ejecta, ranging from a simple parametrization to extensive statistical treatments (ABLA07). The r-process path ends when it reaches an area in the nuclear chart where fission dominates over further neutron captures. The position of this point is determined by the fission barriers and the neutron separation energies of the nuclei involved. As these values both depend on the choice of the nuclear mass model, so does the r-process path. Here we present calculations using the FRDM (Finite Range Droplet Model) and the ETFSI (Extended Thomas Fermi with Strutinsky Integral) mass model with the related TF and ETFSI fission barrier predictions. Utilizing sophisticated fission fragment distribution leads to a highly improved abundance distribution.

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here the current understanding of the mass ejection mechanisms, the properties of the ejecta, and the resulting radioactive transients are reviewed. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, questions where no consensus has been reached yet are also pointed out, and possible directions for the future research are sketched. In an appendix, a publicly available heating rate library based on the WinNet nuclear reaction network is described, and a simple fit formula to alleviate the implementation in hydrodynamic simulations is provided.

Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22248 high-fidelity simulations, we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multi-band inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle.

The interaction of $\beta$-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti (1971) following the seminal formulation of Bethe (1930), in a spectral solver of the Vlasov-Maxwell-Boltzmann equations. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon-matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric kilonova ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for $\beta$-particles in the kilonova are a bulk velocity parameter $\vec{u}=0$, a thermal velocity parameter $\vec{\alpha}=0.5c$, and a 9x9x9 mode velocity basis set (Hermite orders 0 to 8 in each dimension). For ejecta-interior sample zones, we estimate the ratio of thermalization from large-angle ($\gtrsim2.5^{\circ}$) bound excitation scattering to total thermalization is $\sim$0.002-0.003., Comment: Accepted to ApJ, 26 pages, 8 figures

In this study, we introduce a tensor-train (TT) finite difference WENO method for solving compressible Euler equations. In a step-by-step manner, the tensorization of the governing equations is demonstrated. We also introduce \emph{LF-cross} and \emph{WENO-cross} methods to compute numerical fluxes and the WENO reconstruction using the cross interpolation technique. A tensor-train approach is developed for boundary condition types commonly encountered in Computational Fluid Dynamics (CFD). The performance of the proposed WENO-TT solver is investigated in a rich set of numerical experiments. We demonstrate that the WENO-TT method achieves the theoretical $\text{5}^{\text{th}}$-order accuracy of the classical WENO scheme in smooth problems while successfully capturing complicated shock structures. In an effort to avoid the growth of TT ranks, we propose a dynamic method to estimate the TT approximation error that governs the ranks and overall truncation error of the WENO-TT scheme. Finally, we show that the traditional WENO scheme can be accelerated up to 1000 times in the TT format, and the memory requirements can be significantly decreased for low-rank problems, demonstrating the potential of tensor-train approach for future CFD application. This paper is the first study that develops a finite difference WENO scheme using the tensor-train approach for compressible flows. It is also the first comprehensive work that provides a detailed perspective into the relationship between rank, truncation error, and the TT approximation error for compressible WENO solvers.

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

We have performed r-process calculations in neutron star mergers (NSM) and jets of magnetohydrodynamically driven (MHD) supernovae. In these very neutron-rich environments the fission model of heavy nuclei has an impact on the shape of the final abundance distribution and the second r-process peak in particular. We have studied the effect of different fission fragment mass distribution models in calculations of low-Yₑ ejecta, ranging from a simple parametrization to extensive statistical treatments (ABLA07). The r-process path ends when it reaches an area in the nuclear chart where fission dominates over further neutron captures. The position of this point is determined by the fission barriers and the neutron separation energies of the nuclei involved. As these values both depend on the choice of the nuclear mass model, so does the r-process path. Here we present calculations using the FRDM (Finite Range Droplet Model) and the ETFSI (Extended Thomas Fermi with Strutinsky Integral) mass model with the related TF and ETFSI fission barrier predictions. Utilizing sophisticated fission fragment distribution leads to a highly improved abundance distribution.

We have performed r-process calculations in neutron star mergers (NSM) and jets of magnetohydrodynamically driven (MHD) supernovae. In these very neutron-rich environments the fission model of heavy nuclei has an impact on the shape of the final abundance distribution and the second r-process peak in particular. We have studied the effect of different fission fragment mass distribution models in calculations of low-Yₑ ejecta, ranging from a simple parametrization to extensive statistical treatments (ABLA07). The r-process path ends when it reaches an area in the nuclear chart where fission dominates over further neutron captures. The position of this point is determined by the fission barriers and the neutron separation energies of the nuclei involved. As these values both depend on the choice of the nuclear mass model, so does the r-process path. Here we present calculations using the FRDM (Finite Range Droplet Model) and the ETFSI (Extended Thomas Fermi with Strutinsky Integral) mass model with the related TF and ETFSI fission barrier predictions. Utilizing sophisticated fission fragment distribution leads to a highly improved abundance distribution.

We present a nucleosynthesis study on the neutrino-driven wind after neutron star merger. Post-processing tracers from a hydrodynamical simulation, we determine nucleosynthesis yields that depend on both life time of the massive neutron star and polar angle. In comparison with the dynamic ejecta we find a complementary nucleosynthesis, which can explain the r-process pattern beyond the first r-process peak. Additionally, we predict possible forms of light curves for the detection of a kilonova.

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models that are quick to evaluate and are calibrated with more detailed multi-physics simulations are needed to design observational strategies for kilonovae detection and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from binary neutron star (BNS) mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of BNS mergers produce peak i-band absolute magnitudes of −20 ≤ Mi ≤ −11. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution.

The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OT) discovered as gravitational wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-time direct comparison of models with our data. The principal facility is the Hobby-Eberly Telescope. LIGHETR uses the massively-replicated VIRUS array of spectrographs to search for associated OTs and obtain early blue spectra and in a complementary role, the low-resolution LRS-2 spectrograph is used to obtain spectra of viable candidates as well as a densely-sampled series of spectra of true counterparts. Once an OT is identified, the anticipated cadence of spectra would match or considerably exceed anything achieved for GW170817 = AT2017gfo for which there were no spectra in the first 12 hours and thereafter only roughly once daily. We describe special HET-specific software written to facilitate the program and attempts to determine the flux limits to undetected sources. We also describe our campaign to follow up OT candidates during the third observational campaign of the LIGO and Virgo Scientific Collaborations. We obtained VIRUS spectroscopy of candidate galaxy hosts for 5 LVC gravitational wave events and LRS-2 spectra of one candidate for the OT associated with S190901ap. We identified that candidate, ZTF19abvionh = AT2019pip, as a possible Wolf-Rayet star in an otherwise unrecognized nearby dwarf galaxy., Comment: 26 pages, 15 figures

We present the state-of-the-art single-zone nuclear reaction network WinNet that is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication, WinNet is publicly available and open source at GitHub and Zenodo.

In this work, we present a simple interpolation methodology for spectroscopic time series, based on conventional interpolation techniques (random forests) implemented in widely-available libraries. We demonstrate that our existing library of simulations is sufficient for training, producing interpolated spectra that respond sensitively to varied ejecta parameter, post-merger time, and viewing angle inputs. We compare our interpolated spectra to the AT2017gfo spectral data, and find parameters similar to our previous inferences using broadband light curves. However, the spectral observations have significant systematic short-wavelength residuals relative to our models, which we cannot explain within our existing framework. Similar to previous studies, we argue that an additional blue component is required. We consider a radioactive heating source as a third component characterized by light, slow-moving, lanthanide-free ejecta with $M_{\rm th} = 0.003~M_\odot$, $v_{\rm th} = 0.05$c, and $\kappa_{\rm th} = 1$ cm$^2$/g. When included as part of our radiative transfer simulations, our choice of third component reprocesses blue photons into lower energies, having the opposite effect and further accentuating the blue-underluminosity disparity in our simulations. As such, we are unable to overcome short-wavelength deficits at later times using an additional radioactive heating component, indicating the need for a more sophisticated modeling treatment., Comment: 11 pages, 7 figures, submitted to Phys. Rev. Research, presented at April APS session F13.00006

The electromagnetic emission from the non-relativistic ejecta launched in neutron star mergers (either dynamically or through a disk wind) has the potential to probe both the total mass and composition of this ejecta. These observations are crucial in understanding the role of these mergers in the production of r-process elements in the universe. However, many properties of the ejecta can alter the light-curves and we must both identify which properties play a role in shaping this emission and understand the effects these properties have on the emission before we can use observations to place strong constraints on the amount of r-process elements produced in the merger. This paper focuses on understanding the effect of the velocity distribution (amount of mass moving at different velocities) for lanthanide-rich ejecta on the light-curves and spectra. The simulations use distributions guided by recent calculations of disk outflows and compare the velocity-distribution effects to those of ejecta mass, velocity and composition. Our comparisons show that uncertainties in the velocity distribution can lead to factor of 2-4 uncertainties in the inferred ejecta mass based on peak infra-red luminosities. We also show that early-time UV or optical observations may be able to constrain the velocity distribution, reducing the uncertainty in the ejecta mass., Comment: 21 pages, 16 figures, submitted to ApJ

We investigate a variety of short warning time, terminal mitigation scenarios via fragmentation for a hypothetical impact of asteroid 2023 NT1, a Near-Earth Object (NEO) that was discovered on July 15, 2023, two days after its closest approach to Earth on July 13. The asteroid passed by Earth within ~0.25 lunar distances with a closest approach distance of ~10$^{5}$ km and speed of 11.27 km/s. Its size remains largely uncertain, with an estimated diameter range of 26 - 58 m and probable diameter estimate (weighted by the NEO size frequency distribution) of 34 m (JPL Sentry, September 12, 2023). The asteroid approached Earth from the direction of the Sun, as did both the Chelyabinsk asteroid in 2013 and comet NEOWISE in 2021. As a result, 2023 NT1 remained undetected until after its closest approach. If it had been on a collision course, it would have had an impact energy of ~1.5 Mt (assuming a spherical asteroid with the probable diameter estimate of 34 m, 2.6 g/cm$^{3}$ uniform density, and impact speed of 15.59 km/s). 2023 NT1 represents a threat that could have caused significant local damage (~3x Chelyabinsk airburst energy). We utilize the PI ("Pulverize It") method for planetary defense to model potential mitigation scenarios of an object like 2023 NT1 through simulations of hypervelocity asteroid disruption and atmospheric ground effects for the case of a terminal defense mode. Simulations suggest that PI is an effective multimodal approach for planetary defense that can operate in extremely short interdiction modes (with intercepts as short as hours prior to impact), in addition to long interdiction time scales with months to years of warning. Our simulations support the proposition that threats like 2023 NT1 can be effectively mitigated with intercepts of one day (or less) prior to impact, yielding minimal to no ground damage, using modest resources and existing technologies.

Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements., Comment: Presented at the Scientific Visualization Supercomputing Conference 2021, St. Louis, MO

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of Athena with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs).

Studying compact star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, Smoothed Particle Hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, Halted-Pendulum Relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

The electromagnetic emission resulting from neutron star mergers have been shown to encode properties of the ejected material in their light curves. The ejecta properties inferred from the kilonova emission has been in tension with those calculated based on the gravitational wave signal and numerical relativity models. Motivated by this tension, we construct a broad set of surrogate light curve models derived for kilonova ejecta. The four-parameter family of two-dimensional anisotropic simulations and its associated surrogate explore different assumptions about the wind outflow morphology and outflow composition, keeping the dynamical ejecta component consistent. We present the capabilities of these surrogate models in interpolating kilonova light curves across various ejecta parameters and perform parameter estimation for AT2017gfo both without any assumptions on the outflow and under the assumption that the outflow must be representative of solar r-process abundance patterns. Our parameter estimation for AT2017gfo shows these surrogate models help alleviate the ejecta property discrepancy while also illustrating the impact of systematic modeling uncertainties on these properties, urging further investigation., Comment: 15 pages, 6 figures, data available in Zenodo (https://zenodo.org/record/7335961) and GitHub (https://github.com/markoris/surrogate_kne)

Smoothed Particle Hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron-star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches which can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations., Comment: 28 pages, 34 figures

Compact binary mergers involving neutron stars can eject a fraction of their mass to space. Being extremely neutron rich, this material undergoes rapid neutron capture nucleosynthesis, and the resulting radioactivity powers fast, short-lived electromagnetic transients known as kilonova or macronova. Such transients are exciting probes of the most extreme physical conditions and their observation signals the enrichment of the Universe with heavy elements. Here we review our current understanding of the mass ejection mechanisms, the properties of the ejecta and the resulting radioactive transients. The first well-observed event in the aftermath of GW170817 delivered a wealth of insights, but much of today's picture of such events is still based on a patchwork of theoretical studies. Apart from summarizing the current understanding, we also point out questions where no consensus has been reached yet, and we sketch possible directions for the future research. In an appendix, we describe a publicly available heating rate library based on the WinNet nuclear reaction network, and we provide a simple fit formula to alleviate the implementation in hydrodynamic simulations., Comment: 22 pages, updated according to the suggestions of two referees; accepted by Annalen der Physik

Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming the observed source is prototypical, this inferred abundance pattern in turn must match $r$-process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron-star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta ($Y_{\rm{e}} = 0.035$) along with moderately neutron-rich ($Y_{\rm e} = 0.27$) wind ejecta composition yields a wind-to-dynamical mass ratio of $M_{\rm{w}}/M_{\rm{d}}$ = 0.47 which best matches the observed AT2017gfo kilonova light curves while also producing the best-matching abundance of neutron-capture elements in the solar system., Comment: 16 pages, 9 figures, submitted to PRD

The 2017 detection of a kilonova coincident with gravitational-wave emission has identified neutron star mergers as the major source of the heaviest elements, and dramatically constrained alternative theories of gravity. Observing a population of such sources has the potential to transform cosmology, nuclear physics, and astrophysics. However, with only one confident multi-messenger detection currently available, modelling the diversity of signals expected from such a population requires improved theoretical understanding. In particular, models which are quick to evaluate, and are calibrated with more detailed multi-physics simulations, are needed to design observational strategies for kilonovae detection, and to obtain rapid-response interpretations of new observations. We use grey-opacity models to construct populations of kilonovae, spanning ejecta parameters predicted by numerical simulations. Our modelling focuses on wavelengths relevant for upcoming optical surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). In these simulations, we implement heating rates that are based on nuclear reaction network calculations. We create a Gaussian-process emulator for kilonova grey opacities, calibrated with detailed radiative transfer simulations. Using recent fits to numerical relativity simulations, we predict how the ejecta parameters from BNS mergers shape the population of kilonovae, accounting for the viewing-angle dependence. Our simulated population of binary neutron star (BNS) mergers produce peak i-band absolute magnitudes $-20 \leq M_i \leq -11$. A comparison with detailed radiative transfer calculations indicates that further improvements are needed to accurately reproduce spectral shapes over the full light curve evolution., Comment: 14 pages, 12 figures, 2 tables, matches version accepted to MNRAS. Fixes a software bug from previous versions which incorrectly determined the kilonovae ejecta velocities, v4: added link to public code package https://github.com/cnsetzer/kilopop

Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements., Comment: Presented at the Scientific Visualization Supercomputing Conference 2021, St. Louis, MO

A domain decomposition approach is proposed to study the problem of a high-speed vessel which moves in close proximity to the edge of an ice field. A two-dimensional boundary element method which incorporates fully non-linear free surface conditions shall be used in combination with 2D+t theory to model the flow in the open water region. Theory of linear hydro-elasticity shall be used in the ice covered region with the ice modelled as semi-infinite, thin and elastic plate of constant thickness. A non-local boundary is specified on a fictitious boundary between both regions to account for interaction effects., US Office of Naval Research

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of \textsl{Athena} with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs)., Comment: White paper of the multi-messenger Athena synergy excersise

Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best knowledge of nuclear matter in the form of an equation of state as well as effects of general relativity (GR). However, one important aspect of neutron stars is usually ignored: the solid nature of their crust. The solid matter in the crust is the strongest material known in nature which could lead to a multitude of possible observational effects that have not been studied in dynamical simulations yet. The crust could change the way a neutron star deforms during a merger, leaving an imprint in the gravitational wave signal. It could even shatter during the inspiral, producing a potentially observable electromagnetic signal. Here, we present a first study of the dynamical behavior of neutron stars with a solid crust and fixed GR background with FleCSPH. FleCSPH is a general-purpose SPH code, developed at Los Alamos National Laboratory. It features an efficient algorithm for gravitational interactions via the Fast Multipole Method, which, together with the implemented nuclear equation of state, makes it appropriate for astrophysical applications. The solid material dynamics is described via the elastic-perfectly plastic model with maximum-strain breaking. Despite its simplicity, the model reproduces the stress-strain curve of crustal material as extracted from microphysical simulations very well. We present first tests of our implementation via simulations of neutron star oscillations and give an outlook on our study of the dynamical behavior of the solid crust in neutron star merger events., Comment: 9 pages, 9 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-12 2021

Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron stars, black holes, and white dwarfs). However, in astrophysical applications that require the inclusion of gravity, calculating pairwise particle interactions becomes prohibitively expensive. In the Fast Multipole Method (FMM), they are, therefore, replaced with symmetric interactions between distant clusters of particles (contained in the tree nodes) Although such an algorithm is linear momentum-conserving, it introduces spurious torques that violate conservation of angular momentum. We present a modification of FMM that is free of spurious torques and conserves angular momentum explicitly. The new method has practically no computational overhead compared to the standard FMM., Comment: 6 pages, 3 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-11 2021

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M$_{\odot}$ and 0.05 to 0.3$c$, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction ($Y_e$) ejecta with a robust r-process composition and either a spherical or lobed high-$Y_e$ ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817., Comment: 19 pages, 11 figures, submitted to ApJ

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters., Comment: 22 pages, 8 figures, 1 appendix (with 9 pages, 10 figures), submitted to ApJ

Detailed radiative transfer simulations of kilonovae are difficult to apply directly to observations; they only sparsely cover simulation parameters, such as the mass, velocity, morphology, and composition of the ejecta. On the other hand, semianalytic models for kilonovae can be evaluated continuously over model parameters, but neglect important physical details which are not incorporated in the simulations, thus introducing systematic bias. Starting with a grid of 2D anisotropic simulations of kilonova light curves covering a wide range of ejecta properties, we apply adaptive-learning techniques to iteratively choose new simulations and produce high-fidelity surrogate models for those simulations. These surrogate models allow for continuous evaluation across model parameters while retaining the microphysical details about the ejecta. Using a new code for multimessenger inference, we demonstrate how to use our interpolated models to infer kilonova parameters. Comparing to inferences using simplified analytic models, we recover different ejecta properties. We discuss the implications of this analysis which is qualitatively consistent with similar previous work using detailed ejecta opacity calculations and which illustrates systematic challenges for kilonova modeling. An associated data and code release provides our interpolated light-curve models, interpolation implementation which can be applied to reproduce our work or extend to new models, and our multimessenger parameter inference engine., Comment: 18 pages, 15 figures, submitting to Physical Review Research

This groundbreaking report exposes how Border Patrol, an agency within U.S. Customs and Border Protection (CBP), uses racial profiling to target immigrants from Latin America and other people of color throughout Michigan. The report also reveals how Border Patrol colludes with state and local police agencies to target, arrest, and deport immigrants, many of whom are longtime Michigan residents.

The detailed observations of GW170817 proved for the first time directly that neutron star mergers are a major production site of heavy elements. The observations could be fit by a number of simulations that qualitatively agree, but can quantitatively differ (e.g., in total r-process mass) by an order of magnitude. We categorize kilonova ejecta into several typical morphologies motivated by numerical simulations, and apply a radiative transfer Monte Carlo code to study how the geometric distribution of the ejecta shapes the emitted radiation. We find major impacts on both spectra and light curves. The peak bolometric luminosity can vary by two orders of magnitude and the timing of its peak by a factor of five. These findings provide the crucial implication that the ejecta masses inferred from observations around the peak brightness are uncertain by at least an order of magnitude. Mixed two-component models with lanthanide-rich ejecta are particularly sensitive to geometric distribution. A subset of mixed models shows very strong viewing angle dependence due to lanthanide curtaining, which persists even if the relative mass of lanthanide-rich component is small. The angular dependence is weak in the rest of our models, but different geometric combinations of the two components lead to a highly diverse set of light curves. We identify geometry-dependent P Cygni features in late spectra that directly map out strong lines in the simulated opacity of neodymium, which can help to constrain the ejecta geometry and to directly probe the r-process abundances.

Radiography is often used to probe complex, evolving density fields in dynamic systems and in so doing gain insight into the underlying physics. This technique has been used in numerous fields including materials science, shock physics, inertial confinement fusion, and other national security applications. In many of these applications, however, complications resulting from noise, scatter, complex beam dynamics, etc. prevent the reconstruction of density from being accurate enough to identify the underlying physics with sufficient confidence. As such, density reconstruction from static/dynamic radiography has typically been limited to identifying discontinuous features such as cracks and voids in a number of these applications. In this work, we propose a fundamentally new approach to reconstructing density from a temporal sequence of radiographic images. Using only the robust features identifiable in radiographs, we combine them with the underlying hydrodynamic equations of motion using a machine learning approach, namely, conditional generative adversarial networks (cGAN), to determine the density fields from a dynamic sequence of radiographs. Next, we seek to further enhance the hydrodynamic consistency of the ML-based density reconstruction through a process of parameter estimation and projection onto a hydrodynamic manifold. In this context, we note that the distance from the hydrodynamic manifold given by the training data to the test data in the parameter space considered both serves as a diagnostic of the robustness of the predictions and serves to augment the training database, with the expectation that the latter will further reduce future density reconstruction errors. Finally, we demonstrate the ability of this method to outperform a traditional radiographic reconstruction in capturing allowable hydrodynamic paths even when relatively small amounts of scatter are present., Comment: Submitted to Optics Express

In this paper we explore the scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, and CTA for very high energy observations. These science themes encompass pressing issues in astrophysics, cosmology and fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in Super-Massive Binary Black Holes (SMBBHs) and in compact stellar binaries, the equation of state of neutron stars, cosmic accelerators and the origin of Cosmic Rays (CRs), the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and the Standard Model. Observational strategies for implementing the identified science topics are also discussed. A significant part of the sources targeted by multi-messenger facilities is of transient nature. We have thus also discussed the synergy of \textsl{Athena} with wide-field high-energy facilities, taking THESEUS as a case study for transient discovery. This discussion covers all the Athena science goals that rely on follow-up observations of high-energy transients identified by external observatories, and includes also topics that are not based on multi-messenger observations, such as the search for missing baryons or the observation of early star populations and metal enrichment at the cosmic dawn with Gamma-Ray Bursts (GRBs)., Comment: White paper of the multi-messenger Athena synergy excersise

Measurements of the gravitational-wave signals from neutron star mergers allow scientists to learn about the interior of neutron stars and the properties of dense nuclear matter. The study of neutron star mergers is usually performed with computational fluid dynamics codes, mostly in Eulerian but also in Lagrangian formulation such as smoothed particle hydrodynamics (SPH). Codes include our best knowledge of nuclear matter in the form of an equation of state as well as effects of general relativity (GR). However, one important aspect of neutron stars is usually ignored: the solid nature of their crust. The solid matter in the crust is the strongest material known in nature which could lead to a multitude of possible observational effects that have not been studied in dynamical simulations yet. The crust could change the way a neutron star deforms during a merger, leaving an imprint in the gravitational wave signal. It could even shatter during the inspiral, producing a potentially observable electromagnetic signal. Here, we present a first study of the dynamical behavior of neutron stars with a solid crust and fixed GR background with FleCSPH. FleCSPH is a general-purpose SPH code, developed at Los Alamos National Laboratory. It features an efficient algorithm for gravitational interactions via the Fast Multipole Method, which, together with the implemented nuclear equation of state, makes it appropriate for astrophysical applications. The solid material dynamics is described via the elastic-perfectly plastic model with maximum-strain breaking. Despite its simplicity, the model reproduces the stress-strain curve of crustal material as extracted from microphysical simulations very well. We present first tests of our implementation via simulations of neutron star oscillations and give an outlook on our study of the dynamical behavior of the solid crust in neutron star merger events., Comment: 9 pages, 9 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-12 2021

Smoothed particle hydrodynamics (SPH) is positioned as having ideal conservation properties. When properly implemented, conservation of total mass, energy, and both linear and angular momentum is guaranteed exactly, up to machine precision. This is particularly important for some applications in computational astrophysics, such as binary dynamics, mergers, and accretion of compact objects (neutron stars, black holes, and white dwarfs). However, in astrophysical applications that require the inclusion of gravity, calculating pairwise particle interactions becomes prohibitively expensive. In the Fast Multipole Method (FMM), they are, therefore, replaced with symmetric interactions between distant clusters of particles (contained in the tree nodes) Although such an algorithm is linear momentum-conserving, it introduces spurious torques that violate conservation of angular momentum. We present a modification of FMM that is free of spurious torques and conserves angular momentum explicitly. The new method has practically no computational overhead compared to the standard FMM., Comment: 6 pages, 3 figures, proceedings of 2021 international SPHERIC workshop (virtual), June, 8-11 2021

Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of 2-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1 M$_{\odot}$ and 0.05 to 0.3$c$, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction ($Y_e$) ejecta with a robust r-process composition and either a spherical or lobed high-$Y_e$ ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code SuperNu using a full suite of lanthanide and 4th row element opacities. We examine the trends of these models with parameter variation, show how it can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817., Comment: 19 pages, 11 figures, submitted to ApJ

Kilonovae are ultraviolet, optical, and infrared transients powered by the radioactive decay of heavy elements following a neutron star merger. Joint observations of kilonovae and gravitational waves can offer key constraints on the source of Galactic $r$-process enrichment, among other astrophysical topics. However, robust constraints on heavy element production requires rapid kilonova detection (within $\sim 1$ day of merger) as well as multi-wavelength observations across multiple epochs. In this study, we quantify the ability of 13 wide field-of-view instruments to detect kilonovae, leveraging a large grid of over 900 radiative transfer simulations with 54 viewing angles per simulation. We consider both current and upcoming instruments, collectively spanning the full kilonova spectrum. The Roman Space Telescope has the highest redshift reach of any instrument in the study, observing kilonovae out to $z \sim 1$ within the first day post-merger. We demonstrate that BlackGEM, DECam, GOTO, the Vera C. Rubin Observatory's LSST, ULTRASAT, and VISTA can observe some kilonovae out to $z \sim 0.1$ ($\sim$475 Mpc), while DDOTI, MeerLICHT, PRIME, $Swift$/UVOT, and ZTF are confined to more nearby observations. Furthermore, we provide a framework to infer kilonova ejecta properties following non-detections and explore variation in detectability with these ejecta parameters., Comment: 22 pages, 8 figures, 1 appendix (with 9 pages, 10 figures), submitted to ApJ

Detailed radiative transfer simulations of kilonovae are difficult to apply directly to observations; they only sparsely cover simulation parameters, such as the mass, velocity, morphology, and composition of the ejecta. On the other hand, semianalytic models for kilonovae can be evaluated continuously over model parameters, but neglect important physical details which are not incorporated in the simulations, thus introducing systematic bias. Starting with a grid of 2D anisotropic simulations of kilonova light curves covering a wide range of ejecta properties, we apply adaptive-learning techniques to iteratively choose new simulations and produce high-fidelity surrogate models for those simulations. These surrogate models allow for continuous evaluation across model parameters while retaining the microphysical details about the ejecta. Using a new code for multimessenger inference, we demonstrate how to use our interpolated models to infer kilonova parameters. Comparing to inferences using simplified analytic models, we recover different ejecta properties. We discuss the implications of this analysis which is qualitatively consistent with similar previous work using detailed ejecta opacity calculations and which illustrates systematic challenges for kilonova modeling. An associated data and code release provides our interpolated light-curve models, interpolation implementation which can be applied to reproduce our work or extend to new models, and our multimessenger parameter inference engine., Comment: 19 pages, 15 figures

This groundbreaking report exposes how Border Patrol, an agency within U.S. Customs and Border Protection (CBP), uses racial profiling to target immigrants from Latin America and other people of color throughout Michigan. The report also reveals how Border Patrol colludes with state and local police agencies to target, arrest, and deport immigrants, many of whom are longtime Michigan residents.

This groundbreaking report exposes how Border Patrol, an agency within U.S. Customs and Border Protection (CBP), uses racial profiling to target immigrants from Latin America and other people of color throughout Michigan. The report also reveals how Border Patrol colludes with state and local police agencies to target, arrest, and deport immigrants, many of whom are longtime Michigan residents.

Dip coating is a common technique used to cover a solid surface with a thin liquid film, the thickness of which was successfully predicted by the theory developed in the 1940s by Landau & Levich (Acta Physicochem. URSS, vol. 17, 1942, pp. 141-153) and Derjaguin (Acta Physicochem. URSS, vol. 20, 1943, pp. 349-352). In this work, we present an extension of their theory to the case where the dipping bath contains two immiscible liquids, one lighter than the other, resulting in the entrainment of two thin films on the substrate. We report how the thicknesses of the coated films depend on the capillary number, on the ratios of the properties of the two liquids and on the relative thickness of the upper fluid layer in the bath. We also show that the liquid/liquid and liquid/gas interfaces evolve independently from each other as if only one liquid were coated, except for a very small region where their separation falls quickly to its asymptotic value and the shear stresses at the two interfaces peak. Interestingly, we find that the final coated thicknesses are determined by the values of these maximum shear stresses., SCOPUS: ar.j, info:eu-repo/semantics/published

This groundbreaking report exposes how Border Patrol, an agency within U.S. Customs and Border Protection (CBP), uses racial profiling to target immigrants from Latin America and other people of color throughout Michigan. The report also reveals how Border Patrol colludes with state and local police agencies to target, arrest, and deport immigrants, many of whom are longtime Michigan residents.