@article{SposiniChechkinSokolovetal.2022, author = {Sposini, Vittoria and Chechkin, Aleksei and Sokolov, Igor M. and Roldan-Vargas, Sandalo}, title = {Detecting temporal correlations in hopping dynamics in Lennard-Jones liquids}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {32}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac7e0a}, pages = {15}, year = {2022}, abstract = {Lennard-Jones mixtures represent one of the popular systems for the study of glass-forming liquids. Spatio/temporal heterogeneity and rare (activated) events are at the heart of the slow dynamics typical of these systems. Such slow dynamics is characterised by the development of a plateau in the mean-squared displacement (MSD) at intermediate times, accompanied by a non-Gaussianity in the displacement distribution identified by exponential tails. As pointed out by some recent works, the non-Gaussianity persists at times beyond the MSD plateau, leading to a Brownian yet non-Gaussian regime and thus highlighting once again the relevance of rare events in such systems. Single-particle motion of glass-forming liquids is usually interpreted as an alternation of rattling within the local cage and cage-escape motion and therefore can be described as a sequence of waiting times and jumps. In this work, by using a simple yet robust algorithm, we extract jumps and waiting times from single-particle trajectories obtained via molecular dynamics simulations. We investigate the presence of correlations between waiting times and find negative correlations, which becomes more and more pronounced when lowering the temperature.}, language = {en} } @article{BolotovSmirnovBubnovaetal.2021, author = {Bolotov, Maxim I. and Smirnov, Lev A. and Bubnova, E. S. and Osipov, Grigory V. and Pikovskij, Arkadij}, title = {Spatiotemporal regimes in the Kuramoto-Battogtokh system of nonidentical oscillators}, series = {Journal of experimental and theoretical physics}, volume = {132}, journal = {Journal of experimental and theoretical physics}, number = {1}, publisher = {Springer}, address = {Heidelberg [u.a.]}, issn = {1063-7761}, doi = {10.1134/S1063776121010106}, pages = {127 -- 147}, year = {2021}, abstract = {We consider the spatiotemporal states of an ensemble of nonlocally coupled nonidentical phase oscillators, which correspond to different regimes of the long-term evolution of such a system. We have obtained homogeneous, twisted, and nonhomogeneous stationary solutions to the Ott-Antonsen equations corresponding to key variants of the realized collective rotational motion of elements of the medium in question with nonzero mesoscopic characteristics determining the degree of coherence of the dynamics of neighboring particles. We have described the procedures of the search for the class of nonhomogeneous solutions as stationary points of the auxiliary point map and of determining the stability based on analysis of the eigenvalue spectrum of the composite operator. Static and breather cluster regimes have been demonstrated and described, as well as the regimes with an irregular behavior of averaged complex fields including, in particular, the local order parameter.}, language = {en} } @article{KlettCherstvyShinetal.2021, author = {Klett, Kolja and Cherstvy, Andrey G. and Shin, Jaeoh and Sokolov, Igor M. and Metzler, Ralf}, title = {Non-Gaussian, transiently anomalous, and ergodic self-diffusion of flexible dumbbells in crowded two-dimensional environments}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {104}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {6}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.104.064603}, pages = {18}, year = {2021}, abstract = {We employ Langevin-dynamics simulations to unveil non-Brownian and non-Gaussian center-of-mass self-diffusion of massive flexible dumbbell-shaped particles in crowded two-dimensional solutions. We study the intradumbbell dynamics of the relative motion of the two constituent elastically coupled disks. Our main focus is on effects of the crowding fraction phi and of the particle structure on the diffusion characteristics. We evaluate the time-averaged mean-squared displacement (TAMSD), the displacement probability-density function (PDF), and the displacement autocorrelation function (ACF) of the dimers. For the TAMSD at highly crowded conditions of dumbbells, e.g., we observe a transition from the short-time ballistic behavior, via an intermediate subdiffusive regime, to long-time Brownian-like spreading dynamics. The crowded system of dimers exhibits two distinct diffusion regimes distinguished by the scaling exponent of the TAMSD, the dependence of the diffusivity on phi, and the features of the displacement-ACF. We attribute these regimes to a crowding-induced transition from viscous to viscoelastic diffusion upon growing phi. We also analyze the relative motion in the dimers, finding that larger phi suppress their vibrations and yield strongly non-Gaussian PDFs of rotational displacements. For the diffusion coefficients D(phi) of translational and rotational motion of the dumbbells an exponential decay with phi for weak and a power-law variation D(phi) proportional to (phi - phi(star))(2.4) for strong crowding is found. A comparison of simulation results with theoretical predictions for D(phi) is discussed and some relevant experimental systems are overviewed.}, language = {en} } @article{DoerriesLoosKlapp2021, author = {D{\"o}rries, Timo and Loos, Sarah Anna Marie and Klapp, Sabine H. L.}, title = {Correlation functions of non-Markovian systems out of equilibrium}, series = {Journal of statistical mechanics: theory and experiment : JSTAT}, journal = {Journal of statistical mechanics: theory and experiment : JSTAT}, number = {3}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1742-5468}, doi = {10.1088/1742-5468/abdead}, pages = {36}, year = {2021}, abstract = {This paper is concerned with correlation functions of stochastic systems with memory, a prominent example being a molecule or colloid moving through a complex (e.g. viscoelastic) fluid environment. Analytical investigations of such systems based on non-Markovian stochastic equations are notoriously difficult. A common approximation is that of a single-exponential memory, corresponding to the introduction of one auxiliary variable coupled to the Markovian dynamics of the main variable. As a generalization, we here investigate a class of 'toy' models with altogether three degrees of freedom, giving rise to more complex forms of memory. Specifically, we consider, mainly on an analytical basis, the under- and overdamped motion of a colloidal particle coupled linearly to two auxiliary variables, where the coupling between variables can be either reciprocal or non-reciprocal. Projecting out the auxiliary variables, we obtain non-Markovian Langevin equations with friction kernels and colored noise, whose structure is similar to that of a generalized Langevin equation. For the present systems, however, the non-Markovian equations may violate the fluctuation-dissipation relation as well as detailed balance, indicating that the systems are out of equilibrium. We then study systematically the connection between the coupling topology of the underlying Markovian system and various autocorrelation functions. We demonstrate that already two auxiliary variables can generate surprisingly complex (e.g. non-monotonic or oscillatory) memory and correlation functions. Finally, we show that a minimal overdamped model with two auxiliary variables and suitable non-reciprocal coupling yields correlation functions resembling those describing hydrodynamic backflow in an optical trap.}, language = {en} } @article{CherstvyWangMetzleretal.2021, author = {Cherstvy, Andrey G. and Wang, Wei and Metzler, Ralf and Sokolov, Igor M.}, title = {Inertia triggers nonergodicity of fractional Brownian motion}, series = {Physical review : E, Statistical, nonlinear and soft matter physics}, volume = {104}, journal = {Physical review : E, Statistical, nonlinear and soft matter physics}, number = {2}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0045}, doi = {10.1103/PhysRevE.104.024115}, pages = {12}, year = {2021}, abstract = {How related are the ergodic properties of the over- and underdamped Langevin equations driven by fractional Gaussian noise? We here find that for massive particles performing fractional Brownian motion (FBM) inertial effects not only destroy the stylized fact of the equivalence of the ensemble-averaged mean-squared displacement (MSD) to the time-averaged MSD (TAMSD) of overdamped or massless FBM, but also dramatically alter the values of the ergodicity-breaking parameter (EB). Our theoretical results for the behavior of EB for underdamped or massive FBM for varying particle mass m, Hurst exponent H, and trace length T are in excellent agreement with the findings of stochastic computer simulations. The current results can be of interest for the experimental community employing various single-particle-tracking techniques and aiming at assessing the degree of nonergodicity for the recorded time series (studying, e.g., the behavior of EB versus lag time). To infer FBM as a realizable model of anomalous diffusion for a set single-particle-tracking data when massive particles are being tracked, the EBs from the data should be compared to EBs of massive (rather than massless) FBM.}, language = {en} } @article{GrebenkovMetzlerOshanin2022, author = {Grebenkov, Denis S. and Metzler, Ralf and Oshanin, Gleb}, title = {Search efficiency in the Adam-Delbruck reduction-of-dimensionality scenario versus direct diffusive search}, series = {New journal of physics : the open-access journal for physics}, volume = {24}, journal = {New journal of physics : the open-access journal for physics}, number = {8}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {1367-2630}, doi = {10.1088/1367-2630/ac8824}, pages = {32}, year = {2022}, abstract = {The time instant-the first-passage time (FPT)-when a diffusive particle (e.g., a ligand such as oxygen or a signalling protein) for the first time reaches an immobile target located on the surface of a bounded three-dimensional domain (e.g., a hemoglobin molecule or the cellular nucleus) is a decisive characteristic time-scale in diverse biophysical and biochemical processes, as well as in intermediate stages of various inter- and intra-cellular signal transduction pathways. Adam and Delbruck put forth the reduction-of-dimensionality concept, according to which a ligand first binds non-specifically to any point of the surface on which the target is placed and then diffuses along this surface until it locates the target. In this work, we analyse the efficiency of such a scenario and confront it with the efficiency of a direct search process, in which the target is approached directly from the bulk and not aided by surface diffusion. We consider two situations: (i) a single ligand is launched from a fixed or a random position and searches for the target, and (ii) the case of 'amplified' signals when N ligands start either from the same point or from random positions, and the search terminates when the fastest of them arrives to the target. For such settings, we go beyond the conventional analyses, which compare only the mean values of the corresponding FPTs. Instead, we calculate the full probability density function of FPTs for both scenarios and study its integral characteristic-the 'survival' probability of a target up to time t. On this basis, we examine how the efficiencies of both scenarios are controlled by a variety of parameters and single out realistic conditions in which the reduction-of-dimensionality scenario outperforms the direct search.}, language = {en} } @article{MolinaGarciaSandevSafdarietal.2018, author = {Molina-Garcia, Daniel and Sandev, Trifce and Safdari, Hadiseh and Pagnini, Gianni and Chechkin, Aleksei V. and Metzler, Ralf}, title = {Crossover from anomalous to normal diffusion}, series = {New Journal of Physics}, volume = {20}, journal = {New Journal of Physics}, publisher = {IOP Publishing Ltd}, address = {London und Bad Honnef}, issn = {1367-2630}, doi = {10.1088/1367-2630/aae4b2}, pages = {28}, year = {2018}, abstract = {Abstract The emerging diffusive dynamics in many complex systems show a characteristic crossover behaviour from anomalous to normal diffusion which is otherwise fitted by two independent power-laws. A prominent example for a subdiffusive-diffusive crossover are viscoelastic systems such as lipid bilayer membranes, while superdiffusive-diffusive crossovers occur in systems of actively moving biological cells. We here consider the general dynamics of a stochastic particle driven by so-called tempered fractional Gaussian noise, that is noise with Gaussian amplitude and power-law correlations, which are cut off at some mesoscopic time scale. Concretely we consider such noise with built-in exponential or power-law tempering, driving an overdamped Langevin equation (fractional Brownian motion) and fractional Langevin equation motion. We derive explicit expressions for the mean squared displacement and correlation functions, including different shapes of the crossover behaviour depending on the concrete tempering, and discuss the physical meaning of the tempering. In the case of power-law tempering we also find a crossover behaviour from faster to slower superdiffusion and slower to faster subdiffusion. As a direct application of our model we demonstrate that the obtained dynamics quantitatively describes the subdiffusion-diffusion and subdiffusion-subdiffusion crossover in lipid bilayer systems. We also show that a model of tempered fractional Brownian motion recently proposed by Sabzikar and Meerschaert leads to physically very different behaviour with a seemingly paradoxical ballistic long time scaling.}, language = {en} } @article{FeldmannMaduarSanteretal.2016, author = {Feldmann, David and Maduar, Salim R. and Santer, Mark and Lomadze, Nino and Vinogradova, Olga I. and Santer, Svetlana}, title = {Manipulation of small particles at solid liquid interface}, series = {Scientific reports}, volume = {6}, journal = {Scientific reports}, publisher = {Nature Publishing Group}, address = {London}, issn = {2045-2322}, doi = {10.1038/srep36443}, pages = {10}, year = {2016}, abstract = {The strong adhesion of sub-micron sized particles to surfaces is a nuisance, both for removing contaminating colloids from surfaces and for conscious manipulation of particles to create and test novel micro/nano-scale assemblies. The obvious idea of using detergents to ease these processes suffers from a lack of control: the action of any conventional surface-modifying agent is immediate and global. With photosensitive azobenzene containing surfactants we overcome these limitations. Such photo-soaps contain optical switches (azobenzene molecules), which upon illumination with light of appropriate wavelength undergo reversible trans-cis photo-isomerization resulting in a subsequent change of the physico-chemical molecular properties. In this work we show that when a spatial gradient in the composition of trans- and cis- isomers is created near a solid-liquid interface, a substantial hydrodynamic flow can be initiated, the spatial extent of which can be set, e.g., by the shape of a laser spot. We propose the concept of light induced diffusioosmosis driving the flow, which can remove, gather or pattern a particle assembly at a solid-liquid interface. In other words, in addition to providing a soap we implement selectivity: particles are mobilized and moved at the time of illumination, and only across the illuminated area.}, language = {en} } @article{EkhtiariAgarwalMarwanetal.2019, author = {Ekhtiari, Nikoo and Agarwal, Ankit and Marwan, Norbert and Donner, Reik Volker}, title = {Disentangling the multi-scale effects of sea-surface temperatures on global precipitation}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {29}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {6}, publisher = {American Institute of Physics}, address = {Melville}, issn = {1054-1500}, doi = {10.1063/1.5095565}, pages = {12}, year = {2019}, abstract = {The oceans and atmosphere interact via a multiplicity of feedback mechanisms, shaping to a large extent the global climate and its variability. To deepen our knowledge of the global climate system, characterizing and investigating this interdependence is an important task of contemporary research. However, our present understanding of the underlying large-scale processes is greatly limited due to the manifold interactions between essential climatic variables at different temporal scales. To address this problem, we here propose to extend the application of complex network techniques to capture the interdependence between global fields of sea-surface temperature (SST) and precipitation (P) at multiple temporal scales. For this purpose, we combine time-scale decomposition by means of a discrete wavelet transform with the concept of coupled climate network analysis. Our results demonstrate the potential of the proposed approach to unravel the scale-specific interdependences between atmosphere and ocean and, thus, shed light on the emerging multiscale processes inherent to the climate system, which traditionally remain undiscovered when investigating the system only at the native resolution of existing climate data sets. Moreover, we show how the relevant spatial interdependence structures between SST and P evolve across time-scales. Most notably, the strongest mutual correlations between SST and P at annual scale (8-16 months) concentrate mainly over the Pacific Ocean, while the corresponding spatial patterns progressively disappear when moving toward longer time-scales. Published under license by AIP Publishing.}, language = {en} } @article{ToetzkeKardjilovLenoiretal.2019, author = {T{\"o}tzke, Christian and Kardjilov, Nikolay and Lenoir, Nicolas and Manke, Ingo and Oswald, Sascha and Tengattini, Alessandro}, title = {What comes NeXT?}, series = {Optics express : the international electronic journal of optics}, volume = {27}, journal = {Optics express : the international electronic journal of optics}, number = {20}, publisher = {Optical Society of America}, address = {Washington}, issn = {1094-4087}, doi = {10.1364/OE.27.028640}, pages = {28640 -- 28648}, year = {2019}, abstract = {Here, we report on a new record in the acquisition time for fast neutron tomography. With an optimized imaging setup, it was possible to acquire single radiographic projection images with 10 ms and full tomographies with 155 projections images and a physical spatial resolution of 200 mu m within 1.5 s. This is about 6.7 times faster than the current record. We used the technique to investigate the water infiltration in the soil with a living lupine root system. The fast imaging setup will be part of the future NeXT instrument at ILL in Grenoble with a great field of possible future applications. (C) 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement}, language = {en} } @article{KraemerMarwan2019, author = {Kr{\"a}mer, Hauke Kai and Marwan, Norbert}, title = {Border effect corrections for diagonal line based recurrence quantification analysis measures}, series = {Modern physics letters : A, Particles and fields, gravitation, cosmology, nuclear physics}, volume = {383}, journal = {Modern physics letters : A, Particles and fields, gravitation, cosmology, nuclear physics}, number = {34}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0375-9601}, doi = {10.1016/j.physleta.2019.125977}, pages = {16}, year = {2019}, abstract = {Recurrence Quantification Analysis (RQA) defines a number of quantifiers, which base upon diagonal line structures in the recurrence plot (RP). Due to the finite size of an RP, these lines can be cut by the borders of the RP and, thus, bias the length distribution of diagonal lines and, consequently, the line based RQA measures. In this letter we investigate the impact of the mentioned border effects and of the thickening of diagonal lines in an RP (caused by tangential motion) on the estimation of the diagonal line length distribution, quantified by its entropy. Although a relation to the Lyapunov spectrum is theoretically expected, the mentioned entropy yields contradictory results in many studies. Here we summarize correction schemes for both, the border effects and the tangential motion and systematically compare them to methods from the literature. We show that these corrections lead to the expected behavior of the diagonal line length entropy, in particular meaning zero values in case of a regular motion and positive values for chaotic motion. Moreover, we test these methods under noisy conditions, in order to supply practical tools for applied statistical research.}, language = {en} } @article{AgarwalGuntuBanerjeeetal.2022, author = {Agarwal, Ankit and Guntu, Ravikumar and Banerjee, Abhirup and Gadhawe, Mayuri Ashokrao and Marwan, Norbert}, title = {A complex network approach to study the extreme precipitation patterns in a river basin}, series = {Chaos : an interdisciplinary journal of nonlinear science}, volume = {32}, journal = {Chaos : an interdisciplinary journal of nonlinear science}, number = {1}, publisher = {American Institute of Physics}, address = {Woodbury, NY}, issn = {1054-1500}, doi = {10.1063/5.0072520}, pages = {12}, year = {2022}, abstract = {The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.}, language = {en} } @article{ClarkFreyseYashinaetal.2022, author = {Clark, Oliver J. and Freyse, Friedrich and Yashina, L. V. and Rader, Oliver and Sanchez-Barriga, Jaime}, title = {Robust behavior and spin-texture stability of the topological surface state in Bi2Se3 upon deposition of gold}, series = {npj quantum materials}, volume = {7}, journal = {npj quantum materials}, number = {1}, publisher = {Nature Publishing Group}, address = {London}, issn = {2397-4648}, doi = {10.1038/s41535-022-00443-9}, pages = {7}, year = {2022}, abstract = {The Dirac point of a topological surface state (TSS) is protected against gapping by time-reversal symmetry. Conventional wisdom stipulates, therefore, that only through magnetisation may a TSS become gapped. However, non-magnetic gaps have now been demonstrated in Bi2Se3 systems doped with Mn or In, explained by hybridisation of the Dirac cone with induced impurity resonances. Recent photoemission experiments suggest that an analogous mechanism applies even when Bi2Se3 is surface dosed with Au. Here, we perform a systematic spin- and angle-resolved photoemission study of Au-dosed Bi2Se3. Although there are experimental conditions wherein the TSS appears gapped due to unfavourable photoemission matrix elements, our photon-energy-dependent spectra unambiguously demonstrate the robustness of the Dirac cone against high Au coverage. We further show how the spin textures of the TSS and its accompanying surface resonances remain qualitatively unchanged following Au deposition, and discuss the mechanism underlying the suppression of the spectral weight.}, language = {en} } @article{PadashSandevKantzetal.2022, author = {Padash, Amin and Sandev, Trifce and Kantz, Holger and Metzler, Ralf and Chechkin, Aleksei}, title = {Asymmetric Levy flights are more efficient in random search}, series = {Fractal and fractional}, volume = {6}, journal = {Fractal and fractional}, number = {5}, publisher = {MDPI}, address = {Basel}, issn = {2504-3110}, doi = {10.3390/fractalfract6050260}, pages = {23}, year = {2022}, abstract = {We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta. We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times. We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient. When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity.}, language = {en} } @article{GrebenkovKumar2022, author = {Grebenkov, Denis S. and Kumar, Aanjaneya}, title = {First-passage times of multiple diffusing particles with reversible target-binding kinetics}, series = {Journal of physics : A, Mathematical and theoretical}, volume = {55}, journal = {Journal of physics : A, Mathematical and theoretical}, number = {32}, publisher = {IOP Publ.}, address = {Bristol}, issn = {1751-8113}, doi = {10.1088/1751-8121/ac7e91}, pages = {33}, year = {2022}, abstract = {We investigate a class of diffusion-controlled reactions that are initiated at the time instance when a prescribed number K among N particles independently diffusing in a solvent are simultaneously bound to a target region. In the irreversible target-binding setting, the particles that bind to the target stay there forever, and the reaction time is the Kth fastest first-passage time to the target, whose distribution is well-known. In turn, reversible binding, which is common for most applications, renders theoretical analysis much more challenging and drastically changes the distribution of reaction times. We develop a renewal-based approach to derive an approximate solution for the probability density of the reaction time. This approximation turns out to be remarkably accurate for a broad range of parameters. We also analyze the dependence of the mean reaction time or, equivalently, the inverse reaction rate, on the main parameters such as K, N, and binding/unbinding constants. Some biophysical applications and further perspectives are briefly discussed.}, language = {en} } @article{KoelschDietrichUjevicetal.2022, author = {K{\"o}lsch, Maximilian and Dietrich, Tim and Ujevic, Maximiliano and Br{\"u}gmann, Bernd}, title = {Investigating the mass-ratio dependence of the prompt-collapse threshold with numerical-relativity simulations}, series = {Physical review : D, Particles, fields, gravitation, and cosmology}, volume = {106}, journal = {Physical review : D, Particles, fields, gravitation, and cosmology}, number = {4}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0010}, doi = {10.1103/PhysRevD.106.044026}, pages = {27}, year = {2022}, abstract = {The next observing runs of advanced gravitational-wave detectors will lead to a variety of binary neutron star detections and numerous possibilities for multimessenger observations of binary neutron star systems. In this context a clear understanding of the merger process and the possibility of prompt black hole formation after merger is important, as the amount of ejected material strongly depends on the merger dynamics. These dynamics are primarily affected by the total mass of the binary, however, the mass ratio also influences the postmerger evolution. To determine the effect of the mass ratio, we investigate the parameter space around the prompt-collapse threshold with a new set of fully relativistic simulations. The simulations cover three equations of state and seven mass ratios in the range of 1.0 <= q <= 1.75, with five to seven simulations of binary systems of different total mass in each case. The threshold mass is determined through an empirical relation based on the collapse time, which allows us to investigate effects of the mass ratio on the threshold mass and also on the properties of the remnant system. Furthermore, we model effects of mass ratio and equation of state on tidal parameters of threshold configurations.}, language = {en} } @article{EmmaSchianchiPannaraleetal.2022, author = {Emma, Mattia and Schianchi, Federico and Pannarale, Francesco and Sagun, Violetta and Dietrich, Tim}, title = {Numerical simulations of dark matter admixed neutron star binaries}, series = {Particles}, volume = {5}, journal = {Particles}, number = {3}, publisher = {MDPI}, address = {Basel}, issn = {2571-712X}, doi = {10.3390/particles5030024}, pages = {273 -- 286}, year = {2022}, abstract = {Multi-messenger observations of compact binary mergers provide a new way to constrain the nature of dark matter that may accumulate in and around neutron stars. In this article, we extend the infrastructure of our numerical-relativity code BAM to enable the simulation of neutron stars that contain an additional mirror dark matter component. We perform single star tests to verify our code and the first binary neutron star simulations of this kind. We find that the presence of dark matter reduces the lifetime of the merger remnant and favors a prompt collapse to a black hole. Furthermore, we find differences in the merger time for systems with the same total mass and mass ratio, but different amounts of dark matter. Finally, we find that electromagnetic signals produced by the merger of binary neutron stars admixed with dark matter are very unlikely to be as bright as their dark matter-free counterparts. Given the increased sensitivity of multi-messenger facilities, our analysis gives a new perspective on how to probe the presence of dark matter.}, language = {en} } @article{DudiAdhikariBruegmannetal.2022, author = {Dudi, Reetika and Adhikari, Ananya and Br{\"u}gmann, Bernd and Dietrich, Tim and Hayashi, Kota and Kawaguchi, Kyohei and Kiuchi, Kenta and Kyutoku, Koutarou and Shibata, Masaru and Tichy, Wolfgang}, title = {Investigating GW190425 with numerical-relativity simulations}, series = {Physical review : D, Particles, fields, gravitation, and cosmology}, volume = {106}, journal = {Physical review : D, Particles, fields, gravitation, and cosmology}, number = {8}, publisher = {American Physical Society}, address = {College Park}, issn = {2470-0010}, doi = {10.1103/PhysRevD.106.084039}, pages = {11}, year = {2022}, abstract = {The third observing run of the LIGO-Virgo Collaboration has resulted in many gravitational wave detections, including the binary neutron star merger GW190425. However, none of these events have been accompanied with an electromagnetic transient found during extensive follow-up searches. In this article, we perform new numerical-relativity simulations of binary neutron star and black hole-neutron star systems that have a chirp mass consistent with GW190425. Assuming that the GW190425's sky location was covered with sufficient accuracy during the electromagnetic follow-up searches, we investigate whether the nondetection of the kilonova is compatible with the source parameters estimated through the gravitational -wave analysis and how one can use this information to place constraints on the properties of the system. Our simulations suggest that GW190425 is incompatible with an unequal mass binary neutron star merger with a mass ratio q < 0.8 when considering stiff or moderately stiff equations of state if the binary was face on and covered by the observation. Our analysis shows that a detailed observational result for kilonovae will be useful to constrain the mass ratio of binary neutron stars in future events.}, language = {en} } @article{GiegSchianchiDietrichetal.2022, author = {Gieg, Henrique and Schianchi, Federico and Dietrich, Tim and Ujevic, Maximiliano}, title = {Incorporating a Radiative Hydrodynamics Scheme in the Numerical-Relativity Code BAM}, series = {Universe : open access journal}, volume = {8}, journal = {Universe : open access journal}, number = {7}, publisher = {MDPI}, address = {Basel}, issn = {2218-1997}, doi = {10.3390/universe8070370}, pages = {25}, year = {2022}, abstract = {To study binary neutron star systems and to interpret observational data such as gravitational-wave and kilonova signals, one needs an accurate description of the processes that take place during the final stages of the coalescence, for example, through numerical-relativity simulations. In this work, we present an updated version of the numerical-relativity code BAM in order to incorporate nuclear-theory-based equations of state and a simple description of neutrino interactions through a neutrino leakage scheme. Different test simulations, for stars undergoing a neutrino-induced gravitational collapse and for binary neutron stars systems, validate our new implementation. For the binary neutron stars systems, we show that we can evolve stably and accurately distinct microphysical models employing the different equations of state: SFHo, DD2, and the hyperonic BHB Lambda phi. Overall, our test simulations have good agreement with those reported in the literature.}, language = {en} } @article{MatzkaSiddiquiLilienkampetal.2017, author = {Matzka, J{\"u}rgen and Siddiqui, Tarique Adnan and Lilienkamp, Henning and Stolle, Claudia and Veliz, Oscar}, title = {Quantifying solar flux and geomagnetic main field influence on the equatorial ionospheric current system at the geomagnetic observatory Huancayo}, series = {Journal of Atmospheric and Solar-Terrestrial Physics}, volume = {163}, journal = {Journal of Atmospheric and Solar-Terrestrial Physics}, publisher = {Elsevier}, address = {Oxford}, issn = {1364-6826}, doi = {10.1016/j.jastp.2017.04.014}, pages = {120 -- 125}, year = {2017}, abstract = {In order to analyse the sensitivity of the equatorial ionospheric current system, i.e. the solar quiet current system and the equatorial electrojet, to solar cycle variations and to the secular variation of the geomagnetic main field, we have analysed 51 years (1935-1985) of geomagnetic observatory data from Huancayo, Peru. This period is ideal to analyse the influence of the main field strength on the amplitude of the quiet daily variation, since the main field decreases significantly from 1935 to 1985, while the distance of the magnetic equator to the observatory remains stable. To this end, we digitised some 19 years of hourly mean values of the horizontal component (H), which have not been available digitally at the World Data Centres. Then, the sensitivity of the amplitude Ali of the quiet daily variation to both solar cycle variations (in terms of sunspot numbers and solar flux F10.7) and changes of the geomagnetic main field strength (due to secular variation) was determined. We confirm an increase of Delta H for the decreasing main field in this period, as expected from physics based models (Cnossen, 2016), but with a somewhat smaller rate of 4.4\% (5.8\% considering one standard error) compared with 6.9\% predicted by the physics based model.}, language = {en} }