Detection of Earth-mass planets with the radial velocity method requires a precision of about 10cm/s to identify a signal caused by such a planet. At the same time, noise originating in the photospheric and subphotospheric layers of the parent star is of the order of meters per second. Understanding the physical nature of the photospheric noise (so-called stellar jitter) and characterizing it are critical for developing techniques to filter out these unwanted signals. We take advantage of current computational and technological capabilities to create 3D realistic models of stellar subsurface convection and atmospheres to characterize the photospheric jitter. We present 3D radiative hydrodynamic models of several solar-type target stars of various masses and metallicities, discuss how the turbulent surface dynamics and spectral line characteristics depend on stellar properties, and provide stellar jitter estimates for these stars.

Photoelectric effect of dust grains by UV radiation is an important process for disk heating, but as a disk evolves, the amount of dust grains decreases. Photoeaporation is a disk dispersal process, which is caused by high-energy radiation. We perform a set of photoevaporation simulations solving hydrodynamics with radiative transfer and non-equilibrium chemistry in a self-consistent way. We run a series of simulations with varying the dust-to-gas mass ratio in a range . We show that H2 pumping and X-ray heating mainly contribute to the disk heating in case of and photoelectric effect mainly heats the gas in cases. The mass-loss profile changes significantly with respect to the main heating process. The outer disk is more efficiently dispersed when photoelectric effect is the main heating source.

During the recombination of the universe, supersonic relative motion between baryons and dark matter (DM) generally existed. In the presence of such streaming motions, gas clumps can collapse outside of virial radii of their closest dark matter halos. Such baryon dominant objects are thought to be self-gravitating and are called supersonically induced gas objects; SIGOs. We perform three-dimensional hydrodynamical simulations by including H2 chemical reactions and stream velocity and follow SIGO’s formation from z = 200 to z = 25. SIGOs can be formed under the influence of stream velocity, and cooling is effective in contracting gas clouds. We follow its further evolution with higher resolution. We find that there are SIGOs which become Jeans unstable outside of the virial radius of the closest DM halos. Those SIGOs are gravitationally unstable and trigger star formation.

We introduce hydrodynamic simulations in which a protostar captures a cloudlet with a relatively small angular momentum. The cloudlet accretes onto the protostar and perturbs the gas disk rotating around the protostar. This cloudlet capture can reproduce some features observed in the molecular emission lines from TMC-1A. First, the cloudlet can reproduce the blue asymmetry observed in the CS emission. Second, the cloudlet can explain the slow infall observed in the C18O emission. Third, the impact of the cloudlet can explain the offset of the SO emission from the disk center. We also argue that a warm gas should confine the cloudlet through pressure. A protostar may obtain substantial mass by capturing cloudlets.

We present an overview of PION, an open-source software project for solving radiation-magnetohydrodynamics equations on a nested grid, aimed at modelling asymmetric nebulae around massive stars. A new implementation of hybrid OpenMP/MPI parallel algorithms is briefly introduced, and improved scaling is demonstrated compared with the current release version. Three-dimensional simulations of an expanding nebula around a Wolf-Rayet star are then presented and analysed, similar to previous 2D simulations in the literature. The evolution of the emission measure of the gas and the X-ray surface brightness are calculated as a function of time, and some qualitative comparison with observations is made.

The Bisous model is a tool that uses stochastic methods to detect the network of galactic filaments. This model is explicitly developed to detect the structure from observational data, using only galaxy positions as input. This paper shows that the Bisous model gives reliable results and including photometric data improves the resulting filamentary network. We used MultiDark-Galaxies catalogue to create a mock with photometric redshifts and samples with different galaxy number densities. We found that the filaments detected with the Bisous model are reliable; 85% of the detected filaments are unchanged compared to results with more complete input data. Adding photometric data improves the fraction of galaxies in filaments. Using the confusion matrix technique, we found the false discovery rate to always be below 5% when using photometric data.

We analyzed 39 ks NuSTAR data of Cen X-3 through both orbital- and pulse-phase resolved spectroscopy. Orbital-phase resolved spectra show extrinsic fluctuations due to absorption by surrounding plasma, as the spectral fluctuation mainly emerges below 10 keV. Pulse-phase resolved spectra, on the other hand, show intrinsic fluctuations depending on effectiveness of Comptonization, since the spectrum becomes hard above 10 keV at the pulse peak.

Gaia data allows for search for extended stellar structures in phase (coordinates plus velocities) space. We describe a method of using DBSCAN clustering algorithm, which is used to group closely-packed-together data points, to a list of preliminary selected pairs of stars, with parameters expected to be found within stellar streams and comoving groups: loose structures in which stars are not gravitationally bound, but do share motion and evolutionary properties. To test our approach, we construct a model population of background stars, and use pair-constructing and clustering algorithms on it. Results show that transitioning to a list of pairs sharply reveals structures not presented in background model, which then become more apparent targets in coordinate-velocity phase space for DBSCAN algorithm thanks to now increased relative density of the extended stellar structure.

The INTEGRAL satellite, in orbit since October 2002, has significantly contributed to the study of magnetars and, thanks to its unique capabilities for the study of transient gamma-ray phenomena, it is now playing an important role in multimessenger astrophysics. The most recent results include the discovery of a peculiar burst from SGR J1935+2154, which gave the first observational evidence for the connection between magnetars and fast radio bursts, and extensive searches for bursting activity in peculiar sources, such as the repeating FRB 20200120E in M81 and the ultra-long period magnetar candidate GLEAMX J162759.5–523504.3.

In this paper I present a method to enhance the search sensitivity for long transient Gravitational Waves produced by Neutron Star instabilities. This method consists in a selective image filter, called Triangular Filter, to be applied to data spectrograms. It is shown that thanks to this implementation a ∼20% gain in sensitivity is achievable.

In this paper we have inferred the magnetic shielding characteristics and space weather hazards of selected potentially habitable extrasolar planets using a dynamical geophysical model from calculations of internal heat, phases of volcanism and planetary magnetic moments. The space weather hazards on the extrasolar planet Kepler-452b orbiting around a Sun-like star are found to be a minimum which enhances the habitability probability of this planet.

We present Ekster, a new method for simulating the formation and dynamics of individual stars in a relatively low-resolution gas background. Here, we use Ekster to simulate star cluster formation in two different regions from each of two galaxy models with different spiral potentials. We simulate these regions for 3 Myr to study where and how star clusters form. We find that massive GMC regions form more massive clusters than sections of spiral arms. Additionally we find that clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.

The Laser Interferometer Space Antenna (LISA) mission will observe from space gravitational waves emitted by neutron stars and white dwarfs within galactic binaries. These compact stars can have intense magnetic fields. Therefore, the impact of the magnetic fields on the orbital and the spins evolution of binary systems can potentially be detected by LISA through the GW’s strain. Within the magnetic dipole-dipole approximation, we found that magnetism generates a secular drift of the mean longitude which, in turn, shifts all the frequencies contained in the GW signal. For a quasi-circular orbit, the signal is mainly monochromatic and the magnetic shift is proportional to the product of the magnetic moments and is inversely proportional to the 7/2 power of the semi-major axis. Hence, for a highly magnetic binary system in compact orbit, a non-negligible amount of the frequency measured by LISA might have a magnetic origin.

One of the models explaining the high luminosity of pulsing ultra-luminous X-ray sources (pULXs) was suggested by Mushtukov et al. (2015). They showed that the accretion columns on the surfaces of highly magnetized neutron stars can be very luminous due to opacity reduction in the high magnetic field. However, a strong magnetic field leads also to amplification of the electron-positron pairs creation. Therefore, increasing of the electron and positron number densities compensates the cross-section reduction, and the electron scattering opacity does not decrease with the magnetic field magnification. As a result, the maximum possible luminosity of the accretion column does not increase with the magnetic field. It ranges between 1040 − 1041 erg s−1 depending only slightly on the magnetic field strength.

In recent years, an increasing amount of attention is being paid to the gravitational few-body problem and its applications to astrophysical scenarios. Among the main reasons for this renewed interest there is large number of newly discovered exoplanets and the detection of gravitational waves. Here, we present two numerical codes to model three- and few-body systems, called tsunami and okinami. The tsunami code is a direct few-body code with algorithmic regularization, tidal forces and post-Newtonian corrections. okinami is a secular, double-averaged code for stable hierarchical triples. We describe the main methods implemented in our codes, and review our recent results and applications to gravitational-wave astronomy, planetary science and statistical escape theories.

Space interferometry is the inevitable end point of high angular resolution astrophysics, and a key technology that can be leveraged to analyse exoplanet formation and atmospheres with exceptional detail. However, the anticipated cost of large missions, such as Darwin and TPF-I, and inadequate technology readiness levels have resulted in limited developments since the late 2000s. Here, we present a feasibility study into a small-scale formation-flying interferometric array in low Earth orbit, which will aim to prove the technical concepts involved with space interferometry while still making unique astrophysical measurements. We will detail the proposed system architecture and metrology system, as well as present orbital simulations that show that the array should be stable enough to perform interferometry with <50 m s–1 yr–1 delta-v and one thruster per spacecraft. We also conduct observability simulations to identify which parts of the sky are visible for a given orbital configuration. We conclude with optimism that this design is achievable, but a more detailed control simulation factoring in a demonstrated metrology system is the next step to demonstrate full mission feasibility.

We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxy’s stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $\lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpin’s current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy.