Refine
Year of publication
Document Type
- Article (1891)
- Preprint (1261)
- Doctoral Thesis (594)
- Conference Proceeding (248)
- diplomthesis (101)
- Bachelor Thesis (75)
- Master's Thesis (61)
- Contribution to a Periodical (46)
- Book (33)
- Diploma Thesis (33)
Keywords
- Kollisionen schwerer Ionen (47)
- heavy ion collisions (44)
- LHC (25)
- Quark-Gluon-Plasma (25)
- Heavy Ion Experiments (20)
- equation of state (19)
- quark-gluon plasma (19)
- Relativistic heavy-ion collisions (16)
- BESIII (15)
- QCD (15)
Institute
- Physik (4438) (remove)
Neutron stars are unique laboratories for the investigation of the high density properties of bulk matter. In this work, the astrophysical constraints for a phase transition from hadronic matter to deconfined quark matter are examined thoroughly. A scheme for relating known astrophysical observables such as mass, radius and tidal deformability to the parameter space of such a transition is devised and applied to the set of data currently available.
In order to span a wide parameter space, a highly parameterizable relativistic mean field equation in compliance with chiral effective field theory results is used, where the stiffness of the equation of state can be varied via the effective mass at saturation density. The phase transitions are modelled using a Maxwell construction and assumed to be of first order, with a constant speed of sound quark matter model. The resulting equations of state are analyzed and divided into four categories, which can be used to constrain the parameter space that allows phase transition. It is highlighted, that a subset of this parameter space would even be detectable without the need of higher precision measurements. A phase transition at high densities is shown to be particularly promising in this regard. Finally, the groundwork is laid to apply the equation of state used in this work for supernova or merger simulations, by extending it to non-zero temperatures.
Die Wirkungsweise des Protonen-Wendellinearbeschleunigers und sein prinzipieller Aufbau werden kurz beschrieben und Angaben über die erreichbare Parallelimpedanz gemacht. Es wird gezeigt, wie durch sektionsweisen Aufbau eine variable Endenergie erzielt und die Kühl- und Fokussierprobleme gelöst werden können. Zwei Konstruktionsbeispiele für kontinuierlichen und gepulsten Betrieb werden näher ausgeführt.
Light-particle accompanied fission is expected to yield results from which one hopes to learn more about binary scission configurations. As a step in this direction, we present a model that allows the calculation of the probabilities with which a given three-particle setup follows from different binary configurations. First results show the workability of the model.
Über das Absorptionsverhalten verschiedener Quarzglassorten im Spektralbereich von 150 μ bis 275 μ
(1961)
Die Fermi-Verteilung gibt die Verteilung von Elektronen über die verschiedenen Energieniveaus in einem Festkörper unter der Voraussetzung an, daß sich ein thermodynamischer Gleichgewichtszustand eingestellt hat. Bei der Untersuchung von Leitfähigkeitsproblemen werden geringfügige Störungen dieser Gleichgewichtsverteilung durch elektrische Felder, Temperaturgefälle u. dgl. berücksichtigt. Der innere Photoeffekt bewirkt eine kräftige Störung der Gleichgewichtsverteilung, da fortgesetzt lichtelektrisch ausgelöste Elektronen aus einem niederen in einen höheren Energiezustand gehoben werden.
Im folgenden wird ein statistischer Ansatz für die Verteilung von Elektronen über dargebotene Energiestufen versucht für den Fall, daß laufend je Sekunde eine bestimmte Anzahl aus einem tieferen in ein höheres Niveau gehoben wird.
An alternative formulation is presented of the formal theory of multi-channel scattering in nonrelativistic quantum mechanics. We start by defining spaces of state vectors, where two particles either stay together or separate in the limit t →+∞ (or — ∞), when the state vector develops in time by e–i H t (H is the complete Hamiltonian of the n-particle system). A channel is defined as a space of state vectors with the following property: Developing in time by e-i H t they asymptotically describe a state of the n-particle system, where the particles are grouped in fragments. Defining a Hamiltonian Hγ for each channel, in which—compared to H—the interactions acting between particles from different fragments are missing, it is physically plausible that lim eiH e—iHt Ψ exists for vectors Ψ in the channel. Having discussed the limit vectors (asymptotic states), the S-matrix formalism can be introduced as usual. Finally the introduction of the exclusion principle is discussed.
Es wird eine massenspektrometrische Analysenmethode für Germanium angegeben, die den quantitativen Nachweis einiger Elemente in Germanium bis zu Relativkonzentrationen von 10-5—10-6 gestattet. Diese Methode wurde auf die Messung der Verteilungskoeffizienten von B und In in Ge angewandt. Die Analysen wurden an dem mit Hilfe der Zonenschmelze1 angereicherten Ende eines Regulus durchgeführt.
Die massenspektrometrischen Analysen des Germaniumdampfes zeigten einen molekularen Ge-Anteil in der Dampfphase. Es wurden Cluster-Ionen bis zu Ge8 nachgewiesen. Auch bei Zinn und Silicium fand sich ein Cluster-Anteil in der Dampfphase.
Nach einer kurzen Zusammenstellung von Verfahren, die die Herstellung von Thermoelementen mit geringer Wärmekapazität und -ableitung ermöglichen, werden Konstruktion und thermische Daten der vom Verfasser entwickelten mikroskopisch feinen „Thermonadeln“ wiedergegeben. Um die meßtechnischen Möglichkeiten dieser Thermonadeln beurteilen zu können, wird zunächst eine einfache Theorie des Thermoelementes als Wärmekraftmaschine entwickelt: Der Nutzeffekt wird definiert als der Quotient aus der durch Peltier- und Thomson-Effekt erzeugten elektrischen Energie und den durch die an der Meßstelle des Thermoelementes reversibel und irreversibel aufgenommenen Wärmemengen. Mit Hilfe eines „Arbeitsfaktors A“ (der ein Maß für die Qualität des Thermopaares ist) wird ein theoretischer und praktischer „Gütefaktor G bzw. G“ als das Verhältnis des Nutzeffektes eines Thermoelementes zum Nutzeffekt einer Carnot-Maschine definiert. Dieses Verhältnis gibt an, um welchen Faktor der Nutzeffekt eines Thermoelementes „theoretisch“ und „praktisch“ hinter dem thermodynamisch größtmöglichen Nutzeffekt zurückbleibt. Aus der Tatsache, daß G niemals größer als 1 werden kann, wird eine obere Grenze für die Arbeitsfaktoren A ≦ 1 und damit (in Abhängigkeit.
Die Vorstellungen des Wilsonschen Halbleitermodells werden quantitativ durchgeführt und die Abhängigkeit der Fermischen Grenzenergie von Temperatur, Konzentration und Art der Störstellen berechnet. Es ergibt sich, daß die Grenzenergie das elektrische Verhalten des Halbleiters entscheidet. Bisher ungeklärte experimentelle Ergebnisse werden gedeutet.
Nach neueren Messungen besitzt das Auge eine definierte Empfindlichkeit noch bis etwa 1.0 μ. Es wird hier eine Formel mitgeteilt, die die Messungsresultate befriedigend darstellt. Ihre Auslegung deutet darauf hin, daß der Abfall der Augenempfindlichkeit in diesem "ultraroten" Bereich gerade hinreichend rasch erfolgt, daß keine Dunkelreaktion im Auge den Sehakt stört.
Eine Erweiterung des Differenzen-Verfahrens zur Berechnung von kugelsymmetrischen Temperaturfeldern
(1947)
Mit dem Differenzenverfahren können in einfacher Weise beliebige Anheiz- und Abkühlungsvorgänge, soweit sie kugelsymmetrisch verlaufen, behandelt werden. Sowohl der Fall, daß eine Kugel von einem Medium konstanter Außentemperatur umgeben wird, wie auch der Fall, daß die Temperatur ohne Störung durch Konvektion in der Umgebung der geheizten Kugel sich einstellt, wird erörtert und führt auf leicht durchführbare Rechnungen und Konstruktionen. Zum Schluß wird gezeigt, wie sich aus den angegebenen Rekursionsformeln die analytischen Gesetzmäßigkeiten der stationären Temperaturverteilung leicht elementar ableiten lassen.
Based on an e+e− collision data sample corresponding to an integrated luminosity of 2.93 fb−1 collected with the BESIII detector at √s=3.773 GeV, the first amplitude analysis of the singly Cabibbo-suppressed decay D+→K+K0Sπ0 is performed. From the amplitude analysis, the K∗(892)+K0S component is found to be dominant with a fraction of (57.1±2.6±4.2)%, where the first uncertainty is statistical and the second systematic. In combination with the absolute branching fraction B(D+→K+K0Sπ0) measured by BESIII, we obtain B(D+→K∗(892)+K0S)=(8.69±0.40±0.64±0.51)×10−3, where the third uncertainty is due to the branching fraction B(D+→K+K0Sπ0). The precision of this result is significantly improved compared to the previous measurement. This result also differs from most of theoretical predictions by about 4σ, which may help to improve the understanding of the dynamics behind.
In Memoriam: Kurt Kimpel †
(1950)
An Hand neuerer paläoklimatologischer Forschungsergebnisse wird nachgewiesen, daß sowohl die Polwanderungstheorie als auch die astronomische Theorie von Milankovitch in bezug auf die Verhältnisse während des Känozoikums und namentlich während des Alluviums versagen. Es wird weiterhin ausgeführt, daß die Nebelveränderlichkeit nicht etwa durch rein optische Bedeckungseffekte, sondern durch die Bildung einer dichten, den Stern umgebenden Hülle, welche den Energietransport in der Sternatmosphäre stört, hervorgerufen wird. Die kleineren Schwankungen innerhalb der Eiszeiten (Interstadialzeiten) und die kurzdauernden Schwankungen des Alluvialklimas können durch die faserige (Filament-) Struktur der Dunkelwolken erklärt werden.
Im ersten Teil werden kurz die entscheidenden Schwierigkeiten der beiden bekanntesten Versuche zur Erklärung der Eiszeiten, und zwar der Pol-bzw. Kontinentalver-schiebungen (Epeirophorese) und der sog. astronomischen Theorie (Schwankungen der Schiefe der Ekliptik und der Erdbahnelemente) aufgezeigt. - Im zweiten Teil wird wahrscheinlich gemacht, daß eine schon 1921 vorgetragene Hypothese von Shapley, nach der die Sonne in der diluvialen Eiszeit eine von kosmischen Dunkelwolken verursachte schwache Veränderlichkeit ihrer Strahlung zeigte, nach neueren, in erster Linie astrophysikalischen Forschungsergebnissen eine brauchbare Erklärung zu liefern imstande ist, aus folgenden drei Gründen:
1. In allen dichteren kosmischen Dunkelwolken findet man unregelmäßige veränderliche Sterne eines ganz besonderen Typus; im Gegensatz zu fast allen übrigen Typen veränderlicher Sterne handelt es sich bei diesen Sternen um Zwergsterne, wie die Sonne oder noch schwächer.
2. Unser Sonnensystem befindet sich gegenwärtig innerhalb dunkler kosmischer Materie, und es ist höchstwahrscheinlich, daß es in naher Vergangenheit auch dichtere Teile von solchen Dunkelwolken durchquert hat.
3. Die Form der Dunkelwolken, ihre armförmigen, zirrus-oder zirrostratusähnlichen Anordnungen und Verästelungen lassen eine zwanglose Erklärung der mehrfachen Wiederholung der Eiszeiten zu; es wird hier also nicht nur die Eiszeit, sondern aus der gleichen Hauptursache auch deren Gliederung erklärt, was sonst nicht gelungen zu sein scheint.
Die Mitwirkung anderer, in erster Linie astronomischer und geographischer Faktoren, letztere insbesondere in der älteren Erdgeschichte, ist durchaus möglich.
During collisions of heavy ions with heavy targets below the Coulomb barrier, adiabatic molecular orbitals are formed for the inner electrons. Deviations from adiabaticity lead to coupling between various states and can be treated by time-dependent perturbation theory. For high charges ( Z1+Z2 ≧ 60) the molecular electrons are highly relativistic. Therefore, the Dirac equation has to be used to obtain the energies and wave functions. The Dirac Hamiltonian is transformed into the intrinsic rotating coordinate system where prolate spheroidal coordinates are introduced. A set of basis functions is proposed which allows the evaluation of all matrix elements of the Dirac Hamiltonian analytically. The resulting matrix is diagonalized numerically. The finite nuclear charge distribution is also taken into account. Results are presented and discussed for various characteristic systems, e. g. Br-Br, Ni-Ni, I-I, Br-Zr, I-Au, U -U, etc.
We present a unified formulation of the interaction of electrons with the electromagnetic field in heavy ion collisions, based on quantized interacting fields. This reduces the effort in treating many-electron systems substantially, as compared with the usual S-matrix theory. Both formalisms are shown to be equivalent. The simplification achieved by our new approach is demonstrated in detail for the example of quasi-molecular radiation.
A model for the quantum yield of the coloration caused by UV-light in spiropyran layers is described. This model allows to calculate the sensitivity of layers having different compositions. The mechanism concerning the stability of the coloration is essentially clarified. Calculations of the stability for layers of different compositions are possible by a model describing the mechanism approximately.
A scattering theory for reactions with three-particle channels above the two-particle threshold is developed. The S-matrix-technique is used for the calculation of the extended S-matrix. Correlated two-particle wave functions in the exit channels are employed to describe the exact two-particle continuum. For the usual shell model only a few partial waves dominate. The cross section depends on the energy-distribution between the two outgoing nucleons. Numerical results are presented for the model (d, 2n)-reaction exciting 0+ -states in O16 without Coulomb-effects. The treatment is restricted to three (2p2h)-states with the particles in the (sd)-shell and holes in the p-shell.
Two equations for the macroscopic part W of the statistical operator are considered:
1. the master equation W = — MW, t
2. the exact equation W = — J K(t — r) W (r) dr.
It follows from the physical equivalence of the solutions together with a stability assumption and the assumption that there is a time τ* after which also the derivatives of the solutions are equivalent, that τ* is the life-time of the kernel K and that Conversely, the equivalence of the solutions follows from assumptions on the life-time of the kernel K together with a stability assumption and a smoothness assumption on the initial statistical operator W(0).
Coupled channel calculations for Coulomb and nuclear excitation of the systems 136Xe-238U and 238U-238U have been performed using the rotation-vibration model. The impact parameter-, energy-and spin-dependence of the excitation probabilities are discussed for the ground state-, β-and γ-band up to Jπ = 36+. It is shown that the energy levels and quadrupole matrix elements are strongly influenced by the rotation-vibration interaction. Analytic expressions for the elastic and coupling potentials are presented.
Introducing correlated continuum wave functions for the two- and re-particle-continuum a microscopic theory of nuclear reactions based on a method of Fano is developed. The S-matrix-elements are given by the matrix-elements between correlated continuum wave functions and bound state wave functions. The antisymmetrization of the continuum wave functions with more than one particle in the continuum is included. The theory can be straightforwardly applied on the n-nucleon-emission process following photo- and particle excitations.
The master operators B which cause the entropy production dH/dt = - k-1 dS/dt to become extremal for fixed statistical operators W are constructed and discussed. There are boundaries of the set B of master operators, B = {B | Σ B2vu = b} for which the problem is solvable yielding minimal entropy production, while no solution exists in the set B without any constraints. Operators with maximal entropy production must be extremal points of B.
The spike protein (S) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is required for cell entry and is the primary focus for vaccine development. In this study, we combined cryo–electron tomography, subtomogram averaging, and molecular dynamics simulations to structurally analyze S in situ. Compared with the recombinant S, the viral S was more heavily glycosylated and occurred mostly in the closed prefusion conformation. We show that the stalk domain of S contains three hinges, giving the head unexpected orientational freedom. We propose that the hinges allow S to scan the host cell surface, shielded from antibodies by an extensive glycan coat. The structure of native S contributes to our understanding of SARS-CoV-2 infection and potentially to the development of safe vaccines.
The spike (S) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is required for cell entry and is the major focus for vaccine development. We combine cryo electron tomography, subtomogram averaging and molecular dynamics simulations to structurally analyze S in situ. Compared to recombinant S, the viral S is more heavily glycosylated and occurs predominantly in a closed pre-fusion conformation. We show that the stalk domain of S contains three hinges that give the globular domain unexpected orientational freedom. We propose that the hinges allow S to scan the host cell surface, shielded from antibodies by an extensive glycan coat. The structure of native S contributes to our understanding of SARS-CoV-2 infection and the development of safe vaccines. The large scale tomography data set of SARS-CoV-2 used for this study is therefore sufficient to resolve structural features to below 5 Ångstrom, and is publicly available at EMPIAR-10453.
Transport of lipids across membranes is fundamental for diverse biological pathways in cells. Multiple ion-coupled transporters take part in lipid translocation, but their mechanisms remain largely unknown. Major facilitator superfamily (MFS) lipid transporters play central roles in cell wall synthesis, brain development and function, lipids recycling, and cell signaling. Recent structures of MFS lipid transporters revealed overlapping architectural features pointing towards a common mechanism. Here we used cysteine disulfide trapping, molecular dynamics simulations, mutagenesis analysis, and transport assays in vitro and in vivo, to investigate the mechanism of LtaA, a proton-dependent MFS lipid transporter essential for lipoteichoic acid synthesis in the pathogen Staphylococcus aureus. We reveal that LtaA displays asymmetric lateral openings with distinct functional relevance and that cycling through outward- and inward-facing conformations is essential for transport activity. We demonstrate that while the entire amphipathic central cavity of LtaA contributes to lipid binding, its hydrophilic pocket dictates substrate specificity. We propose that LtaA catalyzes lipid translocation by a ‘trap-and-flip’ mechanism that might be shared among MFS lipid transporters.
A key event in cellular physiology is the decision between membrane biogenesis and fat storage. Phosphatidic acid (PA) is an important intermediate at the branch point of these pathways and is continuously monitored by the transcriptional repressor Opi1 to orchestrate lipid metabolism. In this study, we report on the mechanism of membrane recognition by Opi1 and identify an amphipathic helix (AH) for selective binding of PA over phosphatidylserine (PS). The insertion of the AH into the membrane core renders Opi1 sensitive to the lipid acyl chain composition and provides a means to adjust membrane biogenesis. By rational design of the AH, we tune the membrane-binding properties of Opi1 and control its responsiveness in vivo. Using extensive molecular dynamics simulations, we identify two PA-selective three-finger grips that tightly bind the PA phosphate headgroup while interacting less intimately with PS. This work establishes lipid headgroup selectivity as a new feature in the family of AH-containing membrane property sensors.
A key event in cellular physiology is the decision between membrane biogenesis and fat storage. Phosphatidic acid (PA) is an important lipid intermediate and signaling lipid at the branch point of these pathways and constantly monitored by the transcriptional repressor Opi1 to orchestrate lipid metabolism. Here, we report on the mechanism of membrane recognition by Opi1 and identify an amphipathic helix (AH) for the selective binding to membranes containing PA over phosphatidylserine (PS). The insertion of the AH into the hydrophobic core of the membrane renders Opi1 sensitive to the lipid acyl chain composition as an important factor contributing to the regulation of membrane biogenesis. Based on these findings, we rationally designed the membrane binding properties of Opi1 to control its responsiveness in the physiological context. Using extensive molecular dynamics (MD) simulations, we identified two PA-selective three-finger grips that tightly bind the phosphate headgroup, while interacting less intimately and more transiently with PS. This work establishes lipid headgroup selectivity as a new feature in the family of AH-containing membrane property sensors.
The SLC26 family of transporters maintains anion equilibria in all kingdoms of life. The family shares a 7 + 7 transmembrane segments inverted repeat architecture with the SLC4 and SLC23 families, but holds a regulatory STAS domain in addition. While the only experimental SLC26 structure is monomeric, SLC26 proteins form structural and functional dimers in the lipid membrane. Here we resolve the structure of an SLC26 dimer embedded in a lipid membrane and characterize its functional relevance by combining PELDOR/DEER distance measurements and biochemical studies with MD simulations and spin-label ensemble refinement. Our structural model reveals a unique interface different from the SLC4 and SLC23 families. The functionally relevant STAS domain is no prerequisite for dimerization. Characterization of heterodimers indicates that protomers in the dimer functionally interact. The combined structural and functional data define the framework for a mechanistic understanding of functional cooperativity in SLC26 dimers.
The SLC26 family of transporters maintains anion equilibria in all kingdoms of life. The family shares a 7 + 7 transmembrane segments inverted repeat architecture with the SLC4 and SLC23 families, but holds a regulatory STAS domain in addition. While the only experimental SLC26 structure is monomeric, SLC26 proteins form structural and functional dimers in the lipid membrane. Here we resolve the structure of an SLC26 dimer embedded in a lipid membrane and characterize its functional relevance by combining PELDOR distance measurements and biochemical studies with MD simulations and spin-label ensemble refinement. Our structural model reveals a unique interface different from the SLC4 and SLC23 families. The functionally relevant STAS domain exerts a stabilizing effect on regions central in this dimer. Characterization of heterodimers indicates that protomers in the dimer functionally interact. The combined structural and functional data define the framework for a mechanistic understanding of functional cooperativity in SLC26 dimers.
Molecular mechanisms of inorganic-phosphate release from the core and barbed end of actin filaments
(2023)
The release of inorganic phosphate (Pi) from actin filaments constitutes a key step in their regulated turnover, which is fundamental to many cellular functions. However, the molecular mechanisms underlying Pi release from both the core and barbed end of actin filaments remain unclear. Here, we combine cryo-EM with molecular dynamics simulations and in vitro reconstitution to demonstrate how actin releases Pi through a ‘molecular backdoor’. While constantly open at the barbed end, the backdoor is predominantly closed in filament-core subunits and only opens transiently through concerted backbone movements and rotameric rearrangements of residues close to the nucleotide binding pocket. This mechanism explains why Pi escapes rapidly from the filament end and yet slowly from internal actin subunits. In an actin variant associated with nemaline myopathy, the backdoor is predominantly open in filament-core subunits, resulting in greatly accelerated Pi release after polymerization and filaments with drastically shortened ADP-Pi caps. This demonstrates that the Pi release rate from F-actin is controlled by steric hindrance through the backdoor rather than by the disruption of the ionic bond between Pi and Mg2+ at the nucleotide-binding site. Our results provide the molecular basis for Pi release from actin and exemplify how a single, disease-linked point mutation distorts the nucleotide state distribution and atomic structure of the actin filament.
Maximum likelihood estimates of diffusion coefficients from single-particle tracking experiments
(2021)
Single-molecule localization microscopy allows practitioners to locate and track labeled molecules in biological systems. When extracting diffusion coefficients from the resulting trajectories, it is common practice to perform a linear fit on mean-squared-displacement curves. However, this strategy is suboptimal and prone to errors. Recently, it was shown that the increments between the observed positions provide a good estimate for the diffusion coefficient, and their statistics are well-suited for likelihood-based analysis methods. Here, we revisit the problem of extracting diffusion coefficients from single-particle tracking experiments subject to static noise and dynamic motion blur using the principle of maximum likelihood. Taking advantage of an efficient real-space formulation, we extend the model to mixtures of subpopulations differing in their diffusion coefficients, which we estimate with the help of the expectation–maximization algorithm. This formulation naturally leads to a probabilistic assignment of trajectories to subpopulations. We employ the theory to analyze experimental tracking data that cannot be explained with a single diffusion coefficient. We test how well a dataset conforms to the assumptions of a diffusion model and determine the optimal number of subpopulations with the help of a quality factor of known analytical distribution. To facilitate use by practitioners, we provide a fast open-source implementation of the theory for the efficient analysis of multiple trajectories in arbitrary dimensions simultaneously.
The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of the interaction remain elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarse-grained and atomistic models with coevolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted in a novel structural model that rationalizes previous experimental observations. We identify an evolutionarily conserved interaction surface formed by helix II of the DnaJ J-domain and a structurally contiguous region of DnaK, involving lobe IIA of the nucleotide binding domain, the inter-domain linker, and the β-basket of the substrate binding domain.
The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of this fundamental interaction are still elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarsegrained and atomistic models with co-evolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted into a novel structural model that rationalizes previous experimental observations. We identify an evolutionary-conserved interaction surface formed by helix II of the DnaJ J-domain and a groove on lobe IIA of the DnaK nucleotide binding domain, involving the inter-domain linker.
TriMem: A parallelized hybrid Monte Carlo software for efficient simulations of lipid membranes
(2022)
Lipid membranes are integral building blocks of living cells and perform a multitude of biological functions. Currently, molecular simulations of cellular-scale membrane remodeling processes at atomic resolution are extremely difficult, due to their size, complexity, and the large times-scales on which these processes occur. Instead, elastic membrane models are used to simulate membrane shapes and transitions between them and to infer their properties and functions. Unfortunately, an efficiently parallelized open-source simulation code to do so has been lacking. Here, we present TriMem, a parallel hybrid Monte Carlo simulation engine for triangulated lipid membranes. The kernels are efficiently coded in C++ and wrapped with Python for ease-of-use. The parallel implementation of the energy and gradient calculations and of Monte Carlo flip moves of edges in the triangulated membrane enable us to simulate large and highly curved membrane structures. For validation, we reproduce phase diagrams of vesicles with varying surface-to-volume ratios and area difference. We also compute the density of states to verify correct Boltzmann sampling. The software can be used to tackle a range of large-scale membrane remodeling processes as a step toward cell-scale simulations. Additionally, extensive documentation make the software accessible to the broad biophysics and computational cell biology communities.
Abstract
The primary immunological target of COVID-19 vaccines is the SARS-CoV-2 spike (S) protein. S is exposed on the viral surface and mediates viral entry into the host cell. To identify possible antibody binding sites, we performed multi-microsecond molecular dynamics simulations of a 4.1 million atom system containing a patch of viral membrane with four full-length, fully glycosylated and palmitoylated S proteins. By mapping steric accessibility, structural rigidity, sequence conservation, and generic antibody binding signatures, we recover known epitopes on S and reveal promising epitope candidates for structure-based vaccine design. We find that the extensive and inherently flexible glycan coat shields a surface area larger than expected from static structures, highlighting the importance of structural dynamics. The protective glycan shield and the high flexibility of its hinges give the stalk overall low epitope scores. Our computational epitope-mapping procedure is general and should thus prove useful for other viral envelope proteins whose structures have been characterized.
Author summary
The SARS-CoV-2 virus has caused a global health crisis. The spike protein exposed at its surface is key for infection and the primary antibody target. However, spike is covered by highly mobile glycan molecules that could impair antibody binding. To identify accessible epitopes, we performed molecular dynamics simulations of an atomistic model of glycosylated spike embedded in a membrane. By combining extensive simulations with bioinformatics analyses, we recovered known antibody binding sites and identified several epitope candidates as targets for further vaccine development.
More than 75% of surface and secreted proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, influencing protein dynamics and interactions with other molecules. Glycan conformational freedom impairs complete structural elucidation of glycoproteins. Computer simulations may be used to model glycan structure and dynamics. However, such simulations typically require thousands of computing hours on specialized supercomputers, thus limiting routine use. Here, we describe a reductionist method that can be implemented on personal computers to graft ensembles of realistic glycan conformers onto static protein structures in a matter of minutes. Using this open-source pipeline, we reconstructed the full glycan cover of SARS-CoV-2 Spike protein (S-protein) and a human GABAA receptor. Focusing on S-protein, we show that GlycoSHIELD recapitulates key features of extended simulations of the glycosylated protein, including epitope masking, and provides new mechanistic insights on N-glycan impact on protein structural dynamics.
Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ‘t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data.
Seit hundert Jahren ist bekannt, dass die mikroskopische Welt der Atome und Moleküle von den Gesetzen der Quantenphysik regiert wird. Lange Zeit galten Quantenphänomene als verworren und unkontrollierbar. Heute arbeiten Physikerinnen und Physiker daran, unter Nutzung quantenphysikalischer Effekte Materialien mit neuartigen Eigenschaften zu kreieren.
Am Teilchenbeschleuniger in Darmstadt werden die extremen Bedingungen unseres Universums im Labor erforscht. Dabei gelang es den Physikerinnen und Physikern, eine Technologie zu entwickeln, die Energie zur Teilchenbeschleunigung wiederverwendet und einspart. Der Teilchenbeschleuniger ist eingebunden in das Clusterprojekt ELEMENTS, das gemeinsam von der Goethe-Universität Frankfurt und der TU Darmstadt geleitet wird.
Using 10.1 × 109 J/ψ events produced by the Beijing Electron Positron Collider (BEPCII) at a center-of-mass energy √s = 3.097 GeV and collected with the BESIII detector, we present a search for the rare semi-leptonic decay J/ψ → D−e+νe + c.c. No excess of signal above background is observed, and an upper limit on the branching fraction ℬ(J/ψ → D−e+νe + c. c.) < 7.1 × 10−8 is obtained at 90% confidence level. This is an improvement of more than two orders of magnitude over the previous best limit.
We derive the thermal noise spectrum of the longitudinal and transverse electric field operator of a given wave vector starting from the quantum-statistical definitions and relate it to the frequency and wave vector dependent complex conductivity in a homogeneous, isotropic system of electromagnetic interacting charged particles in the frame of the non-relativistic QED. No additional assumptions except the validity of linear response are used in the proof. The Nyquist formula for vanishing frequency, as well as the noise spectral density of Callen-Welton follow as byproduct. Furthermore we discuss also the noise of the photon occupation numbers.
As part of the research for this thesis, a momentum spectrometer was set up and initial measurements on accelerated ions were performed. For this purpose, the necessary hardware for the operation of the spectrometer and for high-precision measurements was were assembled. A control system for remote operation was developed and the spectrometer was installed at the used beamline.
There, measurements of low-energy ion beams in superposition with electrons confined in a Gabor lens can be carried out.
Investigations were made on both the Gabor lens-generated ions and the beam ions, leading to first results regarding the charge changes of beam ions during propagation through an electron atmosphere.
In order to understand the origin of the elements in the universe, one must understand the nuclear reactions by which atomic nuclei are transformed. There are many different astrophysical environments that fulfill the conditions of different nucleosynthesis processes. Even though great progress has been made in recent decades in understanding the origin of the elements in the universe, some questions remain unanswered. In order to understand the processes, it is necessary to measure cross sections of the involved reactions and constrain theoretical model predictions. A variety of methods have been developed to measure nuclear reaction cross sections relevant for nuclear astrophysics. In this thesis, two different experiments and their results, both using the well-established activation method, are presented.
A measurement of the proton capture cross section on the p-nuclide 96Ru was performed at the Institute of Structure and Nuclear Astrophysics ISNAP - Notre Dame, USA. The main goal of this experiment was to compare the results with those obtained by Mei et al. in a pioneering experiment using the method of inverse kinematics at the GSI Helmholtzzentrum für Schwerionenforschung GmbH - Darmstadt, Germany. Therefore, the activations were taken out at the same center of mass energies of 9 MeV, 10 MeV and 11 MeV. Another activation was taken out at an energy of 3.2 MeV to compare the result to a measurement of Bork et al. who also used the activation method. While the results at 3.2 MeV agree quite well with those of Bork et al., the results at higher energies show significantly smaller cross sections than those measured by Mei et al.. Experimental details, the data analysis and sources of uncertainties are discussed.
The second part of this thesis describes a neutron capture cross section experiment. At the Institut für Kernphysik - Goethe Universtität Frankfurt an experimental setup allows to produce quasi maxwell-distributed neutron fields to measure maxwell-averaged cross sections (MACS) relevant for s-process nucleosynthesis. The setup was upgraded by a fast electric linear guide to transport samples from the activation to the detection site. The cyclic activation of the sample allows to increase the signal-to-noise ratio and to measure neutron captures that lead to nuclei with
half-lives on the order of seconds. In a first campaign, MACS of the reactions 51V(n,γ), 107,109Ag(n,γ) and 103Rh(n,γ) were measured. The new components of the setup aswell as the data analysis framework are described and the results of the measurements are discussed.
We study the polarization of relativistic fluids using the relativistic density operator at global and local equilibrium. In global equilibrium, a new technique to compute exact expectation values is introduced, which is used to obtain the exact polarization vector for fields of any spin. The same result has been extended to the case of massless fields. Furthermore, it is demonstrated that at local equilibrium not only the thermal vorticity but also the thermal shear contribute to the polarization vector. It is shown that assuming an isothermal local equilibrium, the new term can solve the polarization sign puzzle in heavy ion collisions.
Using combined data from the Relativistic Heavy Ion and Large Hadron Colliders, we constrain the shear and bulk viscosities of quark-gluon plasma (QGP) at temperatures of ∼150–350 MeV. We use Bayesian inference to translate experimental and theoretical uncertainties into probabilistic constraints for the viscosities. With Bayesian model averaging we propagate an estimate of the model uncertainty generated by the transition from hydrodynamics to hadron transport in the plasma’s final evolution stage, providing the most reliable phenomenological constraints to date on the QGP viscosities.
Tracking influenza a virus infection in the lung from hematological data with machine learning
(2022)
The tracking of pathogen burden and host responses with minimal-invasive methods during respiratory infections is central for monitoring disease development and guiding treatment decisions. Utilizing a standardized murine model of respiratory Influenza A virus (IAV) infection, we developed and tested different supervised machine learning models to predict viral burden and immune response markers, i.e. cytokines and leukocytes in the lung, from hematological data. We performed independently in vivo infection experiments to acquire extensive data for training and testing purposes of the models. We show here that lung viral load, neutrophil counts, cytokines like IFN-γ and IL-6, and other lung infection markers can be predicted from hematological data. Furthermore, feature analysis of the models shows that blood granulocytes and platelets play a crucial role in prediction and are highly involved in the immune response against IAV. The proposed in silico tools pave the path towards improved tracking and monitoring of influenza infections and possibly other respiratory infections based on minimal-invasively obtained hematological parameters.
Intrinsically disordered regions (IDRs) are essential for membrane receptor regulation but often remain unresolved in structural studies. TRPV4, a member of the TRP vanilloid channel family involved in thermo- and osmosensation, has a large N-terminal IDR of approximately 150 amino acids. With an integrated structural biology approach, we analyze the structural ensemble of the TRPV4 IDR and identify a network of regulatory elements that modulate channel activity in a hierarchical lipid-dependent manner through transient long-range interactions. A highly conserved autoinhibitory patch acts as a master regulator by competing with PIP2 binding to attenuate channel activity. Molecular dynamics simulations show that loss of the interaction between the PIP2-binding site and the membrane reduces the force exerted by the IDR on the structured core of TRPV4. This work demonstrates that IDR structural dynamics are coupled to TRPV4 activity and highlights the importance of IDRs for TRP channel function and regulation.
During infection the SARS-CoV-2 virus fuses its viral envelope with cellular membranes of its human host. Initial contact with the host cell and membrane fusion are both mediated by the viral spike (S) protein. Proteolytic cleavage of S at the S2′ site exposes its 40 amino acid long fusion peptide (FP). Binding of the FP to the host membrane anchors the S2 domain of S in both the viral and the host membrane. The reorganization of S2 then pulls the two membranes together. Here we use molecular dynamics (MD) simulations to study the two core functions of the SARS-CoV-2 FP: to attach quickly to cellular membranes and to form an anchor strong enough to withstand the mechanical force during membrane fusion. In eight 10 μs-long MD simulations of FP in proximity to endosomal and plasma membranes, we find that FP binds spontaneously to the membranes and that binding proceeds predominantly by insertion of two short amphipathic helices into the membrane interface. Connected via a flexible linker, the two helices can bind the membrane independently, yet binding of one promotes the binding of the other by tethering it close to the target membrane. By simulating mechanical pulling forces acting on the C-terminus of the FP we then show that the bound FP can bear forces up to 250 pN before detaching from the membrane. This detachment force is more than ten-fold higher than an estimate of the force required to pull host and viral membranes together for fusion. We identify a fully conserved disulfide bridge in the FP as a major factor for the high mechanical stability of the FP membrane anchor. We conclude, first, that the sequential binding of two short amphipathic helices allows the SARS-CoV-2 FP to insert quickly into the target membrane, before the virion is swept away after shedding the S1 domain connecting it to the host cell receptor. Second, we conclude that the double attachment and the conserved disulfide bridge establish the strong anchoring required for subsequent membrane fusion. Multiple distinct membrane-anchoring elements ensure high avidity and high mechanical strength of FP-membrane binding.
Focused ion beam induced deposition (FIBID) is a direct-write technique enabling the growth of individual nanostructures of any shape and dimension with high lateral resolution. Moreover, the fast and reliable writing of periodically arranged nanostructures can be used to fabricate devices for the investigation of collective phenomena and to design novel functional metamaterials. Here, FIBID is employed to prepare dc-Josephson junction arrays (dc-JJA) consisting of superconducting NbC dots coupled through the proximity effect via a granular metal layer. The fabrication is straightforward and allows the preparation of dc-JJA within a few seconds. Microstructure and composition of the arrays are investigated by transmission electron microscopy and energy dispersive X-ray spectroscopy. The superconductor-to-metal transition of the prepared dc-JJA is studied in a direct way, by tuning the Josephson junction resistance in 70 nm-spaced superconducting NbC dots. The observed magnetoresistance oscillations with a period determined by the flux quantum give evidence for the coherent charge transport by paired electrons. Moreover, the measured resistance minima correspond to two fundamental matching configurations of fluxons in the dc-JJA, caused by magnetic frustration. The robust properties of the prepared dc-JJA demonstrate the opportunities for a fast preparation of complex device configurations using direct-write approaches.
Changes in the efficacies of synapses are thought to be the neurobiological basis of learning and memory. The efficacy of a synapse depends on its current number of neurotransmitter receptors. Recent experiments have shown that these receptors are highly dynamic, moving back and forth between synapses on time scales of seconds and minutes. This suggests spontaneous fluctuations in synaptic efficacies and a competition of nearby synapses for available receptors. Here we propose a mathematical model of this competition of synapses for neurotransmitter receptors from a local dendritic pool. Using minimal assumptions, the model produces a fast multiplicative scaling behavior of synapses. Furthermore, the model explains a transient form of heterosynaptic plasticity and predicts that its amount is inversely related to the size of the local receptor pool. Overall, our model reveals logistical tradeoffs during the induction of synaptic plasticity due to the rapid exchange of neurotransmitter receptors between synapses.
Changes in the efficacies of synapses are thought to be the neurobiological basis of learning and memory. The efficacy of a synapse depends on its current number of neurotransmitter receptors. Recent experiments have shown that these receptors are highly dynamic, moving back and forth between synapses on time scales of seconds and minutes. This suggests spontaneous fluctuations in synaptic efficacies and a competition of nearby synapses for available receptors. Here we propose a mathematical model of this competition of synapses for neurotransmitter receptors from a local dendritic pool. Using minimal assumptions, the model produces a fast multiplicative scaling behavior of synapses. Furthermore, the model explains a transient form of heterosynaptic plasticity and predicts that its amount is inversely related to the size of the local receptor pool. Overall, our model reveals logistical tradeoffs during the induction of synaptic plasticity due to the rapid exchange of neurotransmitter receptors between synapses.
Transport of lipids across membranes is fundamental for diverse biological pathways in cells. Multiple ion-coupled transporters participate in lipid translocation, but their mechanisms remain largely unknown. Major facilitator superfamily (MFS) lipid transporters play central roles in cell wall synthesis, brain development and function, lipids recycling, and cell signaling. Recent structures of MFS lipid transporters revealed overlapping architectural features pointing towards a common mechanism. Here we used cysteine disulfide trapping, molecular dynamics simulations, mutagenesis analysis, and transport assays in vitro and in vivo, to investigate the mechanism of LtaA, a proton-dependent MFS lipid transporter essential for lipoteichoic acids synthesis in the pathogen Staphylococcus aureus. We reveal that LtaA displays asymmetric lateral openings with distinct functional relevance and that cycling through outward- and inward-facing conformations is essential for transport activity. We demonstrate that while the entire amphipathic central cavity of LtaA contributes to lipid binding, its hydrophilic pocket dictates substrate specificity. We propose that LtaA catalyzes lipid translocation by a ‘trap-and-flip’ mechanism that might be shared among MFS lipid transporters.
The severity of the COVID-19 pandemic, caused by the SARS-CoV-2 coronavirus, calls for the urgent development of a vaccine. The primary immunological target is the SARS-CoV-2 spike (S) protein. S is exposed on the viral surface to mediate viral entry into the host cell. To identify possible antibody binding sites not shielded by glycans, we performed multi-microsecond molecular dynamics simulations of a 4.1 million atom system containing a patch of viral membrane with four full-length, fully glycosylated and palmitoylated S proteins. By mapping steric accessibility, structural rigidity, sequence conservation and generic antibody binding signatures, we recover known epitopes on S and reveal promising epitope candidates for vaccine development. We find that the extensive and inherently flexible glycan coat shields a surface area larger than expected from static structures, highlighting the importance of structural dynamics in epitope mapping.
Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physicochemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein–protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low-complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS-LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model. For a scaling factor α = 0.65, we recover experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS-LCD condenses. In the region of phase separation, we simulate FUS-LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS-LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01–0.4 mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02 Pa s for the FUS-LCD droplets with scaling factors α in the range of 0.625–0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
The scalar glueball G is the lightest particle of the Yang–Mills sector of QCD, with a lattice predicted mass of about mG≃1.7GeV. It is natural to investigate glueball-glueball scattering and the possible emergence of a bound state, that we call glueballonium. We perform this study in the context of a widely used dilaton potential, that depends on a single dimensionful parameter ΛG. We consider a unitarization prescription that allows us to predict the lowest partial waves in the elastic window. These quantities can be in principle calculated on the lattice, thus offering possibility for testing the validity of the dilaton potential and an independent determination of its parameter. Moreover, we also show that a stable glueballonium exists if ΛG is small enough. In particular, for ΛG compatible with the expectations from the gluon condensate, the glueballonium has a mass of about 3.4GeV.
In physics, the wavefunctions of bosonic particles collapse when the system undergoes a Bose–Einstein condensation. In game theory, the strategy of an agent describes the probability to engage in a certain course of action. Strategies are expected to differ in competitive situations, namely when there is a penalty to do the same as somebody else. We study what happens when agents are interested how they fare not only in absolute terms, but also relative to others. This preference, denoted envy, is shown to induce the emergence of distinct social classes via a collective strategy condensation transition. Members of the lower class pursue identical strategies, in analogy to the Bose–Einstein condensation, with the upper class remaining individualistic.
We present the first very long baseline interferometric (VLBI) observations of the blazar OJ 287 carried out jointly with the Global Millimeter VLBI Array (GMVA) and the phased Atacama Large Millimeter/submillimeter Array (ALMA) at 3.5 mm on 2017 April 2. The participation of phased ALMA has not only improved the GMVA north–south resolution by a factor of ∼3, but has also enabled fringe detections with signal-to-noise ratios up to 300 at baselines longer than 2 Gλ. The high sensitivity has motivated us to image the data with newly developed regularized maximum likelihood imaging methods, revealing the innermost jet structure with unprecedentedly high angular resolution. Our images reveal a compact and twisted jet extending along the northwest direction, with two bends within the inner 200 μas, resembling a precessing jet in projection. The component at the southeastern end shows a compact morphology and high brightness temperature, and is identified as the VLBI core. An extended jet feature that lies at ∼200 μas northwest of the core shows a conical shape, in both total and linearly polarized intensity, and a bimodal distribution of the linear polarization electric vector position angle. We discuss the nature of this feature by comparing our observations with models and simulations of oblique and recollimation shocks with various magnetic field configurations. Our high-fidelity images also enabled us to search for possible jet features from the secondary supermassive black hole (SMBH) and test the SMBH binary hypothesis proposed for this source.
We report the first measurements of the absolute branching fractions of D0 → K0 Lϕ, D0 → K0Lη, D0 → K0Lω, and D0 → K0Lη0, by analyzing 2.93 fb−1 of eþe− collision data taken at a center-of-mass energy of 3.773 GeV with the BESIII detector. Taking the world averages of the branching fractions of D0 → K0Sϕ, D0 → K0Sη, D0 → K0Sω, and D0 → K0Sη0, the K0S − K0L asymmetries RðD0; XÞ in these decay modes are obtained. The CP asymmetries in these decays are also determined. No significant CP violation is observed
In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD)1. These partons subsequently emit further partons in a process that can be described as a parton shower2, which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass mQ and energy E, within a cone of angular size mQ/E around the emitter3. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques4,5 to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.
Modern societies face the challenge that the time scale of opinion formation is continuously accelerating in contrast to the time scale of political decision making. With the latter remaining of the order of the election cycle we examine here the case that the political state of a society is determined by the continuously evolving values of the electorate. Given this assumption we show that the time lags inherent in the election cycle will inevitable lead to political instabilities for advanced democracies characterized both by an accelerating pace of opinion dynamics and by high sensibilities (political correctness) to deviations from mainstream values. Our result is based on the observation that dynamical systems become generically unstable whenever time delays become comparable to the time it takes to adapt to the steady state. The time needed to recover from external shocks grows in addition dramatically close to the transition. Our estimates for the order of magnitude of the involved time scales indicate that socio-political instabilities may develop once the aggregate time scale for the evolution of the political values of the electorate falls below 7–15 months.
We extend the standard solid-state quantum mechanical Hamiltonian containing only Coulomb interactions between the charged particles by inclusion of the (transverse) current-current diamagnetic interaction starting from the non-relativistic QED restricted to the states without photons and neglecting the retardation in the photon propagator. This derivation is supplemented with a derivation of an analogous result along the non-rigorous old classical Darwin-Landau-Lifshitz argumentation within the physical Coulomb gauge.
The rapid spread of the Coronavirus (COVID-19) confronts policy makers with the problem of measuring the effectiveness of containment strategies, balancing public health considerations with the economic costs of social distancing measures. We introduce a modified epidemic model that we name the controlled-SIR model, in which the disease reproduction rate evolves dynamically in response to political and societal reactions. An analytic solution is presented. The model reproduces official COVID-19 cases counts of a large number of regions and countries that surpassed the first peak of the outbreak. A single unbiased feedback parameter is extracted from field data and used to formulate an index that measures the efficiency of containment strategies (the CEI index). CEI values for a range of countries are given. For two variants of the controlled-SIR model, detailed estimates of the total medical and socio-economic costs are evaluated over the entire course of the epidemic. Costs comprise medical care cost, the economic cost of social distancing, as well as the economic value of lives saved. Under plausible parameters, strict measures fare better than a hands-off policy. Strategies based on current case numbers lead to substantially higher total costs than strategies based on the overall history of the epidemic.
Though immensely successful, the standard model of particle physics does not offer any explanation as to why our Universe contains so much more matter than antimatter. A key to a dynamically generated matter–antimatter asymmetry is the existence of processes that violate the combined charge conjugation and parity (CP) symmetry1. As such, precision tests of CP symmetry may be used to search for physics beyond the standard model. However, hadrons decay through an interplay of strong and weak processes, quantified in terms of relative phases between the amplitudes. Although previous experiments constructed CP observables that depend on both strong and weak phases, we present an approach where sequential two-body decays of entangled multi-strange baryon–antibaryon pairs provide a separation between these phases. Our method, exploiting spin entanglement between the double-strange Ξ− baryon and its antiparticle2 Ξ¯+
, has enabled a direct determination of the weak-phase difference, (ξP − ξS) = (1.2 ± 3.4 ± 0.8) × 10−2 rad. Furthermore, three independent CP observables can be constructed from our measured parameters. The precision in the estimated parameters for a given data sample size is several orders of magnitude greater than achieved with previous methods3. Finally, we provide an independent measurement of the recently debated Λ decay parameter αΛ (refs. 4,5). The ΛΛ¯
asymmetry is in agreement with and compatible in precision to the most precise previous measurement.
Living cells constantly remodel the shape of their lipid membranes. In the endo-plasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attraction. At a critical size, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetics barrier for membrane remodeling.
Living cells constantly remodel the shape of their lipid membranes. In the endoplasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attractions. At a critical size, as in a nucleation process, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetic barrier for membrane remodeling.
Gasdermin-D (GSDMD) is the ultimate effector of pyroptosis, a form of programmed cell death associated with pathogen invasion and inflammation. After proteolytic cleavage by caspases activated by the inflammasome, the GSDMD N-terminal domain (GSDMDNT) assembles on the inner leaflet of the plasma membrane and induces the formation of large membrane pores. We use atomistic molecular dynamics simulations to study GSDMDNT monomers, oligomers, and rings in an asymmetric plasma membrane mimetic. We identify distinct interaction motifs of GSDMDNT with phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) and phosphatidylserine (PS) head-groups and describe differential lipid binding between the pore and prepore conformations. Oligomers are stabilized by shared lipid binding sites between neighboring monomers acting akin to double-sided tape. We show that already small GSDMDNT oligomers form stable, water-filled and ion-conducting membrane pores bounded by curled beta-sheets. In large-scale simulations, we resolve the process of pore formation by lipid detachment from GSDMDNT arcs and lipid efflux from partial rings. We find that that high-order GSDMDNT oligomers can crack under the line tension of 86 pN created by an open membrane edge to form the slit pores or closed GSDMDNT rings seen in experiment. Our simulations provide a detailed view of key steps in GSDMDNT-induced plasma membrane pore formation, including sublytic pores that explain nonselective ion flux during early pyroptosis.
Gasdermin-D (GSDMD) is the ultimate effector of pyroptosis, a form of programmed cell death associated with pathogen invasion and inflammation. After proteolytic cleavage by caspases, the GSDMD N-terminal domain (GSDMDNT) assembles on the inner leaflet of the plasma membrane and induces the formation of membrane pores. We use atomistic molecular dynamics simulations to study GSDMDNT monomers, oligomers, and rings in an asymmetric plasma membrane mimetic. We identify distinct interaction motifs of GSDMDNT with phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) and phosphatidylserine (PS) headgroups and describe their conformational dependence. Oligomers are stabilized by shared lipid binding sites between neighboring monomers acting akin to double-sided tape. We show that already small GSDMDNT oligomers support stable, water-filled, and ion-conducting membrane pores bounded by curled beta-sheets. In large-scale simulations, we resolve the process of pore formation from GSDMDNT arcs and lipid efflux from partial rings. We find that high-order GSDMDNT oligomers can crack under the line tension of 86 pN created by an open membrane edge to form the slit pores or closed GSDMDNT rings seen in atomic force microscopy experiments. Our simulations provide a detailed view of key steps in GSDMDNT-induced plasma membrane pore formation, including sublytic pores that explain nonselective ion flux during early pyroptosis.
Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport. Their intricate 120 MDa architecture remains incompletely understood. Here, we report a near-complete structural model of the human NPC scaffold with explicit membrane and in multiple conformational states. We combined AI-based structure prediction with in situ and in cellulo cryo-electron tomography and integrative modeling. We show that linker Nups spatially organize the scaffold within and across subcomplexes to establish the higher-order structure. Microsecond-long molecular dynamics simulations suggest that the scaffold is not required to stabilize the inner and outer nuclear membrane fusion, but rather widens the central pore. Our work exemplifies how AI-based modeling can be integrated with in situ structural biology to understand subcellular architecture across spatial organization levels.
Direct laser acceleration (DLA) of electrons in a plasma of near-critical electron density (NCD) and the associated synchrotron-like radiation are discussed for moderate relativistic laser intensity (normalized laser amplitude a0 ≤ 4.3) and ps length pulse. This regime is typical of kJ PW-class laser facilities designed for high-energy-density (HED) research. In experiments at the PHELIX facility, it has been demonstrated that interaction of a 1019 W/cm2 sub-ps laser pulse with a sub-mm length NCD plasma results in the generation of high-current well-directed super-ponderomotive electrons with an effective temperature ten times higher than the ponderomotive potential [Rosmej et al., Plasma Phys. Controlled Fusion 62, 115024 (2020)]. Three-dimensional particle-in-cell simulations provide good agreement with the measured electron energy distribution and are used in the current work to study synchrotron radiation from the DLA-accelerated electrons. The resulting x-ray spectrum with a critical energy of 5 keV reveals an ultrahigh photon number of 7 × 1011 in the 1–30 keV photon energy range at the focused laser energy of 20 J. Numerical simulations of betatron x-ray phase contrast imaging based on the DLA process for the parameters of a PHELIX laser are presented. The results are of interest for applications in HED experiments, which require a ps x-ray pulse and a high photon flux.
Disordered proteins and nucleic acids can condense into droplets that resemble the membraneless organelles observed in living cells. MD simulations offer a unique tool to characterize the molecular interactions governing the formation of these biomolecular condensates, their physico-chemical properties, and the factors controlling their composition and size. However, biopolymer condensation depends sensitively on the balance between different energetic and entropic contributions. Here, we develop a general strategy to fine-tune the potential energy function for molecular dynamics simulations of biopolymer phase separation. We rebalance protein-protein interactions against solvation and entropic contributions to match the excess free energy of transferring proteins between dilute solution and condensate. We illustrate this formalism by simulating liquid droplet formation of the FUS low complexity domain (LCD) with a rebalanced MARTINI model. By scaling the strength of the nonbonded interactions in the coarse-grained MARTINI potential energy function, we map out a phase diagram in the plane of protein concentration and interaction strength. Above a critical scaling factor of αc ≈ 0.6, FUS LCD condensation is observed, where α = 1 and 0 correspond to full and repulsive interactions in the MARTINI model, respectively. For a scaling factor α = 0.65, we recover the experimental densities of the dilute and dense phases, and thus the excess protein transfer free energy into the droplet and the saturation concentration where FUS LCD condenses. In the region of phase separation, we simulate FUS LCD droplets of four different sizes in stable equilibrium with the dilute phase and slabs of condensed FUS LCD for tens of microseconds, and over one millisecond in aggregate. We determine surface tensions in the range of 0.01 to 0.4mN/m from the fluctuations of the droplet shape and from the capillary-wave-like broadening of the interface between the two phases. From the dynamics of the protein end-to-end distance, we estimate shear viscosities from 0.001 to 0.02Pas for the FUS LCD droplets with scaling factors α in the range of 0.625 to 0.75, where we observe liquid droplets. Significant hydration of the interior of the droplets keeps the proteins mobile and the droplets fluid.
The protein Atg2 has been proposed to form a membrane tether that mediates lipid transfer from the ER to the phagophore in autophagy. However, recent kinetic measurements on the human homolog ATG2A indicated a transport rate of only about one lipid per minute, which would be far too slow to deliver the millions of lipids required to form a phagophore on a physiological time scale. Here, we revisit the analysis of the fluorescence quenching experiments. We develop a detailed kinetic model of the lipid transfer between two membranes bridged by a tether that forms a conduit for lipids. The model provides an excellent fit to the fluorescence experiments, with a lipid transfer rate of about 100 per second and protein. At this rate, Atg2-mediated transfer can supply a significant fraction of the lipids required in autophagosome biogenesis. Our kinetic model is generally applicable to lipid-transfer experiments, in particular to proteins forming organelle contact sites in cells.
Binding of the spike protein of SARS-CoV-2 to the human angiotensin-converting enzyme 2 (ACE2) receptor triggers translocation of the virus into cells. Both the ACE2 receptor and the spike protein are heavily glycosylated, including at sites near their binding interface. We built fully glycosylated models of the ACE2 receptor bound to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. Using atomistic molecular dynamics (MD) simulations, we found that the glycosylation of the human ACE2 receptor contributes substantially to the binding of the virus. Interestingly, the glycans at two glycosylation sites, N90 and N322, have opposite effects on spike protein binding. The glycan at the N90 site partly covers the binding interface of the spike RBD. Therefore, this glycan can interfere with the binding of the spike protein and protect against docking of the virus to the cell. By contrast, the glycan at the N322 site interacts tightly with the RBD of the ACE2-bound spike protein and strengthens the complex. Remarkably, the N322 glycan binds to a conserved region of the spike protein identified previously as a cryptic epitope for a neutralizing antibody. By mapping the glycan binding sites, our MD simulations aid in the targeted development of neutralizing antibodies and SARS-CoV-2 fusion inhibitors.
Binding of the spike protein of SARS-CoV-2 to the human angiotensin converting enzyme 2 (ACE2) receptor triggers translocation of the virus into cells. Both the ACE2 receptor and the spike protein are heavily glycosylated, including at sites near their binding interface. We built fully glycosylated models of the ACE2 receptor bound to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. Using atomistic molecular dynamics (MD) simulations, we found that the glycosylation of the human ACE2 receptor contributes substantially to the binding of the virus. Interestingly, the glycans at two glycosylation sites, N90 and N322, have opposite effects on spike protein binding. The glycan at the N90 site partly covers the binding interface of the spike RBD. Therefore, this glycan can interfere with the binding of the spike protein and protect against docking of the virus to the cell. By contrast, the glycan at the N322 site interacts tightly with the RBD of the ACE2-bound spike protein and strengthens the complex. Remarkably, the N322 glycan binds into a conserved region of the spike protein identified previously as a cryptic epitope for a neutralizing antibody. By mapping the glycan binding sites, our MD simulations aid in the targeted development of neutralizing antibodies and SARS-CoV-2 fusion inhibitors.
The electrical and computational properties of neurons in our brains are determined by a rich repertoire of membrane-spanning ion channels and elaborate dendritic trees. However, the precise reason for this inherent complexity remains unknown. Here, we generated large stochastic populations of biophysically realistic hippocampal granule cell models comparing those with all 15 ion channels to their reduced but functional counterparts containing only 5 ion channels. Strikingly, valid parameter combinations in the full models were more frequent and more stable in the face of perturbations to channel expression levels. Scaling up the numbers of ion channels artificially in the reduced models recovered these advantages confirming the key contribution of the actual number of ion channel types. We conclude that the diversity of ion channels gives a neuron greater flexibility and robustness to achieve target excitability.
We introduce a novel approach based on elas- tic and inelastic scattering rates to extract the hyper-surface of the chemical freeze-out from a hadronic transport model in the energy range from Elab = 1.23 AGeV to √sNN = 62.4 GeV. For this study, the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model combined with a coarse-graining method is employed. The chemical freeze- out distribution is reconstructed from the pions through sev- eral decay and re-formation chains involving resonances and taking into account inelastic, pseudo-elastic and string excita- tion reactions. The extracted average temperature and baryon chemical potential are then compared to statistical model analysis. Finally we investigate various freeze-out criteria suggested in the literature. We confirm within this micro- scopic dynamical simulation, that the chemical freeze-out at all energies coincides with ⟨E⟩/⟨N⟩ ≈ 1 GeV, while other criteria, like s/T 3 = 7 and nB +nB ̄ ≈ 0.12 fm−3 are limited to higher collision energies.
The production of light (anti-)nuclei and (anti-)hypernuclei in ultra-relativistic heavy-ion collisions, but also in more elementary collisions as proton–proton and proton–nucleus collisions, became recently a focus of interest. In particular, the fact that these objects are all loosely bound compared to the temperature and energies, e.g. the kinetic energies involved, is often stressed out to be special for their production. The binding energies of these (anti-)nuclei is between 130 keV (Λ separation energy in the hypertriton) and about 8 MeV per nucleon. Whereas the connected temperatures are of the order of 100 to 160 MeV. This lead to some difficulties in the interpretation of the usually discussed production models, i.e. coalescence and statistical-thermal models, as will be discussed here. In this brief review we discuss selected highlights of the production of light (anti-)nuclei, such as (anti-)deuteron, (anti-)helium and (anti-)alpha nuclei. In addition, we will discuss the current status of the highly debated lifetime of the (anti-)hypertriton and connected measurements and model results.
We investigate hadronic particle spectra and flow characteristics of heavy-ion reactions in the FAIR/NICA energy range of 1 AGeV ≤ Elab ≤ 10 AGeV within a relativistic ideal hydrodynamic one-fluid approach. The particlization is realized by sampling the Cooper-Frye distribution for a grand canonical hadron gas on a hypersurface of constant energy density. Results of the hydrodynamic calculations for different underlying equations of state are presented and compared with experimental data and microscopic transport simulations. The sensitivity of the approach to physical model inputs concerning the initial state and the particlization is studied.
Terahertz (THz) radiation lies between the micro and far-infrared range in the electromagnetic spectrum. Compared with microwave and millimeter waves, it has a larger signal bandwidth and extremely narrow antenna beam. Thus, it is easier to achieve high-resolution for imaging and detection applications. The unique properties, such as penetration for majority non-polar materials, non-ionizing characteristic and the spectral fingerprint of materials, makes THz imaging an appealing artifice in the military, biomedical, astronomical communications, and other areas. However, THz radiation’s current low power level and detection sensitivity block THz imaging system from including fewer optical elements than the visible or infrared range. This leads to imaging resolution, contrast, and imaging field of view degenerate and makes the aberration more serious. THz imaging based on the space Fourier spectrum detection is developed in this thesis to achieve high-quality imaging. The main concept of Fourier imaging is by recording the field distribution in the Fourier plane (focal plane) of the imaging system; the information of the target is obtained. The numerical processing method is needed to extract the amplitude and phase information of the imaged target. With additional process, three-dimensional (3D) information can be obtained based on the phase information. The novel recording and reconstructing ways of the Fourier imaging system enables it to have a higher resolution, better contrast, and broader field of view than conventional imaging systems such as microscopy and plane to plane telescopic imaging system.
The work presented in this thesis consists of two imaging systems, one is working at 300 GHz based on the fundamental heterodyne detection of the THz radiation, the other is operated at 600 GHz by utilizing the sub harmonic heterodyne detection technique. The realization and test of the heterodyne detection are based on the THz antenna-coupled field-effect transistor (TeraFET) detector developed by Dr. Alvydas Lisauskas. Both systems use two synchronized electronic multiplier chains to radiate the THz waves. One radiation works as the local oscillator (LO), the other works as illumination with a slight frequency shift, the radiations are mixed on the detector scanning in the Fourier plane to record the complex Fourier spectrum of the imaged target. The LO has the same frequency range as the illuminating radiation for fundamental heterodyne detection but half the frequency range for the sub-harmonic heterodyne detection. The 2-mm resolution, 60-dB contrast, and 5.5-cm diameter imaging area at 300 GHz and the of 500-μm resolution, 40-dB contrast, and 3.5-cm diameter imaging area at 600 GHz are achieved (the 300-GHz illuminating radiation has the approximate power of 600 μW , the 600-GHz illuminating radiation has the approximate power of 60 μW ).
The thesis consists of 6 parts. After the introduction, the second chapter expands on the topic of Fourier optics from a theoretical point of view and the simulations of the Fourier imaging system. First, the theory of the electromagnetic field propagation in free space and through an optical system are investigated to elicit the Fourier transform function of the imaging system. The simulation is used for theoretical considerations and the implementation of a Fourier optic script that allows for numerical investigations on reconstruction. The preliminary imaging field of view and resolution are also demonstrated. The third chapter describes the Fourier imaging system at 300 GHz based on the fundamental heterodyne detection, including the experimental setup, the 2D, and 3D imaging results. The following fourth chapter reports the integration of the TeraFET detector with two substrate lenses (one is a Si lens on the back-side Si substrate, the other is a wax/PTFE lens on the front side containing the bonding wires) for sub-harmonic heterodyne detection at 600 GHz. The characteristic of the wax/PTFE lens at THz range is presented. After that, the compared imaging results between the detector with and without the wax/PTFE lens are shown. The fifth chapter extends the demonstration on the lateral and depth resolution of the Fourier imaging system in detail and uses the experimental results at 600 GHz to validate the analytical predictions. The comparison of the resolution between the Fourier imaging system and the conventional microscopy system proves that the Fourier imaging system has better imaging quality under the same system configuration. The last chapter in this thesis concludes on the findings of the THz Fourier imaging and gives an outlook for the enhancement of the Fourier imaging system at THz range.
Using 448.1 × 106 ψ(3686) decays collected with the BESIII detector at the BEPCII e+e− storage rings, the branching fractions and angular distributions of the decays χcJ → Ξ−Ξ¯¯¯¯+ and Ξ0Ξ¯¯¯¯0 (J = 0, 1, 2) are measured based on a partial-reconstruction technique. The decays χc1 → Ξ0Ξ¯¯¯¯0 and χc2 → Ξ0Ξ¯¯¯¯0 are observed for the first time with statistical significances of 7σ and 15σ, respectively. The results of this analysis are in good agreement with previous measurements and have significantly improved precision.
Radiative transition of an excited baryon to a nucleon with emission of a virtual massive photon converting to dielectron pair (Dalitz decays) provides important information about baryon-photon coupling at low q2 in timelike region. A prominent enhancement in the respective electromagnetic transition Form Factors (etFF) at q2 near vector mesons ρ/ω poles has been predicted by various calculations reflecting strong baryon-vector meson couplings. The understanding of these couplings is also of primary importance for the interpretation of the emissivity of QCD matter studied in heavy ion collisions via dilepton emission. Dedicated measurements of baryon Dalitz decays in proton-proton and pion-proton scattering with HADES detector at GSI/FAIR are presented and discussed. The relevance of these studies for the interpretation of results obtained from heavy ion reactions is elucidated on the example of the HADES results.
Background: The photon strength functions (PSFs) and nuclear level density (NLD) are key ingredients for calculation of the photon interaction with nuclei, in particular the reaction cross sections. These cross sections are important especially in nuclear astrophysics and in the development of advanced nuclear technologies.
Purpose: The role of the scissors mode in the M1 PSF of (well-deformed) actinides was investigated by several experimental techniques. The analyses of different experiments result in significant differences, especially on the strength of the mode. The shape of the low-energy tail of the giant electric dipole resonance is uncertain as well. In particular, some works proposed a presence of the E1 pygmy resonance just above 7 MeV. Because of these inconsistencies additional information on PSFs in this region is of great interest.
Methods: The γ-ray spectra from neutron-capture reactions on the 234U, 236 U, and 238 U nuclei have been measured with the total absorption calorimeter of the n_TOF facility at CERN. The background-corrected sum-energy and multi-step-cascade spectra were extracted for several isolated s-wave resonances up to about 140 eV.
Results: The experimental spectra were compared to statistical model predictions coming from a large selection of models of photon strength functions and nuclear level density. No combination of PSF and NLD models from literature is able to globally describe our spectra. After extensive search we were able to find model combinations with modified generalized Lorentzian (MGLO) E1 PSF, which match the experimental spectra as well as the total radiative widths.
Conclusions: The constant temperature energy dependence is favored for a NLD. The tail of giant electric dipole resonance is well described by the MGLO model of the E1 PSF with no hint of pygmy resonance. The M1 PSF must contain a very strong, relatively wide, and likely double-resonance scissors mode. The mode is responsible for about a half of the total radiative width of neutron resonances and significantly affects the radiative cross section.
We present the first experimental search for the rare charm decay D0→π0ν¯ν. It is based on an e+e− collision sample consisting of 10.6×10^6 pairs of D0¯D0 mesons collected by the BESIII detector at √s=3.773 GeV, corresponding to an integrated luminosity of 2.93 fb^−1. A data-driven method is used to ensure the reliability of the background modeling. No significant D0→π0ν¯ν signal is observed in data and an upper limit of the branching fraction is set to be 2.1×10^-4 at the 90% confidence level. This is the first experimental constraint on charmed-hadron decays into dineutrino final states.