Institutes
Refine
Year of publication
Document Type
- Article (401)
- Doctoral Thesis (222)
- Preprint (43)
- Conference Proceeding (17)
- Master's Thesis (12)
- Bachelor Thesis (9)
- Contribution to a Periodical (9)
- Book (5)
- Habilitation (3)
- Other (1)
Has Fulltext
- yes (724)
Is part of the Bibliography
- no (724) (remove)
Keywords
- Relativistic heavy-ion collisions (6)
- ALICE (5)
- Black holes (5)
- Heavy-ion collisions (5)
- LHC (5)
- Toroidales Magnetfeld (5)
- Equation of state (4)
- Fluctuations (4)
- HADES (4)
- QCD (4)
Institute
Results on proton and Λ flow, calculated with the UrQMD model that incorporates different realistic density dependent equations of state, are presented. It is shown that the proton and hyperon flow shows sensitivity to the equation of state and especially to the appearance of a phase transition at densities below 4n0. Even though qualitatively hyperons and protons exhibit the same beam energy dependence of the flow, the quantitative results are different. In this context it is suggested that the hyperon measurements can be used to study the density dependence of the hyperon interaction in high density QCD matter.
A single wavelength heterodyne interferometer has been set up to investigate the free electron density integrated axially along the line of sight (line density) in a theta-pinch plasma to determine its applicability as a plasma target for ion beam stripping. The maximal line density reached in this experiment was (3.57 ± 0.28) × 1018 cm−2 at 80 Pa and 20 kV. The findings demonstrate the pulsed character of the line density and its increase by raising the load voltage and the working gas pressure. Additionally, the results were compared with spectroscopic free electron density estimations, which were carried out by Hβ -line broadening and peak separation. The time behavior of the line density indicates that its peak value is delayed by about 10 μs compared to the spectroscopic results. This effect is due to the formation of an extended, magnetically compressed plasma column in the vicinity of the current maximum, although the highest volumetric free electron density is reached near the current zero crossing. Since the line density is an essential parameter in describing the stripping capabilities of the plasma target, the interferometric diagnostic is superior to a spectroscopic diagnostic, because it directly provides integrated values along the line of sight. Furthermore, the measurements of the line density in this experiment partially show nonphysical negative values, which is due to gaseous effects and residual shot vibrations.
This thesis aims to investigate the properties of hadronic matter by analyzing fluctuations of conserved charges. A transport model (SMASH) is used for these studies to achieve this. The first part of this thesis focuses on examining transport coefficients, specifically the diffusion coefficients of conserved charges and the shear viscosity. The second part investigates equal-time correlations of particle numbers in the form of cumulants. The last chapter studies different aspects of the isobar collision systems Ru and Zr.
As a first step, the hadronic medium and interactions between its constituents are introduced, and simultaneously, their impact on transport coefficients is investigated. The methodology is verified by comparing the results of SMASH with Chapman-Enskog calculations, followed by examining 3-to-1 multi-particle reactions, revealing their influence on shear viscosity and electrical diffusion. The analysis of the full hadron gas considers angle-dependent cross-sections and additional elastic cross-sections via the AQM description, showing significant impacts on transport coefficients. The dependency on the number of degrees of freedom is explored, with noticeable effects on diffusion coefficients but a smaller influence on the shear viscosity. At non-zero baryon chemical potential, the diffusion coefficients are strongly influenced, while the shear viscosity remains unaffected. Overall, the study underscores the importance of individual cross-sections and the modeling of interactions on transport coefficients.
The following chapter explores fluctuations of conserved charges, crucial for understanding phase transitions in heavy-ion collision from the quark-gluon plasma to the hadronic phase. Using SMASH, the impact of global charge conservation on particle number cumulants in subvolumes of boxes simulating infinite matter is studied. Comparisons with simpler systems highlights the influence of hadronic interactions on cumulants, especially via charge annihilation processes and the results from SMASH shows agreement with analytical calculations. Calculations at finite baryon chemical potential reveals a transition from a Poisson to Skellam distribution within the net proton cumulants. It is shown that an unfolding procedure to obtain the net baryon fluctuations from the net proton ones deviates from the actual net baryon result, particularly in larger volumes. Finally, net proton correlations at vanishing baryon chemical potential align with ALICE measurements and the net proton cumulants are unaffected by deuteron formation.
In the next step, the goal is to investigate critical fluctuations in the hadronic medium. Therefore, the hadronic system is initialized with critical equilibrium fluctuations by coupling the hadron resonance gas with the 3D Ising model. The single-particle probability distributions are derived from the principle of maximum entropy. Evolving these distributions in SMASH, their development in an expanding sphere adjusted to experimental conditions can be analyzed. It reveals resonance decay and formations as the primary source that affects the particle cumulants. Because of isospin randomization processes, critical fluctuations are better preserved in net nucleon numbers. However, for the strongest coupling investigated in this work, correlations of the critical field are still present in the final state of the evolution in the net proton fluctuations. Examining cumulant dependence on rapidity windows shows a non-monotonic trend.
In the third part, collisions involving the isobars Ru and Zr are studied at a center-of-mass energy of 200 GeV. Initially, SMASH is used to study the initial conditions to hydrodynamical simulations, emphasizing the importance of the nuclear structure of isobars on the geometry of the collision area. It is found that the deformation parameters notably influence the initial state. Correlations between nucleon-nucleon pairs on eccentricity fluctuations yield no significant effect. Subsequently, the hydrodynamic model vHLLE evolves the previously explored initial conditions and for the transition between the hydrodynamic and kinetic descriptions, the Cooper-Frye formula is used. Usage of the canonical ensemble ensures the exact conservation of the conserved charges B, Q, and S. The neutron skin effect, which changes the charge distribution within Ru nuclei, is additionally considered. Fluctuations are assessed, revealing suppression in large rapidity windows due to global charge conservation. The hadronic phase modifies fluctuations of net pions, net kaons, and net protons via annihilation processes, yet fluctuations remain unaffected by the neutron skin effect.
The equation of state (EoS) of matter at extremely high temperatures and densities is currently not fully understood, and remains a major challenge in the field of nuclear physics. Neutron stars harbor such extreme conditions and therefore serve as celestial laboratories for constraining the dense matter EoS. In this thesis, we present a novel algorithm that utilizes the idea of Bayesian analysis and the computational efficiency of neural networks to reconstruct the dense matter equation of state from mass-radius observations of neutron stars. We show that the results are compatible with those from earlier works based on conventional methods, and are in agreement with the limits on tidal deformabilities obtained from the gravitational wave event, GW170817. We also observe that the resulting squared speed of sound from the reconstructed EoS features a peak, indicating a likely convergence to the conformal limit at asymptotic densities, as expected from quantum chromodynamics. The novel algorithm can also be applied across various fields faced with computational challenges in solving inverse problems. We further examine the efficiency of deep learning methods for analyzing gravitational waves from compact binary coalescences in this thesis. In particular, we develop a deep learning classifier to segregate simulated gravitational wave data into three classes: signals from binary black hole mergers, signals from binary neutron star mergers, or white noise without any signals. A second deep learning algorithm allows for the regression of chirp mass and combined tidal deformability from simulated binary neutron star mergers. An accurate estimation of these parameters is crucial to constrain the underlying EoS. Lastly, we explore the effects of finite temperatures on the binary neutron star merger remnant from GW170817. Isentropic EoSs are used to infer the frequencies of the rigidly rotating remnant and are noted to be significantly lower compared to previous estimates from zero temperature EoSs. Overall, this thesis presents novel deep learning methods to constrain the neutron star EoS, which will prove useful in future, as more observational data is expected in the upcoming years.
The appearance of strangeness in the form of hyperons within the inner core of neutron stars is expected to affect its detectable properties, such as its global structure or gravitational wave emission. This work explores the parameter space of hyperonic stars within the framework of the Relativistic Mean Field model allowed by the present uncertainties in the state-of-the-art nuclear and hypernuclear experimental data. We impose multi-physics constraints at different density regimes to restrict the parameter space: Chiral effective field theory, heavy-ion collision data, and multi-messenger astrophysical observations of neutron stars. We investigate possible correlations between empirical nuclear and hypernuclear parameters, particularly the symmetry energy and its slope, with observable properties of neutron stars. We do not find a correlation for the hyperon parameters and the astrophysical data. However, the inclusion of hyperons generates a tension between the astrophysical and heavy-ion data constraining considerably the available parameter space.
Based on the positive results of the 0.63 m unmodulated 325 MHz Ladder-RFQ prototype from 2013 to 2016 [1, 2], a modulated 3.3m Ladder-RFQ (s. Fig. 1) has been designed and built for the acceleration of up to 100 mA protons from 95 keV to 3.0 MeV at the FAIR p-Linac [3, 4]. In this paper, we will show the results of manufacturing as well as low level RF measurements of the Ladder-RFQ including flatness and frequency tuning.
Determining the sound speed cs in compact stars is an important open question with numerous implications on the behavior of matter at large densities and hence on gravitational-wave emission from neutron stars. To this scope, we construct more than 107 equations of state (EOSs) with continuous sound speed and build more than 108 nonrotating stellar models consistent not only with nuclear theory and perturbative QCD, but also with astronomical observations. In this way, we find that EOSs with subconformal sound speeds, i.e., with cs 1 3 2 < within the stars, are possible in principle but very unlikely in practice, being only 0.03% of our sample. Hence, it is natural to expect that cs 1 3 2 > somewhere in the stellar interior. Using our large sample, we obtain estimates at 95% credibility of neutron-star radii for representative stars with 1.4 and 2.0 solar masses, R1.4 12.42 km 0.99 0.52 = - + , R2.0 12.12 km 1.23 1.11 = - + , and for the binary tidal deformability of the GW170817 event, 1.186 485 211 225 L = - ˜ + . Interestingly, our lower bounds on the radii are in very good agreement with the prediction derived from very different arguments, namely, the threshold mass. Finally, we provide simple analytic expressions to determine the minimum and maximum values of L˜ as a function of the chirp mass.
Using more than a million randomly generated equations of state that satisfy theoretical and observational constraints, we construct a novel, scale-independent description of the sound speed in neutron stars, where the latter is expressed in a unit cube spanning the normalized radius, r/R, and the mass normalized to the maximum one, M/MTOV. From this generic representation, a number of interesting and surprising results can be deduced. In particular, we find that light (heavy) stars have stiff (soft) cores and soft (stiff) outer layers, or that the maximum of the sound speed is located at the center of light stars but moves to the outer layers for stars with M/MTOV ≳ 0.7, reaching a constant value of cs = 1 2 2 as M → MTOV. We also show that the sound speed decreases below the conformal limit cs = 1 3 2 at the center of stars with M = MTOV. Finally, we construct an analytic expression that accurately describes the radial dependence of the sound speed as a function of the neutron-star mass, thus providing an estimate of the maximum sound speed expected in a neutron star.
We have investigated the systematic differences introduced when performing a Bayesian-inference analysis of the equation of state (EOS) of neutron stars employing either variable- or constant-likelihood functions. The former has the advantage of retaining the full information on the distributions of the measurements, making exhaustive usage of the data. The latter, on the other hand, has the advantage of a much simpler implementation and reduced computational costs. In both approaches, the EOSs have identical priors and have been built using the sound speed parameterization method so as to satisfy the constraints from X-ray and gravitational waves observations, as well as those from chiral effective theory and perturbative quantum chromodynamics. In all cases, the two approaches lead to very similar results and the 90% confidence levels essentially overlap. Some differences do appear, but in regions where the probability density is extremely small and are mostly due to the sharp cutoff on the binary tidal deformability L˜ 720 set in the constant-likelihood approach. Our analysis has also produced two additional results. First, an inverse correlation between the normalized central number density, nc,TOV/ns, and the radius of a maximally massive star, RTOV. Second, and most importantly, it has confirmed the relation between the chirp mass and the binary tidal deformability. The importance of this result is that it relates chirp, which is measured very accurately, and L˜ , which contains important information on the EOS. Hence, when chirp is measured in future detections, our relation can be used to set tight constraints on L˜ .
A considerable effort has been dedicated recently to the construction of generic equations of state (EOSs) for matter in neutron stars. The advantage of these approaches is that they can provide model-independent information on the interior structure and global properties of neutron stars. Making use of more than 106 generic EOSs, we assess the validity of quasi-universal relations of neutron-star properties for a broad range of rotation rates, from slow rotation up to the mass-shedding limit. In this way, we are able to determine with unprecedented accuracy the quasi-universal maximum-mass ratio between rotating and nonrotating stars and reveal the existence of a new relation for the surface oblateness, i.e., the ratio between the polar and equatorial proper radii. We discuss the impact that our findings have on the imminent detection of new binary neutron-star mergers and how they can be used to set new and more stringent limits on the maximum mass of nonrotating neutron stars, as well as to improve the modeling of the X-ray emission from the surface of rotating stars.
The amplification of magnetic fields plays an important role in explaining numerous astrophysical phenomena associated with binary neutron star mergers, such as mass ejection and the powering of short gamma-ray bursts. Magnetic fields in isolated neutron stars are often assumed to be confined to a small region near the stellar surface, while they are normally taken to fill the whole star in numerical modeling of mergers. By performing high-resolution, global, and high-order general-relativistic magnetohydrodynamic simulations, we investigate the impact of a purely crustal magnetic field and contrast it with the standard configuration consisting of a dipolar magnetic field with the same magnetic energy but filling the whole star. While the crust configurations are very effective in generating strong magnetic fields during the Kelvin–Helmholtz-instability stage, they fail to achieve the same level of magnetic-field amplification of the full-star configurations. This is due to the lack of magnetized material in the neutron-star interiors to be used for further turbulent amplification and to the surface losses of highly magnetized matter in the crust configurations. Hence, the final magnetic energies in the two configurations differ by more than 1 order of magnitude. We briefly discuss the impact of these results on astrophysical observables and how they can be employed to deduce the magnetic topology in merging binaries.
Post-merger gravitational-wave signal from neutron-star binaries: a new look at an old problem
(2023)
The spectral properties of the post-merger gravitational-wave signal from a binary of neutron stars encodes a variety of information about the features of the system and of the equation of state describing matter around and above nuclear saturation density. Characterizing the properties of such a signal is an “old” problem, which first emerged when a number of frequencies were shown to be related to the properties of the binary through “quasiuniversal” relations. Here we take a new look at this old problem by computing the properties of the signal in terms of the Weyl scalar ψ4. In this way, and using a database of more than 100 simulations, we provide the first evidence for a new instantaneous frequency, y f0 4, associated with the instant of quasi-time-symmetry in the dynamics, and which also follows a quasi-universal relation. We also derive a new quasi-universal relation for the merger frequency f h mer, which provides a description of the data that is 4 times more accurate than previous expressions while requiring fewer fitting coefficients. Finally, consistent with the findings of numerous studies before ours, and using an enlarged ensemble of binary systems, we point out that the ℓ = 2, m = 1 gravitational-wave mode could become comparable with the traditional ℓ = 2, m = 2 mode on sufficiently long timescales, with strain amplitudes in a ratio |h21|/|h22| ∼ 0.1–1 under generic orientations of the binary, which could be measured by present detectors for signals with a large signal-to-noise ratio or by third-generation detectors for generic signals should no collapse occur.
We use holography to study the dynamics of a strongly-coupled gauge theory in four-dimensional de Sitter space with Hubble rate H. The gauge theory is non-conformal with a characteristic mass scale M. We solve Einstein’s equations numerically and determine the time evolution of homogeneous gauge theory states. If their initial energy density is high compared with H4 then the early-time evolution is well described by viscous hydrodynamics with a non-zero bulk viscosity. At late times the dynamics is always far from equilibrium. The asymptotic late-time state preserves the full de Sitter symmetry group and its dual geometry is a domain-wall in AdS5. The approach to this state is characterised by an emergent relation of the form P = w ℰ that is different from the equilibrium equation of state in flat space. The constant w does not depend on the initial conditions but only on H/M and is negative if the ratio H/M is close to unity. The event and the apparent horizons of the late-time solution do not coincide with one another, reflecting its non-equilibrium nature. In between them lies an “entanglement horizon” that cannot be penetrated by extremal surfaces anchored at the boundary, which we use to compute the entanglement entropy of boundary regions. If the entangling region equals the observable universe then the extremal surface coincides with a bulk cosmological horizon that just touches the event horizon, while for larger regions the extremal surface probes behind the event horizon.
One-photon and multi-photon absorption, spontaneous and stimulated photon emission, resonance Raman scattering and electron transfer are important molecular processes that commonly involve combined vibrational-electronic (vibronic) transitions. The corresponding vibronic transition profiles in the energy domain are usually determined by Franck-Condon factors (FCFs), the squared norm of overlap integrals between vibrational wavefunctions of different electronic states. FC profiles are typically highly congested for large molecular systems and the spectra usually become not well-resolvable at elevated temperatures. The (theoretical) analyses of such spectra are even more difficult when vibrational mode mixing (Duschinsky) effects are significant, because contributions from different modes are in general not separable, even within the harmonic approximation. A few decades ago Doktorov, Malkin and Man'ko [1979 J. Mol. Spectrosc. 77, 178] developed a coherent state-based generating function approach and exploited the dynamical symmetry of vibrational Hamiltonians for the Duschinsky relation to describe FC transitions at zero Kelvin. Recently, the present authors extended the method to incorporate thermal, single vibronic level, non-Condon and multi-photon effects in energy, time and probability density domains for the efficient calculation and interpretation of vibronic spectra. Herein, recent developments and corresponding generating functions are presented for single vibronic levels related to fluorescence, resonance Raman scattering and anharmonic transition.
Biological as well as advanced artificial intelligences (AIs) need to decide which goals to pursue. We review nature's solution to the time allocation problem, which is based on a continuously readjusted categorical weighting mechanism we experience introspectively as emotions. One observes phylogenetically that the available number of emotional states increases hand in hand with the cognitive capabilities of animals and that raising levels of intelligence entail ever larger sets of behavioral options. Our ability to experience a multitude of potentially conflicting feelings is in this view not a leftover of a more primitive heritage, but a generic mechanism for attributing values to behavioral options that can not be specified at birth. In this view, emotions are essential for understanding the mind. For concreteness, we propose and discuss a framework which mimics emotions on a functional level. Based on time allocation via emotional stationarity (TAES), emotions are implemented as abstract criteria, such as satisfaction, challenge and boredom, which serve to evaluate activities that have been carried out. The resulting timeline of experienced emotions is compared with the “character” of the agent, which is defined in terms of a preferred distribution of emotional states. The long-term goal of the agent, to align experience with character, is achieved by optimizing the frequency for selecting individual tasks. Upon optimization, the statistics of emotion experience becomes stationary.
We analyzed a eukaryotically encoded rubredoxin from the cryptomonad Guillardia theta and identified additional domains at the N- and C-termini in comparison to known prokaryotic paralogous molecules. The cryptophytic N-terminal extension was shown to be a transit peptide for intracellular targeting of the protein to the plastid, whereas a C-terminal domain represents a membrane anchor. Rubredoxin was identified in all tested phototrophic eukaryotes. Presumably facilitated by its C-terminal extension, nucleomorph-encoded rubredoxin (nmRub) is associated with the thylakoid membrane. Association with photosystem II (PSII) was demonstrated by co-localization of nmRub and PSII membrane particles and PSII core complexes and confirmed by comparative electron paramagnetic resonance measurements. The midpoint potential of nmRub was determined as +125 mV, which is the highest redox potential of all known rubredoxins. Therefore, nmRub provides a striking example of the ability of the protein environment to tune the redox potentials of metal sites, allowing for evolutionary adaption in specific electron transport systems, as for example that coupled to the PSII pathway.
The core of this work is represented by the investigation of the chiral phase transition, using Monte Carlo simulations and unimproved staggered fermions, both in the weak and strong coupling regimes of Quantum Chromodynamics. Based on recent results from Monte Carlo simulations, both using unimproved staggered fermions and Wilson fermions, the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf (number of flavours) in range [2:7], at zero baryon chemical potential. This achievement relies on the analytic continuation of Nf to non-integer values on the lattice, which allows to make use of extrapolation techniques to the chiral limit, where simulations are not possible. Furthermore, these results provide a resolution to the ambiguous scenario for Nf = 2 in the chiral limt. The first part of this thesis is devoted to the investigation of the chiral phase transition when a non-zero imaginary baryon chemical potential is involved, whose value corresponds to the 81% of the Roberge-Weiss one. Using the same extrapolation techniques aforementioned, the order of the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf in range [2:6], highlighting a lack of dependence of the order of the chiral phase transition on the imaginary baryon chemical potential value. The second part of this thesis is about the study of the extension of the first-order chiral region in the strong coupling regime, at zero baryon chemical potential. Using Monte Carlo techniques, this can be done by investigating the Z2 boundary on a coarse lattice, whose temporal extent reads Nt = 2, and simulations are realised for Nf = 4, 8. The results in the weak coupling regime show, for $Nt = 8, 6, 4 and fixed Nf value, an inflating first-order chiral region. As in the strong coupling limit a second-order chiral phase transition is expected, the first-order chiral region has to shrink as the strong coupling regime is approached, resulting in a non-monotonic behaviour of the Z2 boundary. For Nf = 8, a critical mass on the Z2 boundary has been obtained, confirming the expected non-monotonic behaviour. For Nf = 4 the results do not provide a unique conclusion: Either a Z2 boundary at extremely low bare quark mass or a second-order chiral phase transition in the O(2) universality class in the chiral limit can take place. In addition to the two main topics, the performances of the second-order minimum norm integrator (2MN) and the fourth-order minimum norm integrator (4MN) have been compared, after implementing the 4MN one in the CL2QCD code used to realise our simulations. The 2MN integrator had already been implemented in the code since the first version was released. The two integrators belong to the class of symplectic integrators and represent an essential component of the RHMC algorithm, involved in our investigation. This step is extremely important, in order to guarantee the best quality when collecting data from simulations, and the results of the comparison suggested to favor the 2MN integrator, for both the topics.
A powerful technique to distinguish the enantiomers of a chiral molecule is the Coulomb Explosion Imaging (CEI). This technique allows us to determine the handedness of a single molecule. In CEI, the molecule becomes charged by losing many electrons in a very short period of time by interacting with the light. The repulsion forces between the positive charged particles of the molecule leads the molecule to break into parts-fragments. By measuring the three vector momentum of (at least) four fragments, the handedness observable can be determined. In this thesis, CEI is induced by absorption of a single high energy photon, which creates an inner-shell hole (K shell) of the molecule. The subsequent cascade of Auger decays lead to fragmentation. We decided to work with the formic acid molecule in this thesis. Two different experiments were conducted. The first experiment focused on exciting electrons to different energy states, while the second experiment focused on extracting directly a photoelectron to the continuum and measure the angular distribution of the photoelectron in the molecular frame. The primary goal was to search for chiral signal in a pure achiral planar molecule under the previous electron processes. The results of these findings were further implemented to two more molecules.
In the framework of the LHC Injectors Upgrade Project (LIU), the CERN Proton Synchrotron Booster (PSB) went through major upgrades resulting in new effects to study, challenges to overcome and new parameter regimes to explore. To assess the achievable beam brightness limit of the machine, a series of experimental and computational studies in the transverse planes were performed. In particular, the new injection scheme induces optics perturbations that are strongly enhanced near the half-integer resonance. In this thesis, methods for dynamically measuring and correcting these perturbations and their impact on the beam performance will be presented. Additionally, the quality of the transverse beam distributions and strategies for improvement will be addressed. Finally, the space charge effects when dynamically crossing the half-integer resonance will be characterized. The results of these studies and their broader significance beyond the PSB will be discussed.
IHMCIF: an extension of the PDBx/mmCIF data standard for integrative structure determination methods
(2024)
IHMCIF (github.com/ihmwg/IHMCIF) is a data information framework that supports archiving and disseminating macromolecular structures determined by integrative or hybrid modeling (IHM), and making them Findable, Accessible, Interoperable, and Reusable (FAIR). IHMCIF is an extension of the Protein Data Bank Exchange/macromolecular Crystallographic Information Framework (PDBx/mmCIF) that serves as the framework for the Protein Data Bank (PDB) to archive experimentally determined atomic structures of biological macromolecules and their complexes with one another and small molecule ligands (e.g., enzyme cofactors and drugs). IHMCIF serves as the foundational data standard for the PDB-Dev prototype system, developed for archiving and disseminating integrative structures. It utilizes a flexible data representation to describe integrative structures that span multiple spatiotemporal scales and structural states with definitions for restraints from a variety of experimental methods contributing to integrative structural biology. The IHMCIF extension was created with the benefit of considerable community input and recommendations gathered by the Worldwide Protein Data Bank (wwPDB) Task Force for Integrative or Hybrid Methods (wwpdb.org/task/hybrid). Herein, we describe the development of IHMCIF to support evolving methodologies and ongoing advancements in integrative structural biology. Ultimately, IHMCIF will facilitate the unification of PDB-Dev data and tools with the PDB archive so that integrative structures can be archived and disseminated through PDB.
The spike protein of SARS-CoV-2 is a highly flexible membrane receptor that triggers the translocation of the virus into cells by attaching to the human receptors. Like other type I membrane receptors, this protein has several extracellular domains connected by flexible hinges. The presence of these hinges results in high flexibility, which consequently results in challenges in defining the conformation of the protein. Here, We developed a new method to define the conformational space based on a few variables inspired by the robotic field’s methods to determine a robotic arm’s forward kinematics. Using newly performed atomistic molecular dynamics (MD) simulations and publicly available data, we found that the Denavit-Hartenberg (DH) parameters can reliably show the changes in the local conformation. Furthermore, the rotational and translational components of the homogenous transformation matrix constructed based on the DH parameters can identify the changes in the global conformation of the spike and also differentiate between the conformation with a similar position of the spike head, which other types of parameters, such as spherical coordinates, fail to distinguish between such conformations. Finally, the new method will be beneficial for looking at the conformational heterogeneity in all other type I membrane receptors.
Highlights
• Sampling the large conformational space of disordered proteins requires extensive molecular dynamics (MD) simulations.
• Fragment assembly complements MD simulations to produce extensive ensembles of disordered proteins with atomic detail.
• Hierarchical chain growth (HCG) ensembles capture key experimental descriptors “out of the box”.
• HCG has revealed local structural characteristics associated with protein dysfunction in neurodegeneration.
Abstract
Disordered proteins and nucleic acids play key roles in cellular function and disease. Here, we review recent advances in the computational exploration of the conformational dynamics of flexible biomolecules. While atomistic molecular dynamics (MD) simulation has seen a lot of improvement in recent years, large-scale computing resources and careful validation are required to simulate full-length disordered biopolymers in solution. As a computationally efficient alternative, hierarchical chain growth (HCG) combines pre-sampled chain fragments in a statistically reproducible manner into ensembles of full-length atomically detailed biomolecular structures. Experimental data can be integrated during and after chain assembly. Applications to the neurodegeneration-linked proteins α-synuclein, tau, and TDP-43, including as condensate, illustrate the use of HCG. We conclude by highlighting the emerging connections to AI-based structural modeling including AlphaFold2.
This paper reports on Monte Carlo simulation results for future measurements of the moduli of time-like proton electromagnetic form factors, |GE | and |GM|, using the ¯pp → μ+μ− reaction at PANDA (FAIR). The electromagnetic form factors are fundamental quantities parameterizing the electric and magnetic structure of hadrons. This work estimates the statistical and total accuracy with which the form factors can be measured at PANDA, using an analysis of simulated data within the PandaRoot software framework. The most crucial background channel is ¯pp → π+π−,due to the very similar behavior of muons and pions in the detector. The suppression factors are evaluated for this and all other relevant background channels at different values of antiproton beam momentum. The signal/background separation is based on a multivariate analysis, using the Boosted Decision Trees method. An expected background subtraction is included in this study, based on realistic angular distribuations of the background contribution. Systematic uncertainties are considered and the relative total uncertainties of the form factor measurements are presented.
How long does it take to emit an electron from an atom? This question has intrigued scientists for decades. As such emission times are in the attosecond regime, the advent of attosecond metrology using ultrashort and intense lasers has re-triggered strong interest on the topic from an experimental standpoint. Here, we present an approach to measure such emission delays, which does not require attosecond light pulses, and works without the presence of superimposed infrared laser fields. We instead extract the emission delay from the interference pattern generated as the emitted photoelectron is diffracted by the parent ion’s potential. Targeting core electrons in CO, we measured a 2d map of photoelectron emission delays in the molecular frame over a wide range of electron energies. The emission times depend drastically on the photoelectrons’ emission directions in the molecular frame and exhibit characteristic changes along the shape resonance of the molecule.
Upon antibiotic stress Gram-negative pathogens deploy resistance-nodulation-cell division-type tripartite efflux pumps. These include a H+/drug antiporter module that recognizes structurally diverse substances, including antibiotics. Here, we show the 3.5 Å structure of subunit AdeB from the Acinetobacter baumannii AdeABC efflux pump solved by single-particle cryo-electron microscopy. The AdeB trimer adopts mainly a resting state with all protomers in a conformation devoid of transport channels or antibiotic binding sites. However, 10% of the protomers adopt a state where three transport channels lead to the closed substrate (deep) binding pocket. A comparison between drug binding of AdeB and Escherichia coli AcrB is made via activity analysis of 20 AdeB variants, selected on basis of side chain interactions with antibiotics observed in the AcrB periplasmic domain X-ray co-structures with fusidic acid (2.3 Å), doxycycline (2.1 Å) and levofloxacin (2.7 Å). AdeABC, compared to AcrAB-TolC, confers higher resistance to E. coli towards polyaromatic compounds and lower resistance towards antibiotic compounds.
DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2
(2023)
More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.
Collective flow phenomena are a sensitive probe for the properties of extreme QCD matter. However, their interpretation relies on the understanding of the initial conditions e.g. the eccentricity of the nuclear overlap region. HADES [1] provides a large acceptance combined with a high mass-resolution and therefore allows to study di-electron and hadron production in heavy-ion collisions with unprecedented precision. In this contribution, the capability of HADES to study flow harmonics by utilizing multi-particle azimuthal correlation techniques is discussed. Due to the high statistics of seven billion Au+Au collisions at 1.23 AGeV collected in 2012, a systematic study of higher-order flow harmonics, the differentiation between collective and non-flow effects, and as well the multi-differential (pt, rapidity, centrality) analysis is possible.
In this thesis, the flow coefficients vn of the orders n = 1 − 6 are studied for protons and light nuclei in Au+Au collisions at Ebeam = 1.23 AGeV, equivalent to a center-of-mass energy in the nucleon-nucleon system of √sNN = 2.4 GeV. The detailed multi-differential measurement is performed with the HADES experiment at SIS18/GSI. HADES, with its large acceptance, covering almost full azimuth angle, combined with its high mass-resolution and good particle-identification capability, is well equipped to study the azimuthal flow pattern not only for protons, deuterons, and tritons but also for charged pions, kaons, the φ-mesons, electrons/positrons, as well as light nuclei like helions and alphas. The high statistics of more than seven billion Au-Au collisions recorded in April/May 2012 with HADES enables for the first time the measurement of higher order flow coefficients up to the 6th harmonic. Since the Fourier coefficient of 7th and 8th order are beyond the statistical significance only an upper bound is given. The Au+Au collision system is the largest reaction system with the highest particle multiplicities, which was measured so far with HADES. A dedicated correction method for the flow measurement had to be developed to cope with the reconstruction in-efficiencies due to occupancies of the detector system. The systematical bias of the flow measurement is studied and several sources of uncertainties identified, which mainly arise from the quality selection criteria applied to the analyzed tracks, the correction procedure for reconstruction inefficiencies, the procedures for particle identification (PID) and the effects of an azimuthally non-uniform detector acceptance. The systematic point-to-point uncertainties are determined separately for each particle type (proton, deuteron and triton), the order of the flow harmonics vn, and the centrality class. Further, the validity of the results is inspected in the range of their evaluated systematic uncertainties with several consistency checks. In order to enable meaningful comparisons between experimental observations and predictions of theoretical models, the classification of events should be well defined and in sufficiently narrow intervals of impact parameter. Part of this work included the implementation of the procedure to determine the centrality and orientation of the reaction.
In the conclusion the experimental results are discussed, including various scaling properties of the flow harmonics. It is found that the ratio v4/v2 for protons and light nuclei (deuterons and tritons) at midrapidity for all centrality classes approaches values close to 0.5 at high transverse momenta, which was suggested to be indicative for an ideal hydrodynamic behaviour. A remarkable scaling is observed in the pt dependence of v2 (v4) at mid-rapidity of the three hydrogen isotopes, when dividing by their nuclear mass number A (A^2) and pt by A. This is consistent with naive expectations from nucleon coalescence, butraises the question whether this mass ordering can also be explained by a hydrodynamical-inspired approach, like the blast-wave model. The relation of v2 and v4 to the shape of the initial eccentricity of the collision system is studied. It is found that v2 is independent of centrality for all three particle species after dividing it by the averaged second order participant eccentricity v2/⟨ε2⟩. A similar scaling is shown for v4 after division by ⟨ε2⟩^2.
Many Polyakov loop models can be written in a dual formulation which is free of sign problem even when a non-vanishing baryon chemical potential is introduced in the action. Here, results of numerical simulations of a dual representation of one such effective Polyakov loop model at finite baryon density are presented. We compute various local observables such as energy density, baryon density, quark condensate and describe in details the phase diagram of the model. The regions of the first order phase transition and the crossover, as well as the line of the second order phase transition, are established. We also compute several correlation functions of the Polyakov loops.
The effect of a single site mutation of Arg-54 to methionine in Paracoccus denitrificans cytochrome c oxidase was studied using a combination of optical spectroscopy, electrochemical and rapid kinetics techniques, and time-resolved measurements of electrical membrane potential. The mutation resulted in a blue-shift of the heme a alpha-band by 15 nm and partial occupation of the low-spin heme site by heme O. Additionally, there was a marked decrease in the midpoint potential of the low-spin heme, resulting in slow reduction of this heme species. A stopped-flow investigation of the reaction with ferrocytochrome c yielded a kinetic difference spectrum resembling that of heme a(3). This observation, and the absence of transient absorbance changes at the corresponding wavelength of the low-spin heme, suggests that, in the mutant enzyme, electron transfer from Cu(A) to the binuclear center may not occur via heme a but that instead direct electron transfer to the high-spin heme is the dominating process. This was supported by charge translocation measurements where Deltapsi generation was completely inhibited in the presence of KCN. Our results thus provide an example for how the interplay between protein and cofactors can modulate the functional properties of the enzyme complex.
We show examples of the impact of the Maxwellian averaged capture cross sections determined at n_TOF over the past 20 years on AGB stellar nucleosynthesis models. In particular, we developed an automated procedure to derive MACSs from evaluated data libraries, which are subsequently used as input to stellar models computed by means of the FuNS code. In this contribution, we present a number of s-process abundances obtained using different data libraries as input to stellar models, with a focus on the role of n_TOF data.
Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.
Present nuclear reaction network computations for astrophysical simulations involve many different types of rates, including neutron-capture reactions of interest for the modeling of heavy-element nucleosynthesis. While for many of them we still have to rely on theoretical calculations, an increasing number of experimentally-determined cross sections have now become available. In this contribution, we present “ASTrophysical Rate and rAw data Library” (ASTRAL), a new online database for neutron-capture cross sections based on experimental results, mainly obtained through activation and timeof-flight measurements. For the evaluation process, cross sections were re-calculated starting from raw data and by considering recent changes in physical properties of the involved isotopes (e.g., half-life and γ-ray intensities). We show the current status of the database, the techniques adopted to derive the recommended Maxwellian-averaged cross sections, and future developments.
Prediction for hyper nuclei multiplicities from GSI to LHC energies from the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model combined with a final state coalescence approach is presented and compared to the thermal model. The influence of the coalescence radius on the collision energy and centrality dependence of the Λ3H/ΛΛ3H/Λ ratio is discussed.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
This article summarizes some of the current theoretical developments and the experimental status of hypernuclei in relativistic heavy-ion collisions and elementary collisions. In particular, the most recent results of hyperhydrogen of mass A = 3 and 4 are discussed. The highlight at SQM2022 in this perspective was the discovery of the anti-hyperhydrogen-4 by the STAR Collaboration, in a large data set consisting of different collision systems. Furthermore, the production yields of hyperhydrogen-4 and hyperhelium-4 from the STAR Collaboration can be described nicely by the thermal model when the excited states of these hypernuclei are taken into account. In contrast, the production measurements in small systems (pp and p–Pb) from the ALICE Collaboration tends to favour the coalescence model over the thermal description. New measurements from STAR, ALICE and HADES Collaborations of the properties, e.g. lifetime, of A = 3 and 4 hypernuclei give similar results of these properties. Also the anti-hyperhydrogen-4 lifetime is in rather good agreement with previous measurements. Interestingly, the new STAR measurement on the R3 value, that is connected to the branching ratio, points to a Λ separation energy that is below 100 keV but definitely consistent with the value of 130 keV assumed since the 70s.
Asymptotic giant branch (AGB) stars are responsible for the production of the main component of the solar s-process distribution. Despite enormous progress in the theoretical modeling of these objects over the last few decades, many uncertainties remain. The still-unknown mechanism leading to the production of 13C neutron source is one example. The nucleosynthetic signature of AGB stars can be examined in a number of stellar sources, from spectroscopic observations of intrinsic and extrinsic stars to the heavy-element isotopic composition of presolar grains found in meteorites. The wealth of available observational data allows for constraining the processes occurring in AGB interiors. In this view, we discuss recent results from new AGB models including the effects of mixing triggered by magnetic fields, and show comparisons of the related s-process nucleosynthesis with available observations.
Lattice QCD and functional methods are making significant progress in constraining the QCD phase diagram. As an important milestone, the chiral phase transition with massless u, d-quarks at zero density is now understood to be of second order for all strange quark masses, and a smooth crossover as soon as mu,d, ≠ 0. Together with information on fluctuations and refined reweighted simulations, this bounds a possible critical point to be at µB/T ≲3. On the other hand, an approximately chiral-spin symmetric temperature window has been discovered above the chiral crossover, Tch<T ≳3Tch, with distinct correlator multiplet patterns and a pion spectral function suggesting resonance-like degrees of freedom, which dissolve graduallly with temperature.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
The phase diagram of the (1+1)-dimensional Gross-Neveu model is reanalyzed for (non-)zero chemical potential and (non-)zero temperature within the mean-field approximation. By investigating the momentum dependence of the bosonic two-point function, the well-known second-order phase transition from the Z2 symmetric phase to the so-called inhomogeneous phase is detected. In the latter phase the chiral condensate is periodically varying in space and translational invariance is broken. This work is a proof of concept study that confirms that it is possible to correctly localize second-order phase transition lines between phases without condensation and phases of spatially inhomogeneous condensation via a stability analysis of the homogeneous phase. To complement other works relying on this technique, the stability analysis is explained in detail and its limitations and successes are discussed in context of the Gross-Neveu model. Additionally, we present explicit results for the bosonic wave-function renormalization in the mean-field approximation, which is extracted analytically from the bosonic two-point function. We find regions -- a so-called moat regime -- where the wave function renormalization is negative accompanying the inhomogeneous phase as expected.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
We explore the phase structure of the 1+1 dimensional Gross-Neveu model at finite number of fermion flavors using lattice field theory. Besides a chirally symmetric phase and a homogeneously broken phase we find evidence for the existence of an inhomogeneous phase, where the condensate is a spatially oscillating function. Our numerical results include a crude μ-T phase diagram.
We explore the phase structure of the 1+1 dimensional Gross-Neveu model at finite number of fermion flavors using lattice field theory. Besides a chirally symmetric phase and a homogeneously broken phase we find evidence for the existence of an inhomogeneous phase, where the condensate is a spatially oscillating function. Our numerical results include a crude μ-T phase diagram.
In this work, the phase diagram of the 2+1-dimensional Gross-Neveu model is investigated with baryon chemical potential as well as chiral chemical potential in the mean-field approximation. We study the theory using two lattice discretizations, which are both based on naive fermions. An inhomogeneous chiral phase is observed only for one of the two discretizations. Our results suggest that this phase disappears in the continuum limit.
In this work, the phase diagram of the 2+1-dimensional Gross-Neveu model is investigated with baryon chemical potential as well as chiral chemical potential in the mean-field approximation. We study the theory using two lattice discretizations, which are both based on naive fermions. An inhomogeneous chiral phase is observed only for one of the two discretizations. Our results suggest that this phase disappears in the continuum limit.
In this work we study the 3+1-dimensional Nambu-Jona-Lasinio (NJL) model in the mean field-approximation. We carry out calculations using five different regularization schemes (two continuum and three lattice regularization schemes) with particular focus on inhomogeneous phases and condensates. The regularization schemes lead to drastically different inhomogeneous regions. We provide evidence that inhomogeneous condensates appear for all regularization schemes almost exclusively at values of the chemical potential and with wave numbers, which are of the order of or even larger than the corresponding regulators. This can be interpreted as indication that inhomogeneous phases in the 3+1-dimensional NJL model are rather artifacts of the regularization and not a consequence of the NJL Lagrangian and its symmetries.
In this work, inhomogeneous chiral phases are studied in a variety of Four-Fermion and Yukawa models in 2+1 dimensions at zero and non-zero temperature and chemical potentials. Employing the mean-field approximation, we do not find indications for an inhomogeneous phase in any of the studied models. We show that the homogeneous phases are stable against inhomogeneous perturbations. At zero temperature, full analytic results are presented.