In recent years the Scalar Field Dark Matter (SFDM), also called ultralight bosonic dark matter, has received considerable attention due to the number of problems it might help to solve. Among these we find the cusp-core problem and the abundance of small structures of the standard Cold Dark Matter (CDM) model. In this letter we show that mixed state solutions of the low energy and weak gravitational field limit of field equations, interpreted as galactic halo density profiles, could explain the anisotropic distribution of satellite galaxies observed in the Milky Way, M31 and Centaurus A, where satellites seem to concentrate more towards the poles of the galaxies instead of following homogeneously distributed trajectories. The core hypothesis is that multi-state solutions of the equations describing the dynamics of this dark matter candidate, namely, the Gross-Pitaevskii-Poisson equations, with monopolar and dipolar contributions, can explain the anisotropy of satellite trajectories.

Using numerical simulations, we study the effects of magnetic resistivity and thermal conductivity in the dynamics and properties of solar jets with characteristics of Type II spicules and cool coronal jets. The dynamic evolution of the jets is governed by the resistive MHD equations with thermal conduction along the magnetic field lines on a 2.5D slice. The magnetic field configuration consists of two symmetric neighboring loops with opposite polarity, used to support reconnection and followed by the plasma jet formation. In total 10 simulations were carried out with different values of resistivity and thermal conductivity, that produce jets with different morphological and thermal properties we quantify. We find that an increase in magnetic resistivity does not produce significant effects on the morphology, velocity and temperature of the jets. However, thermal conductivity affects both temperature and morphology of the jets. In particular, thermal conductivity causes jets to reach greater heights and increases the temperature of the jet-apex. Also, heat flux maps indicate the jet-apex and corona interchange energy more efficiently than the body of the jet. These results could potentially open a new avenue for plasma diagnostics in the Sun's atmosphere.

We present a general strategy to solve the stationary Schr\"odinger-Poisson (SP) system of equations for multistates with axial symmetry. The approach allows us to obtain the well known single and multistate solutions with spherical symmetry, Newtonian multistate $\ell-$boson stars and axially symmetric multistate configurations. For each case we construct particular examples that illustrate the method, whose stability properties are studied by numerically solving the time-dependent SP system. Among the stable configurations there are the mixed-two-state configurations including spherical and dipolar components, which might have an important value as potential anisotropic dark matter halos in the context of ultralight bosonic dark matter scenarios. This is the reason why we also present a possible process of formation of these mixed-two-state configurations that could open the door to the exploration of more general multistate structure formation scenarios.

The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization.

At intermediate to high densities, electron (de-)excitation collisions are the dominant process for populating or depopulating high Rydberg states. In particular, the accurate knowledge of the energy changing ($n$-changing) collisional rates is determinant for predicting the radio recombination spectra of gaseous nebula. The different datasets present in the literature come either from impact parameter calculations or semi-empirical fits and the rate coefficients agree within a factor of two. We show in this paper that these uncertainties cause errors lower than 5% in the emission of radio recombination lines (RRL) of most ionized plasmas of typical nebulae. However, in special circumstances where the transitions between Rydberg levels are amplified by maser effects, the errors can increase up to 20%. We present simulations of the optical depth and H$n\alpha$ line emission of Active Galactic Nuclei (AGN) Broad Line Regions (BLRs) and the Orion Nebula Blister to showcase our findings.

In this work we study some features of head-on mergers of equilibrium solutions of the Gross-Pitaevskii-Poisson system that rules the dynamics of the ultralight bosonic dark matter model. The importance of equilibrium solutions is that they play the role of halo cores in structure formation simulations. We consider a given range of initial conditions in order to sample the parameter space in terms of mass ratio and head-on momentum. In each case we analyze the relaxation process induced by gravitational cooling in the high and low momentum regimes and estimate the relaxation time scales in each case. We detect a low frequency mode in the whole parameters space and it was found that the resulting configuration oscillates under this mode with amplitude that depends on the mass ratio and head-on momentum. In some cases the resulting configuration oscillates with changes in density of two orders of magnitude and with a matter distribution that is far from isotropic. These results could contribute to the collection of possible mass distributions considered in the reconstruction of mass profiles obtained in structure formation simulations.

We present a 3D special-relativistic radiation hydrodynamics code. It uses the radiative inversion scheme with the M1-closure relation for the radiation equations, which allows the treatment of a wide range of optical depth, temperature and opacity. The radiation field is treated in the grey-body approximation. We present the standard 1D and 2D tests that include both optically thin and thick scenarios, as well as hydrodynamical and radiation pressure dominated configurations. As an application in 3D, we show the evolution of a jet driven by radiation-hydrodynamics with a helical perturbation. The code is expected to allow the exploration of scenarios in high-energy astrophysics where the radiation is important, like sources of GRBs.

Within the framework of resistive MHD, implementing the C7 equilibrium atmosphere model and a 3D potential magnetic field realistic configuration, we simulate the formation of a plasm jet with the morphology, upward velocity up to 130 km/s and timescale formation between 60 and 90 s after beginning of simulation, similar to those expected for Type II spicules. Initial results of this simulation were published in Paper (e.g., Gonz\'alez-Avil\'es et al. 2018) and present paper is devoted to the analysis of transverse displacements and rotational type motion of the jet. Our results suggest that 3D magnetic reconnection may be responsible for the formation of the jet in Paper (Gonz\'alez-Avil\'es et al. 2018). In this paper, by calculating times series of the velocity components $v_x$ and $v_y$ in different points near to the jet for various heights we find transverse oscillations in agreement with spicule observations. We also obtain a time-distance plot of the temperature in a cross-cut at the plane $x=$0.1 Mm and find significant transverse displacements of the jet. By analyzing temperature isosurfaces of $10^{4}$ K with the distribution of $v_x$, we find that if the line-of-sight (LOS) is approximately perpendicular to the jet axis then there is both motion towards and away from the observer across the width of the jet. This red-blue shift pattern of the jet is caused by rotational motion, initially clockwise and anti-clockwise afterwards, which could be interpreted as torsional motion. From a nearly vertical perspective of the jet the LOS velocity component shows a central blue-shift region surrounded by red-shifted plasma.

We study the process of magnetic field annihilation and reconnection in simulations of magnetised solar photosphere and chromosphere with magnetic fields of opposite polarities and constant numerical resistivity. Exact analytical solutions for reconnective annihilations are used to interpret the features of magnetic reconnection in simulations of flux cancellation in the solar atmosphere. We use MURaM high-resolution photospheric radiative magneto-convection simulations to demonstrate the presence of magnetic field reconnection consistent with the magnetic flux pile-up models. Also, a simulated data-driven chromospheric magneto-hydrodynamic simulation is used to demonstrate magnetic field and flow structures, which are similar to the ones theoretically predicted. Both simulations demonstrate flow and magnetic field structures roughly consistent with accelerated reconnection with magnetic flux pile-up. The presence of standard Sweet-Parker type reconnection is also demonstrated in stronger photospheric magnetic fields.

Modeling emission lines from the millimeter to the UV and producing synthetic spectra is crucial for a good understanding of observations, yet it is an art filled with hazards. This is the proceedings of "Walking the Line", a 3-day conference held in 2018 that brought together scientists working on different aspects of emission line simulations, in order to share knowledge and discuss the methodology. Emission lines across the spectrum from the millimeter to the UV were discussed, with most of the focus on the interstellar medium, but also some topics on the circumgalactic medium. The most important quality of a useful model is a good synergy with observations and experiments. Challenges in simulating line emission are identified, some of which are already being worked upon, and others that must be addressed in the future for models to agree with observations. Recent advances in several areas aiming at achieving that synergy are summarized here, from micro-physical to galactic and circum-galactic scale.

We present the three-dimensional numerical simulation of jets modeled with Relativistic Radiation Hydrodynamics (RRH), that evolve across two environments: i) a stratified surrounding medium and ii) a 16TI progenitor model. We consider opacities consistent with various processes of interaction between the fluid and radiation, specifically, free-free, bound-free, bound-bound and electron scattering. We explore various initial conditions, with different radiation energy densities of the beam in hydrodynamical and radiation pressure dominated scenarios, considering only ultra-relativistic jets. In order to investigate the impact of the radiation field on the evolution of the jets, we compare our results with purely hydrodynamical jets. Comparing among jets driven by RRH, we find that radiation pressure dominated jets propagate slightly faster than gas pressure dominated ones. Finally, we construct the luminosity Light Curves (LCs) associated with all these cases. The construction of LCs uses the fluxes of the radiation field which is fully coupled to the hydrodynamics equations during the evolution. The main properties of the jets propagating on the stratified surrounding medium are that the LCs show the same order of magnitude as the gamma-ray luminosity of typical Long Gamma-Ray Bursts $10^{50}-10^{54}$erg/s and the difference between the radiation and gas temperatures is of nearly one order of magnitude. The properties of jets breaking out from the progenitor star model are that the LCs are of the order of magnitude of low-luminosity GRBs $10^{46}-10^{49}$ erg/s, and in this scenario the difference between the gas and radiation temperature is of four orders of magnitude, which is a case far from thermal equilibrium.

Using simulated data-driven three-dimensional resistive MHD simulations of the solar atmosphere, we show that magnetic reconnection can be responsible of the formation of jets with characteristic of Type II spicules. For this, we numerically model the photosphere-corona region using the C7 equilibrium atmosphere model. The initial magnetic configuration is a 3D potential magnetic field, extrapolated up to the solar corona region from a dynamic realistic simulation of solar photospheric magnetoconvection model which is mimicking quiet-Sun. In this case we consider a uniform and constant value of the magnetic resistivity of $12.56 ~\Omega~{\rm m}$. We have found that formation of the jets depends on the Lorentz force, which helps to accelerate the plasma upwards. Analyzing various properties of the jet dynamics, we found that the jet structure shows Doppler shift near to regions with high vorticity. The morphology, upward velocity, covering a range up to 100 $\rm km$ $\rm s^{-1}$, and life-time of the estructure, bigger than 100 ${\rm s}$, are similar to those expected for Type II spicules.

We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.

We present a method to classify initial conditions of a long gamma ray bursts model sourced by a single relativistic shock. It is based on the use of artificial neural networks (ANNs) that are trained with light curves (LC) generated with radiation relativistic hydrodynamics simulations. The model we use consists in a single shock with a highly relativistic injected beam into a stratified surrounding medium with profile $1/r^2$. In the process we only consider the bremsstrahlung radiation and Thomson scattering process. The initial conditions we use to train the ANN are three: the rest mass density, Lorentz factor and radiation energy density of the beam that produces the relativistic shock, together with the LC generated during the process. The classification selects the location of a box in the 3d parameter space that better fits a given LC, and in order to decrease the uncertainty of the parameters this box is refined and the classification selects a new box of smaller size.

We compare the results of the semi-classical (SC) and quantum-mechanical (QM) formalisms for angular-momentum changing transitions in Rydberg atom collisions given by Vrinceanu & Flannery, J. Phys. B 34, L1 (2001), and Vrinceanu, Onofrio & Sadeghpour, ApJ 747, 56 (2012), with those of the SC formalism using a modified Monte Carlo realization. We find that this revised SC formalism agrees well with the QM results. This provides further evidence that the rates derived from the QM treatment are appropriate to be used when modelling recombination through Rydberg cascades, an important process in understanding the state of material in the early universe. The rates for $\Delta\ell=\pm1$ derived from the QM formalism diverge when integrated to sufficiently large impact parameter, $b$. Further to the empirical limits to the $b$ integration suggested by Pengelly & Seaton, MNRAS 127, 165 (1964), we suggest that the fundamental issue causing this divergence in the theory is that it does not fully cater for the finite time taken for such distant collisions to complete.

Accurate rates for energy-degenerate l-changing collisions are needed to determine cosmological abundances and recombination. There are now several competing theories for the treatment of this process, and it is not possible to test these experimentally. We show that the H I two-photon continuum produced by astrophysical nebulae is strongly affected by l-changing collisions. We perform an analysis of the different underlying atomic processes and simulate the recombination and two-photon spectrum of a nebula containing H and He. We provide an extended set of effective recombination coefficients and updated l-changing 2s-2p transition rates using several competing theories. In principle, accurate astronomical observations could determine which theory is correct.

Using numerical simulations, we show that jets with features of type II spicules and cold coronal jets corresponding to temperatures $10^{4}$ K can be formed due to magnetic reconnection in a scenario in presence of magnetic resistivity. For this we model the low chromosphere-corona region using the C7 equilibrium solar atmosphere model and assuming Resistive MHD rules the dynamics of the plasma. The magnetic filed configurations we analyze correspond to two neighboring loops with opposite polarity. The separation of the loops' feet determines the thickness of a current sheet that triggers a magnetic reconnection process, and the further formation of a high speed and sharp structure. We analyze the cases where the magnetic filed strength of the two loops is equal and different. In the first case, with a symmetric configuration the spicules raise vertically whereas in an asymmetric configuration the structure shows an inclination. With a number of simulations carried out under a 2.5D approach, we explore various properties of excited jets, namely, the morphology, inclination and velocity. The parameter space involves magnetic field strength between 20 and 40 G, and the resistivity is assumed to be uniform with a constant value of the order $10^{-2}\Omega\cdot m$

Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.

We present an algorithm based on artificial neural networks (ANNs), that estimates the mass ratio in a binary black hole collision out of given Gravitational Wave (GW) strains. In this analysis, the ANN is trained with a sample of GW signals generated with numerical simulations. The effectiveness of the algorithm is evaluated with GWs generated also with simulations for given mass ratios unknown to the ANN. We measure the accuracy of the algorithm in the interpolation and extrapolation regimes. We present the results for noise free signals and signals contaminated with Gaussian noise, in order to foresee the dependence of the method accuracy in terms of the signal to noise ratio.

Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction.

During the On-Station Thermal Test campaign of the LISA Pathfinder the data and diagnostics subsystem was tested in nearly space conditions for the first time after integration in the satellite. The results showed the compliance of the temperature measurement system, obtaining temperature noise around $10^{-4}\,{\rm K}\, {\rm Hz}^{-1/2}$ in the frequency band of $1-30\;{\rm mHz}$. In addition, controlled injection of heat signals to the suspension struts anchoring the LISA Technology Package (LTP) Core Assembly to the satellite structure allowed to experimentally estimate for the first time the phase noise contribution through thermo-elastic distortion of the LTP interferometer, the satellite's main instrument. Such contribution was found to be at $10^{-12}\,{\rm m}\, {\rm Hz}^{-1/2}$, a factor of 30 below the measured noise at the lower end of the measurement bandwidth ($1\,{\rm mHz}$).

Deep phase modulation interferometry was proposed as a method to enhance homodyne interferometers to work over many fringes. In this scheme, a sinusoidal phase modulation is applied in one arm while the demodulation takes place as a post-processing step. In this contribution we report on the development to implement this scheme in a fiber coupled interferometer controlled by means of a FPGA, which includes a LEON3 soft-core processor. The latter acts as a CPU and executes a custom made application to communicate with a host PC. In contrast to usual FPGA-based designs, this implementation allows a real-time fine tuning of the parameters involved in the setup, from the control to the post-processing parameters.

LISA Pathfinder is a mission planned by the European Space Agency to test the key technologies that will allow the detection of gravitational waves in space. The instrument on-board, the LISA Technology package, will undergo an exhaustive campaign of calibrations and noise characterisation campaigns in order to fully describe the noise model. Data analysis plays an important role in the mission and for that reason the data analysis team has been developing a toolbox which contains all the functionalities required during operations. In this contribution we give an overview of recent activities, focusing on the improvements in the modelling of the instrument and in the data analysis campaigns performed both with real and simulated data.

This paper presents a quantitative assessment of the performance of the upcoming LISA Pathfinder geodesic explorer mission. The findings are based on the results of extensive ground testing and simulation campaigns using flight hardware and flight control and operations algorithms. The results show that, for the central experiment of measuring the stray differential acceleration between the LISA test masses, LISA Pathfinder will be able to verify the overall acceleration noise to within a factor two of the LISA requirement at 1 mHz and within a factor 6 at 0.1 mHz. We also discuss the key elements of the physical model of disturbances, coming from LISA Pathfinder and ground measurement, that will guarantee the LISA performance.

We numerically solve the mass-less test scalar field equation on the space-time background of boson stars and black holes. In order to do so, we use a numerical domain that contains future null infinity. We achieve this construction using a scri-fixing conformal compactification technique based on hyperboloidal constant mean curvature foliations of the space-time and solve the conformally invariant wave equation. We present two results: the scalar field shows oscillations of the quasi- normal-mode type found for black holes only for boson star configurations that are compact, and no signs of tail decay is found in the parameter space we explored. Even though our results do not correspond to the master equation of perturbations of boson star solutions, they indicate that the parameter space of boson stars as black hole mimickers is restricted to compact configurations.

We present and constrain a cosmological model which component is a pressureless fluid with bulk viscosity as an explanation for the present accelerated expansion of the universe. We study the particular model of a constant bulk viscosity coefficient \zeta_m. The possible values of \zeta_m are constrained using the cosmological tests of SNe Ia Gold 2006 sample, the CMB shift parameter R from the three-year WMAP observations, the Baryon Acoustic Oscillation (BAO) peak A from the Sloan Digital Sky Survey (SDSS) and the Second Law of Thermodynamics (SLT). It was found that this model is in agreement with the SLT using only the SNe Ia test. However when the model is submitted to the three cosmological tests together (SNe+CMB+BAO) the results are: 1.- the model violates the SLT, 2.- predicts a value of H_0 \approx 53 km sec^{-1} Mpc^{-1} for the Hubble constant, and 3.- we obtain a bad fit to data with a \chi^2_{min} \approx 400 (\chi^2_{d.o.f.} \approx 2.2). These results indicate that this model is ruled out by the observations.