Search

The (re)analysis of data on the X-ray emitting pulsars PSR J0030+0451 and PSR J0740+6620, as well as new results on PSR J0437-4715, are confronted with the predictions of the equation of state (EoS) models allowing for strong first-order phase transition for the mass-radius ($M$-$R$) diagram. We use models that are based on a covariant density functional (CDF) EoS for nucleonic matter at low densities and a quark matter EoS, parameterized by the speed of sound, at higher densities. To account for the variations in the ellipses for PSR J0030+0451 obtained from different analyses, we examined three scenarios to assess their consistency with our models, focusing particularly on the potential formation of twin stars. We found that in two scenarios, where the ellipses for PSR J0030+0451 and PSR J0437-4715 with masses close to the canonical mass $\sim 1.4\,M_{\odot}$ are significantly separated, our models allow for the presence of twin stars as a natural explanation for potential differences in the radii of these stars., Comment: 10 pages, 4 figures

In this study, we examine thermal conductivity and the thermal Hall effect in electron-ion plasmas relevant to hot neutron stars, white dwarfs, and binary neutron star mergers, focusing on densities found in the outer crusts of neutron stars and the interiors of white dwarfs. We consider plasma consisting of single species of ions, which could be either iron $\isotope[56]{Fe}$ or carbon $\isotope[12]{C}$ nuclei. The temperature range explored is from the melting temperature of the solid $T\sim10^9$~K up to $10^{11}$~K. This covers both degenerate and non-degenerate electron regimes. We find that thermal conductivity increases with density and temperature for which we provide analytical scaling relations valid in different regimes. The impact of magnetic fields on thermal conductivity is also analyzed, showing anisotropy in low-density regions and the presence of the thermal Hall effect characterized by the Righi--Leduc coefficient. The transition from a degenerate to non-degenerate regime is characterized by a minimum ratio of thermal conductivity to temperature, which is analogous to the minimum observed already in the case of electrical conductivity. We provide also formulas fit to our numerical results, which can be used in dissipative magneto-hydrodynamics simulations of warm compact stars., Comment: 17 pages, 10 figures

We demonstrate that the interface between $S$-wave and $P$-wave paired superfluids in neutron stars induces a neutron supercurrent, akin to the Josephson junction effect in electronic superconductors. The proton supercurrent entrainment by the neutron superfluid generates, in addition to the neutral supercurrent, a charged current across the interface. Beyond the stationary limit, the motions of the neutron vortex line and proton flux tube arrays, responding to secular changes in the neutron star's rotation rate, induce a time-dependent oscillating Josephson current across this interface. We show that such motion produces radiation from the interface, which is phenomenologically significant enough to heat the star and alter its cooling rate during the photon cooling era., Comment: 6 pages, 3 figures. Comments welcome

We study the bulk viscosity of moderately hot and dense, neutrino-transparent color superconducting quark matter arising from weak-interaction-driven direct URCA processes. The quark matter is modeled using the Nambu--Jona-Lasinio model improved to account for vector and 't Hooft interactions as well as antisymmetric pairing among the red/green up and down quarks. The unpaired excitations are the strange quarks and the blue up or down quarks. We compute the relaxation rates associated with $d$ and $s$-quark decay and electron capture processes on $u$ quark for blue color. The resulting bulk viscosity for density oscillations in the 1--10\,kHz range shows a resonant peak at $T\sim 5\,\MeV$, and the damping time may drop below 10 ms. This is short enough to affect the postmerger evolution and is very similar to the damping predicted in nuclear matter., Comment: v2: 6 pages + 4 figures main text, 2 pages + 2 figures supplemental material. Matches version published in Phys. Rev. D. Letters

Recent measurements of high-momentum correlated neutron-proton pairs at JLab suggest that the dense nucleonic component of the compact stars contains a fraction of high-momentum neutron-proton pairs that is not accounted for in the familiar Fermi-liquid theory of the neutron-proton fluid mixture. We compute the rate of the Urca process in compact stars taking into account the non-Fermi liquid contributions to the proton's spectral widths induced by short-range correlations. The Urca rate differs strongly from the Fermi-liquid prediction at low temperatures, in particular, the high threshold on the proton fraction precluding the Urca process in neutron stars is replaced by a smooth increase with the proton fraction. This observation may have a profound impact on the theories of cooling of compact stars., Comment: v3: added references and fixed typos. Matches the published version, 6 pages, 2 figures. v2: Typos fixed; v1: 5 pages, 2 figures

We construct tables of finite temperature equation of state (EoS) of hypernuclear matter in the range of densities, temperatures, and electron fractions that are needed for numerical simulations of supernovas, proto-neutron stars, and binary neutron star mergers and cast them in the format of {\sc CompOSE} database. The tables are extracted from a model that is based on covariant density functional (CDF) theory that includes the full $J^P=1/2^+$ baryon octet in a manner that is consistent with the current astrophysical and nuclear constraints. We employ a parameterization with three different values of the slope of the symmetry energy $L_{\rm sym}=30$, 50 and 70 MeV and fixed skewness $Q_{\rm sat}= 400$ MeV for above saturation matter. A model for the EoS of inhomogeneous matter is matched at sub-saturation density to the high-density hypernuclear EoS. We discuss the generic features of the resulting EoS and the composition of matter as a function of density, temperature, and electron fraction. The nuclear characteristics and strangeness fraction of these models are compared to the alternatives from the literature. The integral properties of static and rapidly rotating compact stars in the limit of zero temperature are discussed and confronted with the multimessenger astrophysical constraints., Comment: 18 pages, 9 figures; v2 matches published version

Recent astrophysical mass inferences of compact stars HESS J1731-347 and PSR J0952-0607, with extremely small and large masses respectively, as well as the measurement of the neutron skin of Ca in the CREX experiment challenge and constrain the models of dense matter. We examine the concept of hybrid stars - objects containing quark cores surrounded by nucleonic envelopes - as models that account for these new data along with other inferences. We employ a family of 81 nucleonic equations of state (EoSs) with variable skewness and slope of symmetry energy at saturation density and a constant speed-of-sound EoS for quark matter. For each nucleonic EoS, a family of hybrid EoSs is generated by varying the transition density, the energy jump, and the speed of sound. These models are tested against the data from GW170817 and J1731-347, which favor low-density soft EoS and J0592-0607 and J0740+6620, which require high-density stiff EoS. The addition of J0592-0607's mass measurement to the constraints has no significant impact on the parameter space of the admissible EoS, but allows us to explore the potential effect of pulsars more massive than J0740+6620, if such exists. We then examine the occurrence of twin configurations and quantify the ranges of masses and radii that they can possess. It is shown that including J1731-347 data favors EoSs that predict low-mass twins with $M \lesssim 1.3\,M_{\odot}$ that can be realized if the deconfinement transition density is low. If combined with large speed of sound in quark matter such models allow for maximum masses of hybrid stars in $2.0$--$2.6\,M_{\odot}$., Comment: v2: 8 pages, 6 figures, 1 table, expanded data and discussion, matches published version

We show analytically that there exist compact stellar objects akin to neutron stars whose radius is smaller than the Schwarzschild radius defined by Arnowitt-Deser-Misner (ADM) mass. The radius of the compact object is defined by the radius where the energy density and the pressure of ordinary matter vanish, while clouds of scalar(s) can extend beyond this radius -- a situation that is often encountered in modified gravity theories, like $F(R)$ gravity and the scalar--Einstein--Gauss-Bonnet gravity. The clouds of scalar mode(s) give additional contributions to the ADM mass and as a result, the corresponding Schwarzschild radius given by the ADM mass can be larger than that of the compact object., Comment: LaTeX, 13 pages, no figure, version to appear in Nuclear Physics B

This article presents the development of a new wind turbine simulation software to study wake flow physics. To this end, the design and development of waLBerla-wind, a new simulator based on the lattice-Boltzmann method that is known for its excellent performance and scaling properties, will be presented. Here it will be used for large eddy simulations (LES) coupled with actuator wind turbine models. Due to its modular software design, waLBerla-wind is flexible and extensible with regard to turbine configurations. Additionally it is performance portable across different hardware architectures, another critical design goal. The new solver is validated by presenting force distributions and velocity profiles and comparing them with experimental data and a vortex solver. Furthermore, waLBerla-wind's performance is \revision{compared to a theoretical peak performance}, and analysed with weak and strong scaling benchmarks on CPU and GPU systems. This analysis demonstrates the suitability for large-scale applications and future cost-effective full wind farm simulations.

Neutron stars may experience differential rotation on short, dynamical timescales following extreme astrophysical events like binary neutron star mergers. In this work, the masses and radii of differentially rotating neutron star models are computed. We employ a set of equations of states for dense hypernuclear and $\Delta$-admixed-hypernuclear matter obtained within the framework of CDF theory in the relativistic Hartree-Fock (RHF) approximation. Results are shown for varying meson-$\Delta$ couplings, or equivalently the $\Delta$-potential in nuclear matter. A comparison of our results with those obtained for non-rotating stars shows that the maximum mass difference between differentially rotating and static stars is independent of the underlying particle composition of the star. We further find that the decrease in the radii and increase in the maximum masses of stellar models when $\Delta$-isobars are added to hyperonuclear matter (as initially observed for static and uniformly rotating stars) persist also in the case of differentially rotating neutron stars., Comment: 7 pages, 6 figures

We compute the anisotropic electrical conductivity tensor of the inner crust of a compact star at non-zero temperature by extending a previous work on the conductivity of the outer crust. The physical scenarios, where such crust is formed, involve proto-neutron stars born in supernova explosions, binary neutron star mergers and accreting neutron stars. The temperature-density range studied covers the transition from a non-degenerate to a highly degenerate electron gas and assumes that the nuclei form a liquid, i.e., the temperature is above the melting temperature of the lattice of nuclei. The electronic transition probabilities include (a) the dynamical screening of electron-ion interaction in the hard-thermal-loop approximation for the QED plasma, (b) the correlations of the ionic component in a one-component plasma, and (c) finite nuclear size effects. The conductivity tensor is obtained from the Boltzmann kinetic equation in relaxation time approximation accounting for the anisotropies introduced by a magnetic field. The sensitivity of the results towards the matter composition of the inner crust is explored by using several compositions of the inner crust which were obtained using different nuclear interactions and methods of solving the many-body problem. The standard deviation of relaxation time and components of the conductivity tensor from the average are below $\le 10\%$ except close to crust-core transition, where non-spherical nuclear structures are expected. Our results can be used in dissipative magneto-hydrodynamics (MHD) simulations of warm compact stars., Comment: 13 pages, 13 figures, 5 tables

We generate three families of extended covariant density functionals of nuclear matter that have varying slope of symmetry energy and skewness at nuclear saturation density, but otherwise share the same basic parameters (symmetry energy, compressibility, saturation parameters, etc.) with the standard DDME2, DD2, and MPE functionals. Tables of the parameters of these new density functionals are given, which can be straightforwardly used in DDME2, DD2, and MPE parameterization-based codes. Furthermore, we provide tables of a large number of equations of state (81 for each family) that can be used in astrophysical simulations to assess the impact of variations of not-well-known slope of symmetry energy and skewness of nuclear systems on the astrophysics of compact objects. We also provide tables of computed integral parameters (mass, radius, and tidal deformability) that can be used, e.g., for modeling gravitational waveforms. Finally, for the extended DDME2-based parameterization, we implement a first-order phase transition to quark matter to obtain a family of equations of state that accommodates a phase transition to quark matter. Analogous tables of the equations of state and integral parameters are provided for this case as well., Comment: v2: matches published version. v1: 8 pages, 4 figures, 5 tables

A set of hadronic equations of state derived from covariant density functional theory and constrained by terrestrial experiments, and astrophysical observations, in particular by the NICER experiment inferences is used to explore the universal relations among the global properties of compact stars containing heavy baryons at high densities. We confirm the validity of universal $I$-Love-$Q$ relations connecting the moment of inertia $(I)$, the tidal deformability ($\Lambda$), and the spin-induced quadrupole moment ($Q$) for isolated non-rotating stars. We further confirm the validity of the $I$-$C$-$Q$ relations connecting the moment of inertia, compactness $(C)$, and quadrupole moment for uniformly and slowly rotating stars, and extend the validity of these relations to maximally rotating sequences. We then investigate the relations between integral parameters of maximally rotating and static compact stars. The universalities are shown to persist for equations of state and compositions containing hyperons and $\Delta$ degrees of freedom. When heavy baryons are included, however, the radial profiles of integrands in expressions of global properties exhibit ``bumps", which are not present in the case of nucleonic stars in which case the profiles are smooth. We determine the coefficients entering the universal relations in the case of hyperonic and $\Delta$-resonance containing stars., Comment: 19 pages, 15 figures; matches published version

The recent attempt on mass and radius inference of the central compact object within the supernova remnant HESS J1731-347 suggests for this object an unusually low mass of $M = 0.77^{+0.20}_{-0.17}\,M_{\odot}$ and a small radius of $R = 10.4^{+0.86}_{-0.78}$\,km. We explore the ways such a result can be accommodated within models of dense matter with heavy baryonic degrees of freedom which are constrained by the multi-messenger observations. We find that to do so using only purely nucleonic models, one needs to assume a rather small value of the slope of symmetry energy $L_{\rm sym}$. Once heavy baryons are included higher values of the slope $L_{\rm sym}$ become acceptable at the cost of a slightly reduced maximum mass of static configuration. These two scenarios are distinguished by the particle composition and will undergo different cooling scenarios. In addition, we show that the universalities of the $I$-Love-$Q$ relations for static configurations can be extended to very low masses without loss in their accuracy., Comment: v2: matches published version, 8 pages, 6 figures

This review covers several recent developments in the physics of dense QCD with an emphasis on the impact of multiple phase transitions on astrophysical manifestations of compact stars. It is conjectured that pair-correlated quark matter in $\beta$-equilibrium is within the same universality class as spin-imbalanced cold atoms and the isospin asymmetrical nucleonic matter. This then implies the emergence of phases with broken space symmetries and tri-critical (Lifshitz) points. We construct an equation of state (EoS) that extends the two-phase EoS of dense quark matter within the constant speed of sound parameterization by adding a conformal fluid with a speed of sound $c_{\rm conf.}=1/\sqrt{3}$ at densities $\ge 10~n_{\rm sat}$, where $n_{\rm sat}$ is the saturation density. With this input, we construct static, spherically symmetrical compact hybrid stars in the mass--radius diagram, recover such features as the twins and triplets, and show that the transition to conformal fluid leads to the spiraling-in of the tracks in this diagram. Stars on the spirals are classically unstable with respect to the radial oscillations but can be stabilized if the conversion timescale between quark and nucleonic phases at their interface is larger than the oscillation period. Finally, we review the impact of a transition from high-temperature gapped to low-temperature gapless two-flavor phase on the thermal evolution of hybrid stars., Comment: 18 pages, 7 figures. Invited contribution to a Special Issue of PARTICLES: "Strong Interactions in the Standard Model: Massless Bosons to Compact Stars". Eds. M. Ding, C. D. Roberts, S. M. Schmidt; v2: minor typos corrected, matches the published version

We study the bulk viscosity of moderately hot and dense, neutrino-transparent relativistic $npe\mu$ matter arising from weak-interaction direct Urca processes. This work parallels our recent study of the bulk viscosity of $npe\mu$ matter with a trapped neutrino component. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations -- DDME2 (which does not allow for the low-temperature direct-Urca process at any density) and NL3 (which allows for low-temperature direct-Urca process above a low-density threshold). We compute the equilibration rates of Urca processes of neutron decay and lepton capture, as well as the rate of the muon decay, and find that the muon decay process is subdominant to the Urca processes at temperatures $T\geq 3$MeV in the case of DDME2 model and $T\geq 1$MeV in the case of NL3 model. Thus, the Urca-process-driven bulk viscosity is computed with the assumption that pure leptonic reactions are frozen. As a result the electronic and muonic Urca channels contribute to the bulk viscosity independently and at certain densities the bulk viscosity of $npe\mu$ matter shows instead of the standard one-peak (resonant) form a "flattened" shape. In the final step, we estimate the damping timescales of density oscillations by the bulk viscosity. We find that, e.g., at a typical oscillation frequency $f=1$kHz, the damping of oscillations is most efficient at temperatures $3\leq T\leq 5$MeV and densities $n_B\leq 2n_0$ where they can affect the evolution of the post-merger object., Comment: v2: 28 pages, 18 figures, matches published version. v1: 27 pages, 19 figures

This paper focuses on investigating high-order harmonic generation (HHG) in graphene quantum dots (GQDs) under intense near-infrared laser fields. To model the GQD and its interaction with the laser field, we utilize a mean-field approach. Our analysis of the HHG power spectrum reveals fine structures and a noticeable enhancement in cutoff harmonics due to the long-range correlations. We also demonstrate the essential role of Coulomb interaction in determining of harmonics intensities and cutoff position. Unlike atomic HHG, where the cutoff energy is proportional to the pump wave intensity, in GQDs the cutoff energy scales with the square root of the field strength amplitude. A detailed time-frequency analysis of the entire range of HHG spectrum is presented using a wavelet transform. The analysis reveals intricate details of the spectral and temporal fine structures of HHG, offering insights into the various HHG mechanisms in GQDs., Comment: 10 pages, 15 figures

We consider coherent dynamics of graphene charged carriers exposed to an intense few-cycle linearly polarized laser pulse. The results, obtained by solving the generalized semiconductor Bloch equations numerically in the Hartree-Fock approximation, taking into account many-body Coulomb interaction, demonstrate strong dependence of the valley polarization on the carrier-envelope phase (CEP), which is interpolated by the simple sinusoidal law. Then we consider harmonic generation in multi-cycle laser field by graphene preliminary exposed to an intense few-cycle laser pulse. We show that the second harmonic's intensity is a robust observable quantity that provides a gauge of CEP for pulse durations up to two optical cycles, corresponding to 40 $\mathrm{fs}$ at the wavelength of 6.2 $\mathrm{\mu m}$., Comment: 9 pages, 10 figures

As a second-order method, the Natural Gradient Descent (NGD) has the ability to accelerate training of neural networks. However, due to the prohibitive computational and memory costs of computing and inverting the Fisher Information Matrix (FIM), efficient approximations are necessary to make NGD scalable to Deep Neural Networks (DNNs). Many such approximations have been attempted. The most sophisticated of these is KFAC, which approximates the FIM as a block-diagonal matrix, where each block corresponds to a layer of the neural network. By doing so, KFAC ignores the interactions between different layers. In this work, we investigate the interest of restoring some low-frequency interactions between the layers by means of two-level methods. Inspired from domain decomposition, several two-level corrections to KFAC using different coarse spaces are proposed and assessed. The obtained results show that incorporating the layer interactions in this fashion does not really improve the performance of KFAC. This suggests that it is safe to discard the off-diagonal blocks of the FIM, since the block-diagonal approach is sufficiently robust, accurate and economical in computation time., Comment: Under Review

Computing accurate rate constants for catalytic events occurring at the surface of a given material represents a challenging task with multiple potential applications in chemistry. To address this question, we propose an approach based on a combination of the rare event sampling method called Adaptive Multilevel Splitting (AMS) and ab initio molecular dynamics (AIMD). The AMS method requires a one dimensional reaction coordinate to index the progress of the transition. Identifying a good reaction coordinate is difficult, especially for high dimensional problems such a those encountered in catalysis. We probe various approaches to build reaction coordinates such as Support Vector Machine and path collective variables. The AMS is implemented so as to communicate with a DFT-plane wave code. A relevant case study in catalysis: the change of conformation and the dissociation of a water molecule chemisorbed on the (100) $\gamma$-alumina surface is used to evaluate our approach. The calculated rate constants and transition mechanisms are discussed and compared to those obtained by a conventional static approach based on the Eyring-Polanyi equation with harmonic approximation. It is revealed that the AMS method may provide rate constants which are smaller than the static approach by up to two orders of magnitude due to entropic effects involved in the chemisorbed water., Comment: 33 pages, 9 figures

In this work, we perform a phenomenological derivation of the first- and second-order relativistic hydrodynamics of dissipative fluids. To set the stage, we start with a review of the ideal relativistic hydrodynamics from energy-momentum and particle number conservation equations. We then go on to discuss the matching conditions to local thermodynamical equilibrium, symmetries of the energy-momentum tensor, decomposition of dissipative processes according to their Lorentz structure, and finally, the definition of the fluid velocity in the Landau and Eckart frames. With this preparatory work, we first formulate the first-order (Navier-Stokes) relativistic hydrodynamics from the entropy flow equation, keeping only the first-order gradients of thermodynamical forces. A generalized form of diffusion terms is found with a matrix of diffusion coefficients describing the relative diffusion between various flavors. The procedure of finding the dissipative terms is then extended to the second order to obtain the most general form of dissipative function for multiflavor systems up to the second order in dissipative fluxes. The dissipative function now includes in addition to the usual second-order transport coefficients of Israel-Stewart theory also second-order diffusion between different flavors. The relaxation-type equations of second-order hydrodynamics are found from the requirement of positivity of the dissipation function, which features the finite relaxation times of various dissipative processes that guarantee the causality and stability of the fluid dynamics. These equations contain a complete set of nonlinear terms in the thermodynamic gradients and dissipative fluxes arising from the entropy current, which are not present in the conventional Israel-Stewart theory., Comment: 20 pages

In this work, we consider the properties of compact stars in which quark matter has low- and high-density phases that are separated by a first-order phase transition. Thus, unlike the commonly considered case of a single phase transition from hadronic to quark matter, our models of hybrid stars contain sequential phase transitions from hadronic matter to low- and then to high-density quark matter phases. We extend our previous study of the parameter space of hybrid stars with a single phase transition to those with sequential phase transitions, taking into account the constraints on the mass and radius of neutron stars from the NICER experiment, the experimental inferences of the neutron skin thickness of the lead nucleus by the PREX-II experiment, and constraints on the tidal deformability from the gravitational-wave event GW170817. We determine the range of the masses for which both twin and triplet configurations, i.e., identical-mass stars with two and three different values of radii, arise., Comment: v2: matches published version, 11 pages, 11 figures

We study compact stars with hybrid equations of state consisting of a nuclear outer region and two nested quark phases, each separated from the lower density phase by a strong first-order phase transition. The stability of these models is determined by calculating their radial oscillation modes with different conversion rates between adjacent phases and hence junction conditions for the modes at the phase separation interface between them. In the case when the timescale of transition is faster than the period of oscillations, we recover the traditional stability criterion implying that $\partial M/\partial\rho_c>0$ on the stable branch(es), where $M$ is the mass and $\rho_c$ is the central density. In the opposite limit of slow conversion, we find stable stellar multiplets beyond triplets consisting of stars that are stable by the usual criterion plus slow-conversion (denoted by $s$) hybrid stars with $\partial M/\partial\rho_c<0$ that are stabilized due to an alternative junction condition on the fluid displacement field at the interface reflecting the slow rate of conversion. We also study the properties of the reaction mode, the radial mode that only exists for stars with rapid (abbreviated by $r$) phase transitions, in stars with either two rapid phase transitions or alternating rapid and slow phase transitions for the first time. The implications of alternative junction conditions are also examined, with these conditions generally being found to provide stability properties similar to those for a slow conversion rate., Comment: 11 pages, 14 figures

We review the physics of hyperons and $\Delta$-resonances in dense matter in compact stars. The covariant density functional approach to the equation of state and composition of dense nuclear matter in the mean-field Hartree and Hartree-Fock approximation is presented, with regimes covering cold $\beta$-equilibrated matter, hot and dense matter with and without neutrinos relevant for the description of supernovas and binary neutron star mergers, as well as dilute expanding nuclear matter in collision experiments. We discuss the static properties of compact stars with hyperons and $\Delta$-resonances in light of constraints placed in recent years by the multimessenger astrophysics of compact stars on the compact stars' masses, radii, and tidal deformabilities. The effects of kaon condensation and strong magnetic fields on the composition of hypernuclear stars are also discussed. The properties of rapidly rotating compact hypernuclear stars are discussed and confronted with the observations of 2.5-2.8 solar mass compact objects in gravitational wave events. We further discuss the cooling of hypernuclear stars, the neutrino emission reactions, hyperonic pairing, and the mass hierarchy in the cooling curves that arises due to the onset of hyperons. The effects of hyperons and $\Delta$-resonances on the equation of state of hot nuclear matter in the dense regime, relevant for the transient astrophysical event and in the dilute regime relevant to the collider physics is discussed. The review closes with a discussion of universal relations among the integral parameters of hot and cold hypernuclear stars and their implications for the analysis of binary neutron star merger events., Comment: v3 and v2: 79 pages, 26 figures, final version (minor changes and additions, typos corrected). v1: 76 pages, 26 figures. arXiv admin note: text overlap with arXiv:2105.14050

The multiphoton excitation and high harmonic generation (HHG) processes are considered using the microscopic quantum theory of nonlinear interaction of strong coherent electromagnetic (EM) radiation with rectangular graphene quantum dot (RGQD). The dynamic Hartree-Fock approximation is developed for the consideration of the quantum dot-laser field nonlinear interaction at the nonadiabatic multiphoton excitation regime. The many-body Coulomb interaction is described in the extended Hubbard approximation. By numerical results, we show the significance of the RGQD lateral size, shape, and EM wavefield orientation in RGQD of the zigzag edge compear to the armchair edge in the HHG process allowing for increasing the cutoff photon energy and the quantum yield of higher harmonics., Comment: 17 pages, 14 figures

We construct a set of hyperonic equations of state (EoS) by assuming SU(3) symmetry within the baryon octet and by using a covariant density functional (CDF) theory approach. The low-density regions of our EoS are constrained by terrestrial experiments, while the high-density regime is modeled by systematically varying the nuclear matter skewness coefficient $Q_{\rm sat}$ and the symmetry energy slope $L_{\rm sym}$. The sensitivity of the EoS predictions is explored in terms of $z$ parameter of the SU(3) symmetric model that modifies the meson-hyperon coupling constants away from their SU(6) symmetric values. Our results show that model EoS based on our approach can support static Tolman-Oppenheimer-Volkof (TOV) masses in the range $2.3$-$2.5\,M_{\odot}$ in the large-$Q_{\rm sat}$ and small-$z$ regime, however, such stars contain only a trace amount of hyperons compared to SU(6) models. We also construct uniformly rotating Keplerian configurations for our model EoS for which the masses of stellar sequences may reach up to $3.0\,M_{\odot}$. These results are used to explore the systematic dependence of the ratio of maximum masses of rotating and static stars, the lower bound on the rotational frequency of the models that will allow secondary masses in the gravitational waves events to be compact stars with $M_2 \lesssim 3.0\,M_{\odot}$ and the strangeness fraction on the model parameters. We conclude that very massive stellar models can be, in principle, constructed within the SU(3) symmetric model, however, they are nucleonic-like as their strangeness fraction drops below 3\%., Comment: 11 pages, 7 figures

We discuss the bulk viscosity of hot and dense $npe\mu$ matter arising from weak-interaction direct Urca processes. We consider two regimes of interest: (a) the neutrino-transparent regime with $T\leq T_{\rm tr}$ ($T_{\rm tr}\simeq 5\div 10$ MeV is the neutrino-trapping temperature); and (b) the neutrino-trapped regime with $T\geq T_{\rm tr}$. Nuclear matter is modeled in relativistic density functional approach with density-dependent parametrization DDME2. The maximum of the bulk viscosity is achieved at temperatures $T \simeq 5\div 6$ MeV in the neutrino-transparent regime, then it drops rapidly at higher temperatures where neutrino-trapping occurs. As an astrophysical application, we estimate the damping timescales of density oscillations by the bulk viscosity in neutron star mergers and find that, e.g., at the oscillation frequency $f=10$ kHz, the damping will be very efficient at temperatures $4\leq T\leq 7$ MeV where the bulk viscosity might affect the evolution of the post-merger object., Comment: 16 pages, 9 figures, matches published version

In this work, we consider the consequences of phase transition in dense QCD on the properties of compact stars and implications for the observational program in gravitational wave and X-ray astrophysics. The key underlying assumption of our modeling is a strong first-order phase transition past the point where the hadronic branch of compact stars reaches the two-solar mass limit. Our analysis predicts ultracompact stars with very small radii - in the range of 6-9 km - living on compact star sequences that are entirely consistent with the current multimessenger data. We show that sequences featuring two-solar mass hadronic stars consistent with radio-pulsar observations are also consistent with the inferences of large radii for massive neutron stars by NICER X-ray observations of neutron stars and the small radii predicted by gravitational waves analysis of the binary neutron star inspiral event GW170817 for our models that feature a strong first-order QCD phase transition., Comment: v2: matches published version, 6 pages, 5 figures

The equation of state (EoS) and composition of dense and hot $\Delta$-resonance admixed hypernuclear matter is studied under conditions that are characteristic of neutron star binary merger remnants and supernovas. The cold, neutrino free regime is also considered as a reference for the astrophysical constraints on the EoS of dense matter. Our formalism uses the covariant density functional (CDF) theory successfully adapted to include the full $J^P=1/2^+$ baryon octet and non-strange members of $J^P=3/2^+$ decouplet with density-dependent couplings that have been suitably adjusted to the existing laboratory and astrophysical data. The effect of $\Delta$-resonances at finite temperatures is to soften the EoS of hypernuclear matter at intermediate densities and stiffen it at high densities. At low temperatures, the heavy baryons $\Lambda$, $\Delta^-$,$\Xi^-$, $\Xi^0$ and $\Delta^0$ appear in the given order if the $\Delta$-meson couplings are close to those for the nucleon-meson couplings. As is the case for hyperons, the thresholds of $\Delta$-resonances move to lower densities with the increase of temperature indicating a significant fraction of $\Delta$'s in the low-density subnuclear regime. We find that the $\Delta$-resonances comprise a significant fraction of baryonic matter, of the order of $10\%$ at temperatures of the order of several tens of MeV in the neutrino-trapped regime and, thus, may affect the supernova and binary neutron star dynamics by providing, for example, a new source for neutrino opacity or a new channel for bulk viscosity via the direct Urca processes. The mass-radius relation of isentropic static, spherically symmetric hot compact stars is discussed., Comment: v3: matches published version. v2: Version accepted in Euro. Phys. Journ. A, 15 pages, 11 figures. v1: 12 pages, 7 figures. Contribution to the EPJ A topical issue "CompOSE: a repository for Neutron Star Equations of State and Transport Properties"

Several studies have shown the ability of natural gradient descent to minimize the objective function more efficiently than ordinary gradient descent based methods. However, the bottleneck of this approach for training deep neural networks lies in the prohibitive cost of solving a large dense linear system corresponding to the Fisher Information Matrix (FIM) at each iteration. This has motivated various approximations of either the exact FIM or the empirical one. The most sophisticated of these is KFAC, which involves a Kronecker-factored block diagonal approximation of the FIM. With only a slight additional cost, a few improvements of KFAC from the standpoint of accuracy are proposed. The common feature of the four novel methods is that they rely on a direct minimization problem, the solution of which can be computed via the Kronecker product singular value decomposition technique. Experimental results on the three standard deep auto-encoder benchmarks showed that they provide more accurate approximations to the FIM. Furthermore, they outperform KFAC and state-of-the-art first-order methods in terms of optimization speed.

Several global parameters of compact stars are related via empirical relations, which are (nearly) independent of the underlying equation of state (EoS) of dense matter and, therefore, are said to be universal. We investigate the universality of relations that express the maximum mass and the radius of non-rotating and maximally rapidly rotating configurations, as well as their moment of inertia, in terms of the compactness of the star. For this, we first utilize a collection of cold (zero-temperature) and hot (isentropic) nucleonic EoS and confirm that the universal relations are holding for our collection of EoS. We then go on, to add to our collection and test for the same universality models of EoS that admit a strong first-order phase transition from nucleonic to deconfined quark matter. Also in this case we find that the universal relations hold, in particular for hot, isentropic hybrid stars. By fitting the universal relations to our computed data, we determine the coefficients entering these relations and the accuracy to which they hold., Comment: 16 pages, 8 figures, published in MNRAS, Vol. 515, p. 3539 (2022)

Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations., Comment: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

We present a new derivation of relativistic second-order dissipative hydrodynamics for quantum systems using Zubarev's non-equilibrium statistical-operator formalism. This is achieved by a systematic expansion of the energy-momentum tensor and the charge current to second order in deviations from equilibrium. As a concrete example, we obtain the relaxation equations for the shear-stress tensor, the bulk-viscous pressure, and the charge-diffusion currents required to close the set of equations of motion for relativistic second-order dissipative hydrodynamics. We also identify new transport coefficients which describe the relaxation of dissipative processes to second order and express them in terms of equilibrium correlation functions, thus establishing new Kubo-type formulas for second-order transport coefficients., Comment: 54 pages

Finite-temperature equation of state (EoS) and the composition of dense nuclear and hypernuclear matter under conditions characteristic of neutron star binary merger remnants and supernovas are discussed. We consider both neutrino free-streaming and trapped regimes which are separated by a temperature of a few MeV. The formalism is based on covariant density functional (CDF) theory for the full baryon octet with density-dependent couplings, suitably adjusted in the hypernuclear sector. The softening of the EoS with the introduction of the hyperons is quantified under various conditions of lepton fractions and temperatures. We find that $\Lambda$, $\Xi^-$, and $\Xi^0$ hyperons appear in the given order with a sharp density increase at zero temperature at the threshold being replaced by an extended increment over a wide density range at high temperatures. The $\Lambda$ hyperon survives in the deep subnuclear regime. The triplet of $\Sigma$s is suppressed in cold hypernuclear matter up to around seven times the nuclear saturation density, but appears in significant fractions at higher temperatures, $T\geq 20$ MeV, in both supernova and merger remnant matter. We point out that a special isospin degeneracy point exists where the baryon abundances within each of the three isospin multiplets are equal to each other as a result of (approximate) isospin symmetry. At that point, the charge chemical potential of the system vanishes. We find that under the merger remnant conditions, the fractions of electron and $\mu$-on neutrinos are close and are about 1\%, whereas in the supernova case, we only find a significant fraction ($\sim$10\%) of electron neutrinos, given that in this case, the $\mu$-on lepton number is zero., Comment: v2: 16 pages, 7 figure, minor clarifications, matches published version. v1: 17 pages, 6 figures, submitted to "Universe"

Higher harmonic generation in plane graphene quantum dots initiated by intense coherent radiation is investigated, using dynamical Hartree-Fock mean-field theory. A microscopic theory describing the extreme nonlinear optical response of plane graphene quantum dots is developed. The closed set of differential equations for the single-particle density matrix at the graphene quantum dots-strong laser field multiphoton interaction is solved numerically. The obtained solutions indicate the significance of the type of edge and lateral size, and bandgap and laser field strength in the high harmonic generation process on the triangular graphene quantum dot., Comment: 8 pages, 9 figures

We explore the implications of the recent radius determination of PSR J0740+6620 by the NICER experiment combined with the neutron skin measurement by the PREX-II experiment and the associated inference of the slope of symmetry energy, for the structure of hybrid stars with a strong first-order phase transition from nucleonic to quark matter. We combine a covariant density-functional nucleonic equation of state (EOS) with a constant-speed-of-sound EOS for quark matter. We show that the radius and tidal deformability ranges obtained from GW170817 can be reconciled with the implication of the PREX-II experiment if there is a phase transition to quark matter in the low-mass compact star. In the high-mass segment, the EoS needs to be stiff to comply with the large-radius inference for PSR J0740+6620 and J0030+0451 with masses $M\simeq 2M_{\odot}$ and $M\simeq 1.4M_{\odot}$. We show that twin stars are not excluded, but the mass and radius ranges (with $M \geq M_\odot$) are restricted to narrow domains $\Delta M_{\rm twin} \lesssim 0.05 M_\odot$ and $\Delta R_{\rm twin} \sim 1.0$~km. We also show that the existence of twin configurations is compatible with the light companion in the GW190814 event being a hybrid star in the case of values of the sound-speed square $s=0.6$ and $s=1/3$., Comment: 7 pages and 4 figures; matches published version

In this work, we extend our previous study of the bulk viscosity of hot and dense $npe$ matter induced by the Urca processes in the neutrino trapped regime to $npe\mu$ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. The nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and antineutrino dominated matter (NL3 model) at temperatures for which neutrinos/antineutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca-process-driven bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as $\zeta\sim T^{-2}$ at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range $\gtrsim$ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas., Comment: 26 pages, 13 figures

We study the oscillations modes of differential rotating remnants of binary neutron star inspirals by modeling them as incompressible Riemann ellipsoids parametrized by the ratio $f$ of their internal circulation to the rotation frequency. The effects of viscosity and gravitational wave radiation on the modes are studied and it is shown that these bodies exhibit generic instabilities towards gravitational wave radiation akin to the Chandrasekhar--Friedman--Schutz instabilities for uniformly rotating stars. The odd-parity modes are unstable for all values of $f$ (except for the spherical model) and deformations, whereas the even parity unstable modes appear only in highly eccentric ellipsoids. We quantify the modification of the modes with varying mass of the model and the magnitude of the viscosity. The modes are weakly dependent on the range of the masses relevant to the binary neutron star mergers. Large turbulent viscosity can lead to a suppression of the gravitational wave instabilities, whereas kinematical viscosity has a negligible influence on the modes and their damping timescales., Comment: 16 pages, 11 figures. Submitted to MNRAS

We review the covariant density functional approach to the equation of state of the dense nuclear matter in compact stars. The main emphasis is on the hyperonization of the dense matter, and the role played by the $\Delta$-resonances. The implications of hyperonization for the astrophysics of compact stars, including the equation of state, composition, and stellar parameters are examined. The mass-radius relation and tidal deformabilities of static and rapidly rotating (Keplerian) configurations are discussed in some detail. We briefly touch upon some other recent developments involving hyperonization in hot hypernuclear matter at high- and low-densities., Comment: v3: matches final version, references updated; v2: 47 pages, 11 figures, minor corrections; v1: 45 pages, 11 figures, chapter to appear in "Astrophysics in the XXI century with compact stars", World Scientific, eds. Cesar Zen Vasconcellos and Fridolin Weber

In this work, we study the effect of (anti)kaon condensation on the properties of compact stars that develop hypernuclear cores with and without an admixture of $\Delta$-resonances. We work within the covariant density functional theory with the parameters adjusted to $K$-atomic and kaon-nucleon scattering data in the kaonic sector. The density-dependent parameters in the hyperonic sector are adjusted to the data on $\Lambda$ and $\Xi^-$ hypernuclei data. The $\Delta$-resonance couplings are tuned to the data obtained from their scattering off nuclei and heavy-ion collision experiments. We find that (anti)kaon condensate leads to a softening of the equation of state and lower maximum masses of compact stars than in the absence of the condensate. Both the $K^-$ and $\bar K^0$-condensations occur through a second-order phase transition, which implies no mixed-phase formation. For large values of (anti)kaon and $\Delta$-resonance potentials in symmetric nuclear matter, we observe that condensation leads to an extinction of $ \Xi^{-,0}$ hyperons. We also investigate the influence of inclusion of additional hidden-strangeness $\sigma^{*}$ meson in the functional and find that it leads to a substantial softening of the equation of state and delay in the onset of (anti)kaons., Comment: 14 pages, 11 figures, accepted for publication in Phys. Rev. D

We conjecture and verify a set of universal relations between global parameters of hot and fast-rotating compact stars, including a relation connecting the masses of the mass-shedding (Kepler) and static configurations. We apply these relations to the GW170817 event by adopting the scenario in which a hypermassive compact star remnant formed in a merger evolves into a supramassive compact star that collapses into a black hole once the stability line for such stars is crossed. We deduce an upper limit on the maximum mass of static, cold neutron stars $ 2.15^{+0.10}_{-0.07}\le M^\star_{\mathrm{TOV}} \le 2.24^{+0.12}_{-0.10} $ for the typical range of entropy per baryon $2 \le S/A \le 3$ and electron fraction $Y_e = 0.1$ characterizing the hot hypermassive star. Our result implies that accounting for the finite temperature of the merger remnant relaxes previously derived constraints on the value of the maximum mass of a cold, static compact star., Comment: 17 pages, 11 figures

Binary neutron star mergers, such as the multimessenger GW170817 event, may produce hypermassive compact objects which are supported against collapse by the internal circulation of the fluid within the star. We compute their unstable modes of oscillations driven by gravitational wave radiation and shear viscosity, modeling them as triaxial Riemann ellipsoids. We work in a perturbative regime, where the gravitational radiation reaction force is taken into account at 2.5-post-Newtonian order and find unstable modes with dissipation timescales $\gtrsim 1$ ms which are relevant to the transient state of a hypermassive remnant of a merger. We show that the secular instabilities are dominated by gravitational wave radiation. If the shear viscosity is included, it can increase the growth times or even stabilize the unstable modes, but it must have values several orders of magnitude larger than predicted for cold neutron stars., Comment: 7 pages, 3 figures, 1 table. Accepted by ApJ Letters