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
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.
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.
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.
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.
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.
Polarization of Λ and ¯Λ hyperons along the beam direction in Pb-Pb collisions at √sNN=5.02 TeV
(2022)
The polarization of the Λ and ¯Λ hyperons along the beam (z) direction, Pz, has been measured in Pb-Pb collisions at √sNN=5.02 TeV recorded with ALICE at the Large Hadron Collider (LHC). The main contribution to Pz comes from elliptic flow-induced vorticity and can be characterized by the second Fourier sine coefficient Pz,s2=⟨Pzsin(2φ−2Ψ2)⟩, where φ is thhyperon azimuthal emission angle and Ψ2 is the elliptic flow plane angle. We report the measurement of Pz,s2 for different collision centralities and in the 30%–50% centrality interval as a function of the hyperon transverse momentum and rapidity. The Pz,s2 is positive similarly as measured by the STAR Collaboration in Au-Au collisions at √sNN=200 GeV, with somewhat smaller amplitude in the semicentral collisions. This is the first experimental evidence of a nonzero hyperon Pz in Pb-Pb collisions at the LHC. The comparison of the measured Pz,s2 with the hydrodynamic model calculations shows sensitivity to the competing contributions from thermal and the recently found shear-induced vorticity, as well as to whether the polarization is acquired at the quark-gluon plasma or the hadronic phase.
Purpose: A clinical implementation of ion-beam radiography (iRad) is envisaged to provide a method for on-couch verification of ion-beam treatment plans. The aim of this work is to introduce and evaluate a method for quantitative water-equivalent thickness (WET) measurements for a specific helium-ion imaging system for WETs that are relevant for imaging thicker body parts in the future.
Methods: Helium-beam radiographs (αRads) are measured at the Heidelberg Ion-beam Therapy Center with an initial beam energy of 239.5 MeV/u. An imaging system based on three pairs of thin silicon pixel detectors is used for ion path reconstruction and measuring the energy deposition (dE) of each particle behind the object to be imaged. The dE behind homogeneous plastic blocks is related to their well-known WETs between 280.6 and 312.6 mm with a calibration curve that is created by a fit to measured data points. The quality of the quantitative WET measurements is determined by the uncertainty of the measured WET of a single ion (single-ion WET precision) and the deviation of a measured WET value to the well-known WET (WET accuracy). Subsequently, the fitted calibration curve is applied to an energy deposition radiograph of a phantom with a complex geometry. The spatial resolution (modulation transfer function at 10 % —MTF10%) and WET accuracy (mean absolute percentage difference—MAPD) of the WET map are determined.
Results: In the optimal imaging WET-range from ∼280 to 300 mm, the fitted calibration curve reached a mean single-ion WET precision of 1.55
0.00%. Applying the calibration to an ion radiograph (iRad) of a more complex WET distribution, the spatial resolution was determined to be MTF10% = 0.49 0.03 lp/mm and the WET accuracy was assessed as MAPD to 0.21 %.
Conclusions: Using a beam energy of 239.5 MeV/u and the proposed calibration procedure, quantitative αRads of WETs between ∼280 and 300 mm can be measured and show high potential for clinical use. The proposed approach with the resulting image qualities encourages further investigation toward the clinical application of helium-beam radiography.
Die vorliegende Arbeit befasst sich mit der Untersuchung der Transporteigenschaften inklusive Ladungsträgerdynamik von quasi-zweidimensionalen organischen Ladungstransfersalzen. Diese Materialien besitzen eine Schichtstruktur und weisen eine hohe Anisotropie der elektrischen Leitfähigkeit auf. Aufgrund der geringen Bandbreite und der niedrigen Ladungsträgerkonzentration gehören die Materialien zu den stark-korrelierten Elektronensystemen, wobei sich die elektronischen Eigenschaften leicht durch chemische Modifikationen oder äußere Parameter beeinflussen lassen. Die starken Korrelationen resultieren in Metall-Isolator-Übergängen, die sich beim Mott-isolierenden Zustand in einer homogenen Verteilung und beim ladungsgeordneten Zustand in einer periodischen Anordnung der lokalisierten Ladungsträger manifestieren.
Mithilfe der Fluktuationsspektroskopie, die sich mit der Analyse der zeitabhängigen Widerstandsfluktuationen befasst, konnten im Rahmen dieser Arbeit neue Erkenntnisse über die Ladungsträgerdynamik in den verschiedenen elektronischen Zuständen gewonnen werden. Die Metall-Isolator-Übergänge in den untersuchten Systemen, die auf den Molekülen BEDT-TTF (kurz: ET) bzw. BEDT-TSF (kurz: BETS) basieren, sind von der Stärke der strukturellen Dimerisierung abhängig und wurden durch die Kühlrate, eine Zugbelastung sowie durch die Ausnutzung des Feldeffekts beeinflusst.
In den Systemen κ-(BETS)₂Mn[N(CN)₂]₃, κ-(ET)₂Hg(SCN)₂Cl und κ-(ET)₂Cu[N(CN)₂]Br sind die Donormoleküle als Dimere angeordnet, sodass aufgrund der effektiv halben Bandfüllung bei genügender Korrelationsstärke häufig ein Mott-Übergang auftritt. In κ-(ET)₂Hg(SCN)₂Cl führt eine schwächere Dimerisierung jedoch zu einem Ladungsordnungsübergang, der mit elektronischer Ferroelektrizität einhergeht. Dabei wird die polare Ordnung durch eine Ladungsdisproportionierung innerhalb der Dimere verursacht. Die Widerstandsfluktuationen zeigen am ferroelektrischen Übergang einen starken Anstieg der spektralen Leistungsdichte, eine Abhängigkeit vom angelegten elektrischen Feld sowie Zeitabhängigkeiten, die auf räumliche Korrelationen der fluktuierenden Prozesse hindeuten. Diese Eigenschaften wurden ebenfalls für das System κ-(BETS)₂Mn[N(CN)₂]₃ beobachtet. Hierbei wurden mithilfe der dielektrischen Spektroskopie ebenfalls Hinweise auf Ferroelektrizität gefunden, während durch die Analyse der stromabhängigen Widerstandsfluktuationen die Größe der polaren Regionen abgeschätzt werden konnte. Das System κ-(ET)₂Cu[N(CN)₂]Br, das in einer Feldeffekttransistor-Struktur vorliegt, erlaubt neben der Untersuchung des Bandbreiten-getriebenen Mott-Übergangs durch die Zugbelastung eines Substrats auch die Beeinflussung der elektronischen Eigenschaften durch die Änderung der Bandfüllung mittels elektrostatischer Dotierung. Hierbei wurden starke Abhängigkeiten des Widerstands von der Gatespannung beobachtet und Ähnlichkeiten der Ladungsträgerdynamik zu herkömmlichen Volumenproben gefunden.
Bei den Systemen θ-(ET)₂MM'(SCN)₄ mit MM'=CsCo, RbZn, TlZn tritt ein Ladungsordnungsübergang auf, der eine starke Abhängigkeit von der Kühlrate zeigt. Durch schnelles Abkühlen lässt sich der Phasenübergang erster Ordnung kinetisch vermeiden, wodurch ein Ladungsglaszustand realisiert wird. Dieser metastabile Zustand zeigt neuartige physikalische Eigenschaften mit Ähnlichkeiten zu herkömmlichen Gläsern und wurde als Folge der geometrischen Frustration der Ladung auf einem Dreiecksgitter diskutiert. Im Rahmen dieser Arbeit konnte die Ladungsträgerdynamik in den verschiedenen Ladungszuständen von unterschiedlich frustrierten Systemen verglichen werden. Zur Realisierung sehr schneller Abkühlraten wurde dafür eine Heizpulsmethode verwendet und weiterentwickelt. Der Ladungsglaszustand zeigte dabei für verschiedene Systeme ein deutlich niedrigeres Rauschniveau als der ladungsgeordnete Zustand. In Kombination mit Messungen der thermischen Ausdehnung und kühlratenabhängiger Transportmessungen wurde in den Systemen mit der stärksten Frustration die Existenz eines strukturellen Glasübergangs nachgewiesen, der von einer starken Verlangsamung der Ladungsträgerdynamik begleitet wird. Diese Erkenntnisse werfen ein neues Licht auf die bisherige rein elektronische Interpretation des Ladungsglaszustands und heben den Einfluss der strukturellen Freiheitsgrade hervor.
This work ties in with the investigation of the intermediate valent states and valence fluctuations in certain europium based intermetallic systems. Valence fluctuations are a property of the electronic system of a compound that is possibly accompanied by structural effects, which, in some cases, are quite noticable. By assuming how the changes in the electronic system and in the crystal lattice are connected, valence _uctuations of europium are believed to be a possible probe for the theory of quantum critical elasticity, which is investigated on by the SFB TRR 288 (Frankfurt, Mainz, Karlsruhe, Bochum, Dresden).
Here, the proceedings in growing single crystals of di_erent compounds related to this _eld of research are reported. This includes the ThCr2Si2 (122) type compounds EuPd2Si2 as well as the doping series EuPd2(Si1-xGex)2, the Europium based ternary Phosphides EuFe2P2, EuCo2P2, EuNi2P2 and EuRu2P2, and attempts to grow compounds of a derived 1144 structure by ordered substitution of half the Europium, EuKRu4P4.
The largest part of this work focusses on the EuPd2Si2 system, which exhibits intermediate valent europium and a temperature dependent transition between two di_erent intermediate valent states of europium. Crystals of this system were grown using the Czochralski method with a levitating melt and an europium excess flux after a two step prereaction process. Also, explorations of a PdSi-rich flux and external flux methods are reported. Ten Czochralski grown experiments, in six generations iteratevely seeded by the previous generation, were prepared.
Thermodynamical and structural analyses of the crystals located the transition between the di_erent intermediate valent states of europium between 140K and 165 K, transitioning from a high temperature Eu2.3+ state to a low temperature Eu2.7+ state, and classified it as a second order transition. To this transition a lattice anomaly of the a-parameter collapsing about 2% is connected, while the c-parameter remains largely unaffected. Large differences between individual samples can be explained by combining thermodynamical and structural analyses with compositional analysis, revealing the valence transition temperature as strongly dependent on the sample composition and Pd-Si site interchanges.
Searching to change the character of the valence transition to first order, silicon was substituted by germanium to introduce negative pressure. Germanium substituted samples of EuPd2(Si1-xGex)2 were grown using the Czochralski method with the optimized parameters from the growth experiments for the undoped compound. Samples were prepared with a nominal substitution of x = 0.05, x = 0.10, x = 0.15, x = 0.20 (twice) and x = 0.30. For the EuPd2(Si1-xGex)2 system, a phase diagram for the europium valence states is derived from chemical and thermodynamical characterizations.
n ternary europium phosphides EuT2P2, the position of the compounds in the generalized phase diagram and the question of long range magnetic order or valence transition appear connected to an isostructural transition of the tetragonal crystal structure, drastically decreasing the length of the c-parameter while establishing covalent bonds between phosphorus atoms of different interlayers of the structure, the so called ‚collapse‘. While EuFe2P2, EuT2P2 and EuCo2P2 display both long range magnetic order and a non-collapsed crystal structure, EuNi2P2 shows both a valence transition between two intermediate valent states at a characteristic temperature of 36K - accompanied by a small lattice anomaly of the a-parameter shrinking about 0.2% - and a collapsed crystal structure. Samples of EuFe2P2, EuCo2P2 and EuNi2P2 were grown in tin flux and using solid-solid sintering approaches.
Single crystals of EuFe2P2, EuCo2P2 and EuRu2P2 were investigated at ESRF in Grenoble with single crystal X-ray di_ractometry on a pressure range up to 15GPa and at temperatures down to 15K to investigate the nature of the structural transitions in the compounds. While in EuCo2P2 the structural transition occurs as a transition of first order at all temperatures (e.g. at 2GPa for 15 K), in EuFe2P2 and EuRu2P2 the structural collapse evolves over a broad pressure range up to 8GPa and as a transition of second order troughout the temperature ranges, albeit seeming to sharpen at lower temperatures. From the crystallographic data, elastic constants of the compounds could be derived, revealing EuFe2P2 and EuRu2P2 as unexpectedly elastic materials.
In order to probe the structural collapse at more accessible pressures, crystals with a sturcture derived from the 122 structure, but with ordered 50% substitution of europium and hence altering the symmetry from I4/mmm to P4/mmm in a 1144 structure, were exploratively pursued. Different experiments to obtain EuAT4P4 (with A = K, Rb, Cs and T = Fe, Ru) from binary or ternary prereactants or directly from the elements remained largely unsuccessful.
Exciting new scientific opportunities are presented for the PANDA detector at the High Energy Storage Ring in the redefined p¯¯¯p(A) collider mode, HESR-C, at the Facility for Antiproton and Ion Research (FAIR) in Europe. The high luminosity, L∼1031 cm−2 s−1, and a wide range of intermediate and high energies, sNN−−−√ up to 30 GeV for p¯¯¯p(A) collisions will allow to explore a wide range of exciting topics in QCD, including the study of the production of excited open charm and bottom states, nuclear bound states containing heavy (anti)quarks, the interplay of hard and soft physics in the dilepton production, probing short-range correlations in nuclei, and the exploration of the early, complete p¯¯¯-p- annihilation phase, where an initially pure Yang–Mills gluon plasma is formed.
High-resolution, compactness, scalability, efficiency – these are the critical requirements which imaging radar systems have to fulfil in applications such as environmental monitoring, cloud mapping, body sensing or autonomous driving. This thesis presents a modular millimetre-wave frequency modulated continuous-wave (FMCW) radar front-end solution intended for such applications. High-resolution is achieved by enlarging the operating frequency band of the radar system. This can be realized at millimetre-wave frequencies due to the large spectrum availability. Furthermore, the size of components decreasing with increasing frequency makes millimetre-wave systems a good candidate for compactness. However, the full integration of radar front-ends is a challenge at millimetre-wave frequencies due to poor signal integrity and spectral purity, which are essential for imaging applications. The proposed radar uses an alternative technique and tackles this limitation by featuring highly-integrable architectures, specifically the Hartley architecture for signal conversion and enhanced push-pull amplifier for harmonic suppression. The resolution of imaging radars can be further improved by increasing the number of transmitters and receivers. This has spurred the investigation of spectrum, time and energy-efficient multiplexing techniques for multi-input multi-output (MIMO) radar systems. The FMCW radar architecture proposed in this thesis is based on code-division technique using intra-pulse, also called intra-chirp modulation. This advanced scalable and non-complex solution, made possible by the latest achievements on direct digital synthesis for signal generation, guarantees signal integrity and compact size implementation. The proposed architecture is investigated by a thorough system analysis. A transmitter module and a receiver module for a 35 GHz imaging radar prototype are designed, fabricated and fully characterized to validate the feasibility of our novel approach for high-resolution highly-integrated MIMO front-ends.
Cysteine cross-linking in native membranes establishes the transmembrane architecture of Ire1
(2021)
The ER is a key organelle of membrane biogenesis and crucial for the folding of both membrane and secretory proteins. Sensors of the unfolded protein response (UPR) monitor the unfolded protein load in the ER and convey effector functions for maintaining ER homeostasis. Aberrant compositions of the ER membrane, referred to as lipid bilayer stress, are equally potent activators of the UPR. How the distinct signals from lipid bilayer stress and unfolded proteins are processed by the conserved UPR transducer Ire1 remains unknown. Here, we have generated a functional, cysteine-less variant of Ire1 and performed systematic cysteine cross-linking experiments in native membranes to establish its transmembrane architecture in signaling-active clusters. We show that the transmembrane helices of two neighboring Ire1 molecules adopt an X-shaped configuration independent of the primary cause for ER stress. This suggests that different forms of stress converge in a common, signaling-active transmembrane architecture of Ire1.
Simulations of conformational changes and enzyme-substrate interactions in protein drug targets
(2022)
Finding new drugs is a difficult, time-consuming, and costly challenge, with only a small success rate along the drug discovery pipeline of far less than 10%. The high failure rate of drug discovery projects motivates the integration of computational tools throughout the whole drug discovery pipeline, from target identification to clinical trials. Target identification is the first step in the process. A biological target, e.g., a protein that plays a role in disease, is identified and its molecular mechanism in the disease is studied. Further, a potential binding site on the target, where therapeutic molecules can bind and modulate the target’s activity, needs to be characterized. Computational tools can contribute to improving the initial molecular target elucidation and assessment.
In this thesis, I use computational, physics-based approaches to characterize binding sites of drug targets and to decipher enzyme-substrate interactions, which play a role in disease mechanisms. Molecular dynamics (MD) simulations were applied to study the dynamics of molecules in solution at high temporal and spatial resolution. The method generates time-resolved trajectories of the particles in a system of interest by integrating Newton’s equations of motion numerically, starting from a set of coordinates and velocities. In MD simulations, all atoms of a chosen system, including solvent, are represented explicitly. Atomistic simulations are especially well-suited to study detailed interactions that depend on intermolecular interactions, such as hydration effects, hydrogen bonding, hydrophobic interactions, or subtle chemical differences. System properties are inferred from the trajectories, provided that the force fields, describing the interactions between the particles in the system, have a high accuracy. The bonded and non-bonded interactions are parametrized on experimental and quantum chemical data. The purpose of MD simulations can be to gain insight into the behavior of complex biological systems at molecular level, which often cannot be observed in experiments at the same resolution. With recent advances in computer hardware and simulation software, molecular systems of increasing size and simulation length can be investigated.
In the first part of the thesis, I investigated the conformational ensemble of various protein drug targets. Proteins are dynamic biomacromolecules that can have diverse and nearly isoenergetic conformational states. Ligand binding can shift the equilibrium of this conformational ensemble and can uncover binding sites, called cryptic sites. Cryptic sites only emerge upon small molecule binding and are often flat and featureless, and thus not easily recognized in crystal structures without bound ligands. If new binding sites including cryptic sites are detected, they can potentially be exploited for binding to ligands and enable a druggable target. Druggability is the ability of a protein to bind small, drug-like molecules, which is the basis for rational drug design. In this thesis, I used state-of-the-art physics-based, computational approaches to investigate the conformational ensembles of binding sites. In all studied systems, it is known from experiment that a specific group of ligands can induce conformational changes. The aim is to sample the conformational space made accessible upon ligand binding, yet without using the specific ligand structures or details about their interactions. We are interested in sampling the
pocket conformational states and identifying the respective pocket opening mechanism. For some cases, I additionally assessed whether the observed flexibility is a feature of the protein family, or specific to the protein under consideration.
The first studied system is factor VIIa (FVIIa). FVIIa is an essential part of the coagulation cascade and hence a potential drug target for thrombotic diseases. In addition, I investigated various other trypsin-like serine proteases from the same protein family. The binding pocket of trypsin-like serine proteases is called S1 pocket. An X-ray crystal structure solved by our collaborators reveals that a b-sheet structure in the S1 pocket is distorted by a bound ligand. I resolved the conformational change with MD simulations, starting from the unbound protein structure solvated in water and ions. I observed multiple spontaneous transition events. In 7 out of 22 simulations with the b-sheet as starting structure, the S1 pocket eventually rearranged into a distorted loop structure. These transitions occurred spontaneously and were mediated by water molecules probing the backbone hydrogen bonds. The conformational change studied here controls the onset of substrate binding and catalysis. Furthermore, I used metadynamics simulation, an enhanced-sampling method, to estimate the free energy barrier of this conformational change..
Background: Organoids are morphologically heterogeneous three-dimensional cell culture systems and serve as an ideal model for understanding the principles of collective cell behaviour in mammalian organs during development, homeostasis, regeneration, and pathogenesis. To investigate the underlying cell organisation principles of organoids, we imaged hundreds of pancreas and cholangiocarcinoma organoids in parallel using light sheet and bright-field microscopy for up to 7 days.
Results: We quantified organoid behaviour at single-cell (microscale), individual-organoid (mesoscale), and entire-culture (macroscale) levels. At single-cell resolution, we monitored formation, monolayer polarisation, and degeneration and identified diverse behaviours, including lumen expansion and decline (size oscillation), migration, rotation, and multi-organoid fusion. Detailed individual organoid quantifications lead to a mechanical 3D agent-based model. A derived scaling law and simulations support the hypotheses that size oscillations depend on organoid properties and cell division dynamics, which is confirmed by bright-field microscopy analysis of entire cultures.
Conclusion: Our multiscale analysis provides a systematic picture of the diversity of cell organisation in organoids by identifying and quantifying the core regulatory principles of organoid morphogenesis.
In dieser Bachelorarbeit werden verschiedene Methoden zur Bestimmung der Betriebsfrequenz von CH-Kavitäten untersucht. Aufgrund der geometrisch komplexen Form der Beschleunigungsstruktur, können die Eigenfrequenzen nicht mithilfe von analytischen Methoden bestimmt werden. Üblicherweise werden die Eigenfrequenzen, ihre Ladungsund Stromdichten, sowie die elektromagnetischen Felder über numerische Methoden der Computational Electromagnetics (CEM) ermittelt. Die CEM ist eine junge Disziplin, deren Performanz und Anwendungsgebiete in den letzten 20 Jahren rapide gewachsen sind. Hauptverantwortlich hierfür ist zum einen das exponentielle Wachstum der Rechenleitung bei gleichbleibenden Kosten, zum anderen die Entwicklung und Verbesserung der Algorithmen. Bis zum Ende des letzten Jahrhunderts wurden elektronische Komponenten hauptsächlich dadurch entwickelt, indem Prototypen angefertigt und analysiert wurden. Diese zeitaufwendige und kostspielige Herangehensweise ist heutzutage nahezu vollständig durch CEM-Simulationen ersetzt worden. Die Hauptmethoden der CEM sind die Finite-Differenzen-Methode (FDM), die Momenten-Methode (MoM) und die Finite-Elemente-Methode (FEM). Für die Bestimmung der Eigenwerte und Eigenvektoren der Beschleunigungsstrukturen eignet sich aufgrund der Stabilität von diesen Dreien am besten die Methode der finiten Elemente. Da die FEM ein rechen- und speicherintensives Verfahren ist, wurde in dieser Arbeit nach einer schnelleren Methode gesucht, um die Betriebsfrequenz von CH-Kavitäten zu bestimmen. Hierfür wurden 84 CH-Kavitäten mithilfe von CST Studio Suite erstellt und simuliert. Es handelt sich hierbei um vier Grundtypen, drei wurden bei einer fixierten Sollfrequenz von 300 MHz konstruiert; die Sollfrequenz des vierten Grundtyps betrug 175 MHz. Die Teilchengeschwindigkeit wurde jeweils in 0,01er-Schrtitten von 0,05 c bis 0,25 c variiert. Aus den Untersuchungen der EM-Felder wurde anschließend ein semi-analytisches Modell entwickelt, das aufgrund der Geometrie der CH-Kavität die Betriebsfrequenz liefern soll.
We experimentally investigated the quasifree mechanism (QFM) in one-photon double ionization of He and H2 at 800 eV photon energy and circular polarization with a COLTRIMS reaction microscope. Our work provides new insight into this elusive photoionization mechanism that was predicted by Miron Amusia more than four decades ago. We found the distinct four-fold symmetry in the angular emission pattern of QFM electrons from H2 double ionization that has previously only been observed for He. Furthermore, we provide experimental evidence that the photon momentum is not imparted onto the center of mass in quasifree photoionization, which is in contrast to the situation in single ionization and in double ionization mediated by the shake-off and knock-out mechanisms. This finding is substantiated by numerical results obtained by solving the system’s full-dimensional time-dependent Schrödinger equation beyond the dipole approximation.
This thesis has two main parts.
The first part is based on our publication [1], where we use perturbation theory to calculate decay rates of magnons in the Kitaev-Heisenberg-Γ (KHΓ) model. This model describes the magnetic properties of the material α-RuCl 3 , which is a candidate for a Kitaev spin liquid. Our motivation is to validate a previous calculation from Ref. [2]. In this thesis, we map out the classical phase diagram of the KHΓ model. We use the Holstein-Primakoff
transformation and the 1/S expansion to describe the low temperature dynamics of the Kitaev-Heisenberg-Γ model in the experimentally relevant zigzag phase by spin waves. By parametrizing the spin waves in terms of hermitian fields, we find a special parameter region within the KHΓ model where the analytical expressions simplify. This enables us to construct the Bogoliubov transformation analytically. For a representative point in the special parameter region, we use these results to numerically calculate the magnon damping, which is to leading order caused by the decay of single magnons into two. We also calculate the dynamical structure factor of the magnons.
The second part of this thesis is based on our publication [3], where we use the functional renormalization group to analyze a discontinuous quantum phase transition towards a non-Fermi liquid phase in the Sachdev-Ye-Kitaev (SYK) model. In this thesis, we perform a disorder average over the random interactions in the SYK model. We argue that in the thermodynamic limit, the average renormalization group (RG) flow of the SYK model is identical to the RG flow of an effective disorder averaged model. Using the functional RG, we find a fixed point describing the discontinuous phase transition to the non-Fermi liquid phase at zero temperature. Surprisingly, we find a finite anomalous dimension of the fermions, which indicates critical fluctuations and is unusual for a discontinuous transition. We also determine the RG flow at zero temperature, and relate it to the phase diagram known from the literature.
The direct study of transcription or DNA–protein-binding events, requires imaging of individual genes at molecular resolution. Electron microscopy (EM) can show local detail of the genome. However, direct visualization and analysis of specific individual genes is currently not feasible as they cannot be unambiguously localized in the crowded, landmark-free environment of the nucleus. Here, we present a method for the genomic insertion of gene clusters that can be localized and imaged together with their associated protein complexes in the EM. The method uses CRISPR/Cas9 technology to incorporate several genes of interest near the 35S rRNA gene, which is a frequently occurring, easy-to-identify genomic locus within the nucleolus that can be used as a landmark in micrographs. As a proof of principle, we demonstrate the incorporation of the locus-native gene RDN5 and the locus-foreign gene HSX1. This led to a greater than 7-fold enrichment of RNA polymerase III (Pol III) complexes associated with the genes within the field of view, allowing for a significant increase in the analysis yield. This method thereby allows for the insertion and direct visualization of gene clusters for a range of analyses, such as changes in gene activity upon alteration of cellular or external factors.
Ion stopping in warm dense matter is a process of fundamental importance for the understanding of the properties of dense plasmas, the realization and the interpretation of experiments involving ion-beam-heated warm dense matter samples, and for inertial confinement fusion research. The theoretical description of the ion stopping power in warm dense matter is difficult notably due to electron coupling and degeneracy, and measurements are still largely missing. In particular, the low-velocity stopping range, that features the largest modelling uncertainties, remains virtually unexplored. Here, we report proton energy-loss measurements in warm dense plasma at unprecedented low projectile velocities. Our energy-loss data, combined with a precise target characterization based on plasma-emission measurements using two independent spectroscopy diagnostics, demonstrate a significant deviation of the stopping power from classical models in this regime. In particular, we show that our results are in closest agreement with recent first-principles simulations based on time-dependent density functional theory.
The stellar nucleosynthesis of elements heavier than iron can primarily be attributed to neutron capture reactions in the s and r process. While the s process is considered to be well understood with regards to the stellar sites, phases and conditions where it occurs, nucleosynthesis networks still need accurate neutron capture cross sections
with low uncertainties as input parameters. Their quantitative outputs for the isotopic abundances produced in the s process, coupled with the observable solar abundances, can be used to indirectly infer the expected r process abundances. The two stable gallium isotopes, 69Ga and 71Ga, have been shown in sensitivity studies to have considerable impact on the weak s process in massive stars. The available experimental data, mostly derived from neutron activation measurements for quasi-stellar neutron spectra at kBT = 25 keV, show disagreements up to a factor of three.
Determining the differential neutron capture cross section can provide input data for the whole range of astrophysically relevant energies. To that end, a neutron time of flight experimental campaign at the n_TOF facility at CERN was performed for three months, using isotopically enriched samples of both isotopes. The data taken at the EAR1 experimental area covered a wide neutron energy range from thermal to several hundred keV. The respective differential and spectrum averaged neutron capture cross sections for 69Ga and 71Ga were determined in this thesis. They show good agreement with the evaluated cross sections for 71Ga, but reproduce the deviations from the evaluated data that other, more recent activation measurements showed for 69Ga.
Correlations between the harmonic flow coefficients v1, v2, v3 and v4 of nucleons in semi-peripheral Au+Au collisions at a beam energy of 1.23 AGeV are investigated within the hadronic transport approach ultra-relativistic quantum molecular dynamics (UrQMD). In contrast to ultra-relativistic collision energies (where the flow coefficients are evaluated with respect to the respective event plane), we predict strong correlations between the flow harmonics with respect to the reaction plane. Based on an event-by-event selection of the midrapidity final state elliptic flow of nucleons we show that as a function of rapidity, (I) the sign of the triangular flow changes, (II) that the shape of v4 changes from convex to concave, and (III) that v3∝v1v2 and v4∝v22 for all different event classes, indicating strong correlations between all investigated harmonic flow coefficients.
Die vorliegende Arbeit präsentiert Forschungsarbeiten basierend auf nanoskopischen Oberflächenmessungen an plasmonischen Metaoberflächen und zweidimensionalen Materialien, insbesondere dem halbleitenden Übergangsmetal-Dichalcogenid (TMDC) WS_2. Die Thesis ist in sieben Kapitel untergegliedert. Die Einleitung vermittelt einen Überblick über die treibenden Kräfte hinter der Forschung im Bereich der Nanophotonik an zweidimensionalen Materialsystemen. Die Untersuchung der Licht-Materie-Wechselwirkung an dünnen Materialgrenzflächen zieht sich als roter Faden durch die gesamte Arbeit.
Das zweite Kapitel beschreibt den experimentellen Aufbau, der für die Durchführung der nanoskopischen Messungen in dieser Arbeit implementiert wurde. Es werden theoretische Grundlagen, das Messprinzip und die Implementierung des optischen Rasternahfeldmikroskops (s-SNOM) skizziert. Außerdem wird ein Strom-Spannungs-Rasterkraftmikroskop (c-AFM) im Kontaktmodus genutzt, um elektrische Ströme auf mikroskopischen zweidimensionalen TMDC-Terrassen zu messen. In den darauffolgenden vier Kapiteln werden die Beiträge dieser Arbeit zur Untersuchung der Licht-Materie-Wechselwirkung auf der Nanoskala aus verschiedenen Perspektiven vorgestellt. Jedes Kapitel enthält eine kurze Einleitung, einen Theorieteil, Messdaten oder Simulationsergebnisse sowie eine Analyse; vervollständigt durch einen Schlussteil.
Die zentrale Arbeit an einer metallischen Metaoberfläche aus elliptischen Goldscheiben wird in Kapitel 3 vorgestellt. Der zugehörige Theorieteil führt in das Konzept von Oberflächen-Plasmon-Polaritonen (SPP) ein, das für den Forschungsbereich der Plasmonik im Allgemeinen wesentlich ist. Verschiedene Methoden zur Berechnung der Dispersionsrelation dieser Oberflächenmoden an ein- und mehrschichtigen Grenzflächen werden auf die untersuchte Metaoberflächenprobe angewendet. Das Modell sagt drei verschiedene Moden voraus, die sich an der Grenzfläche ausbreiten. Eine teil-gebundene ins Substrat abstrahlende Oberflächenmode sowie zwei vergrabene stark gebundene anisotrope Moden. Eine auf der Probe platzierte Nanokugel aus Silizium wird als radiale Anregungsquelle verwendet.
Der Vergleich mit s-SNOM-Nahfeldbildern zeigt, dass nur die schwach gebundene geführte Modenresonanz ausreichend angeregt wurde, um durch s-SNOM-Bildgebung nachgewiesen werden zu können. Die schwache Oberflächenbindung erklärt die scheinbar isotrope Ausbreitung auf der anisotropen Oberfläche. Die Beobachtung der verbleibenden stark eingegrenzten anisotropen vergrabenen Moden würde eine verbesserte tiefenempfindliche Auflösung des Systems erfordern, die im Prinzip für Schichtdicken von 20 nm möglich sein sollte. Darüber hinaus wirft die Beobachtung die Frage auf, ob die durch Impuls- und Modenvolumenanpassung der Nanokugel gegebene Anregungseffizienz einen ausreichenden Anregungsquerschnitt erzeugt, um nachweisbare vergrabene SPP-Moden zu erzeugen.
In Kapitel 4 wird die Idee der Visualisierung vergrabener elektrischer Felder mit s-SNOM fortgesetzt. Hier wird es auf die Untersuchung von WS_2 angewendet, einem zweidimensionalen TMDC-Material, welches Photolumineszenz zeigt. Durch die Strukturierung des Galliumphosphid-Substrats unter der hängenden Monolage, die von einer dünnen Schicht aus hBN getragen wird, wird die Photolumineszenzausbeute um den Faktor 10 erhöht. Dies wird durch den Entwurf einer lateralen DBR-Mikrokavität mit zusätzlich optimierter vertikaler Tiefe erreicht, die in das Substrat geätzt wurde.
Die hochauflösende Abbildung der elektrischen Feldverteilung im Resonator wird durch den Einsatz von s-SNOM ermöglicht, um die Verbesserung der Einkopplung durch diese beiden Ansätze zu bewerten. Es konnte festgestellt werden, dass die laterale Struktur überwiegend zur verstärkten Photolumineszenzausbeute beiträgt, während für die Einkopplung keine offensichtliche Verstärkung auf die vertikale Strukturoptimierung zurückgeführt werden konnte.
Das zweidimensionale Material WS_2 wird in Kapitel 5 erneut mit Hilfe von c-AFM untersucht. Unterschiedlich dicke Multilagen auf Graphen und Gold dienen als Tunnelbarrieren für vertikale Ströme zwischen Substrat und leitender c-AFM-Messpitze. Die Daten können mit einem Fowler-Nordheim-Modell mit Parametern für die Tunnelbreite und Schottky-Barrierenhöhen der beiden Grenzflächen erklärt werden. Die Messungen zeigen jedoch eine schwache Reproduzierbarkeit, was eine detailliertere Zusammenfassung der relevanten Fehlerquellen erfordert. In der Schlussfolgerung des Kapitels werden mehrere Schlüsselaspekte vorgeschlagen, die bei künftigen Messungen berücksichtigt werden sollten. Entscheidend ist, dass c-AFM sehr empfindlich auf die Adsorption von Wasserfilmen an der Probenoberfläche reagiert, worunter WS_2-Oberflächen unter Umgebungsbedingungen leiden...
The Born cross sections of the e+e− → D*+D*− and e+e− → D*+D− processes are measured using e+e− collision data collected with the BESIII experiment at center-of-mass energies from 4.085 to 4.600 GeV, corresponding to an integrated luminosity of 15.7 fb−1. The results are consistent with and more precise than the previous measurements by the Belle, Babar and CLEO collaborations. The measurements are essential for understanding the nature of vector charmonium and charmonium-like states.
Neural networks have been recently proposed as variational wave functions for quantum many-body systems [G. Carleo and M. Troyer, Science 355, 602 (2017)]. In this work, we focus on a specific architecture, known as Restricted Boltzmann Machine (RBM), and analyse its accuracy for the spin-1/2 J1−J2 antiferromagnetic Heisenberg model in one spatial dimension. The ground state of this model has a non-trivial sign structure, especially for J2/J1>0.5, forcing us to work with complex-valued RBMs. Two variational Ans\"atze are discussed: one defined through a fully complex RBM, and one in which two different real-valued networks are used to approximate modulus and phase of the wave function. In both cases, translational invariance is imposed by considering linear combinations of RBMs, giving access also to the lowest-energy excitations at fixed momentum k. We perform a systematic study on small clusters to evaluate the accuracy of these wave functions in comparison to exact results, providing evidence for the supremacy of the fully complex RBM. Our calculations show that this kind of Ans\"atze is very flexible and describes both gapless and gapped ground states, also capturing the incommensurate spin-spin correlations and low-energy spectrum for J2/J1>0.5. The RBM results are also compared to the ones obtained with Gutzwiller-projected fermionic states, often employed to describe quantum spin models [F. Ferrari, A. Parola, S. Sorella and F. Becca, Phys. Rev. B 97, 235103 (2018)]. Contrary to the latter class of variational states, the fully-connected structure of RBMs hampers the transferability of the wave function from small to large clusters, implying an increase of the computational cost with the system size.
The 124Xe(p,γ) reaction has been measured for the first time at energies around the Gamow window by using stored ions at the ESR facility. The desired beam energies below 10 MeV/u introduce new experimental challenges like windowless ions detection under UHV conditions, extremely short beam lifetimes and efficient beam deceleration and cooling, all of which have been successfully met.
The neutron-unbound isotope 13Be has been studied in several experiments using different reactions, different projectile energies, and different experimental setups. There is, however, no real consensus in the interpretation of the data, in particular concerning the structure of the low-lying excited states. Gathering new experimental information, which may reveal the 13Be structure, is a challenge, particularly in light of its bridging role between 12Be, where the N = 8 neutron shell breaks down, and the Borromean halo nucleus 14Be. The purpose of the present study is to investigate the role of bound excited states in the reaction product 12Be after proton knockout from 14B, by measuring coincidences between 12Be, neutrons, and γ rays originating from de-excitation of states fed by neutron decay of 13Be. The 13Be isotopes were produced in proton knockout from a 400 MeV/nucleon 14B beam impinging on a CH2 target. The 12 Be-n relative-energy spectrum d σ /d Ef n was obtained from coincidences between 12Be(g.s.) and a neutron, and also as threefold coincidences by adding γ rays, from the de-excitation of excited states in 12Be. Neutron decay from the first 5/2+ state in 13Be to the 2+ state in 12Be at 2.11 MeV is confirmed. An energy independence of the proton-knockout mechanism is found from a comparison with data taken with a 35 MeV/nucleon 14B beam. A low-lying p-wave resonance in 13Be(1/2−) is confirmed by comparing proton- and neutron-knockout data from 14B and 14Be.
Accurate spectroscopy of highly-charged high-Z ions in a storage ring is demonstrated to be feasible by the use of specially adapted crystal optics. The method has been applied for the measurement of the 1s Lamb shift in hydrogen-like gold (Au+78) in a storage ring through spectroscopy of the Lyman x-rays. This measurement represents the first result obtained for a high-Z element using high-resolution wavelength-dispersive spectroscopy in the hard x-ray regime, paving the way for sensitivity to higher- order QED effects.
Attempts to extract the order of the chiral transition of QCD at zero chemical potential, with two dynamical flavors of massless quarks, from simulations with progressively decreasing pion mass, have remained inconclusive because of their increasing numerical cost. In an alternative approach to this problem, we consider the path integral as a function of continuous number Nf of degenerate quarks. If the transition in the chiral limit is first order for Nf≥3, a second-order transition for Nf=2 then requires a tricritical point in between. This, in turn, implies tricritical scaling of the critical boundary line between the first-order and crossover regions as the chiral limit is approached. Noninteger numbers of fermion flavors are easily implemented within the staggered fermion discretization. Exploratory simulations at μ=0 and Nf=2.8, 2.6, 2.4, 2.2, 2.1, on coarse Nτ=4 lattices, indeed show a smooth variation of the critical mass mapping out a critical line in the (m, Nf) plane. For the smallest masses, the line appears consistent with tricritical scaling, allowing for an extrapolation to the chiral limit.
Quasifree one-proton knockout reactions have been employed in inverse kinematics for a systematic study of the structure of stable and exotic oxygen isotopes at the R3B/LAND setup with incident beam energies in the range of 300–450 MeV/u. The oxygen isotopic chain offers a large variation of separation energies that allows for a quantitative understanding of single-particle strength with changing isospin asymmetry. Quasifree knockout reactions provide a complementary approach to intermediate-energy one-nucleon removal reactions. Inclusive cross sections for quasifree knockout reactions of the type AO(p,2p)A−1N have been determined and compared to calculations based on the eikonal reaction theory. The reduction factors for the single-particle strength with respect to the independent-particle model were obtained and compared to state-of-the-art ab initio predictions. The results do not show any significant dependence on proton-neutron asymmetry.
The radiative electron capture (REC) into the K shell of bare Xe ions colliding with a hydrogen gas target has been investigated. In this study, the degree of linear polarization of the K-REC radiation was measured and compared with rigorous relativistic calculations as well as with the previous results recorded for U92+. Owing to the improved detector technology, a significant gain in precision of the present polarization measurement is achieved compared to the previously published results. The obtained data confirms that for medium-Z ions such as Xe, the REC process is a source of highly polarized x rays which can easily be tuned with respect to the degree of linear polarization and the photon energy. We argue, in particular, that for relatively low energies the photons emitted under large angles are almost fully linear polarized.
The proton-removal mechanism of the 12CB reaction induced on a carbon target via elementary nucleon-nucleon scattering is investigated in exclusive triple-coincidence measurements. The observed two-nucleon angular correlations are found to be consistent with quasi-free scattering of a projectile-like proton off a target-like nucleon. Exclusive cross sections for one-step pp and pn interactions are determined as [formula] and [formula], respectively. The extracted quasi-free component amounts up to 58(4)% of the total proton-removal cross section. The results are compared to total proton-removal cross sections obtained from the experiment and eikonal reaction theory.
The Coulomb Dissociation (CD) cross sections of the stable isotopes 92,94,100Mo and of the unstable isotope 93Mo were measured at the LAND/R3B setup at GSI Helmholtzzentrum für Schwerionenforschung in Darmstadt, Germany. Experimental data on these isotopes may help to explain the problem of the underproduction of 92,94Mo and 96,98Ru in the models of p-process nucleosynthesis. The CD cross sections obtained for the stable Mo isotopes are in good agreement with experiments performed with real photons, thus validating the method of Coulomb Dissociation. The result for the reaction 93Mo(γ,n) is especially important since the corresponding cross section has not been measured before. A preliminary integral Coulomb Dissociation cross section of the 94Mo(γ,n) reaction is presented. Further analysis will complete the experimental database for the (γ,n) production chain of the p-isotopes of molybdenum.
The process of electron-loss to the continuum (ELC) has been studied for the collision systems U28++H2 at a collision energy of 50 MeV/u, U28++N2 at 30 MeV/u, and U28++Xe at 50 MeV/u. The energy distributions of cusp electrons emitted at an angle of 0∘ with respect to the projectile beam were measured using a magnetic forward-angle electron spectrometer. For these collision systems far from equilibrium charge state, a significantly asymmetric cusp shape is observed. The experimental results are compared to calculations based on first-order perturbation theory, which predict an almost symmetric cusp shape. Some possible reasons for this discrepancy are discussed.
The decay properties of the Pygmy Dipole Resonance (PDR) have been investigated in the semi-magic N=82 nucleus 140Ce using a novel combination of nuclear resonance fluorescence and γ–γ coincidence techniques. Branching ratios for transitions to low-lying excited states are determined in a direct and model-independent way both for individual excited states and for excitation energy intervals. Comparison of the experimental results to microscopic calculations in the quasi-particle phonon model exhibits an excellent agreement, supporting the observation that the Pygmy Dipole Resonance couples to the ground state as well as to low-lying excited states. A 10% mixing of the PDR and the [21+ x PDR] is extracted.
Navigating a complex environment is assumed to require stable cortical representations of environmental stimuli. Previous experimental studies, however, show substantial ongoing remodeling at the level of synaptic connections, even under behaviorally and environmentally stable conditions. It remains unclear, how these changes affect sensory representations on the level of neuronal populations during basal conditions and how learning influences these dynamics.
Our approach is a joint effort between the analysis of experimental data and theory. We analyze chronic neuronal population activity data – acquired by out collaborators in Mainz – to describe population activity dynamics during basal dynamics and during learning (fear conditioning). The data analysis is complemented by the analysis of a circuit model investigating the link between a neural network’s activity and changes in its underlying structure.
Using chronic two-photon imaging data recorded in awake mouse auditory cortex, we reproduce previous findings that responses of neuronal populations to short complex sounds typically cluster into a near discrete set of possible responses. This means that different stimuli evoke basically the same response and are thus grouped together into one of a small set of possible response modes. The near discrete set of response modes can be utilized as a sensitive and robust means to detect and track changes in population activity over time. Doing so we find that sound representations are subject to a significant ongoing remodeling across the time span of days under basal conditions. Auditory cued fear conditioning introduces a bias into these ongoing dynamics, resulting in a differential generalization both on the level of neuronal populations and on the behavioral level. This means that sounds that are perceived similar to the conditioned stimulus (CS+) show an increased co-mapping to the same response mode the CS+ is mapped to. This differential generalization is also observed in animal behavior, where sounds similar to the CS+ result in the same freezing behavior as the CS+, whereas dissimilar sounds do not. These observations could provide a potential mechanism of stimulus generalization, which is one of the most common phenomena associated with post-traumatic stress disorder, on the level of neuronal populations.
To investigate how the aforementioned changes in neuronal population activity are linked to changes in the underlying synaptic connectivity, we devised a circuit model of excitatory and inhibitory neurons. We studied this firing rate model to investigate the effect of gradual changes in the network’s connectivity on its activity. Apart from an input dominated uni-stable regime (one response per stimulus independent of the network) and a network dominated uni-stable regime (one response per network independent of the stimulus), we also find a multi-stable regime for strong recurrent connectivity and a high ratio of inhibition to excitation. In this regime the model reproduces properties of neural population activity in mouse auditory cortex, including sparse activity, a broad distribution of firing rates, and clustering of stimuli into a near discrete set of response modes. This clustering in the multi-stable regime means that, not only can identical stimuli evoke different responses, depending on the network’s initial condition, but different stimuli can also evoke the same response.
Applying gradual drift to the network connectivity we find periods of stable responses, interrupted by abrupt transitions altering the stimulus response mapping. We study the mechanism underlying these transitions by analyzing changes in the fixed points of this network model, employing a method to numerically find all the fixed points of the system. We find that such abrupt transitions typically cannot be explained by the mere displacement of existing fixed points, but involve qualitative changes in the fixed point structure in the vicinity of the response trajectory. We conclude that gradual synaptic drift can lead to abrupt transitions in stimulus responses and that qualitative changes in the network’s fixed point topology underlie such transitions.
In summary we find that cortical networks display ongoing representational drift under basal conditions that is biased towards a differential generalization during fear conditioning. A circuit model is able to reproduce key characteristics of auditory cortex, including a clustering of stimulus responses into a near discrete set of response modes. Implementing synaptic drift into this model leads to periods of stable responses interrupted by abrupt transitions towards new responses.
Impact parameter sensitive study of inner-shell atomic processes in Xe54+, Xe52+ → Xe collisions
(2020)
In this work, we present a pilot experiment in the experimental storage ring (ESR) at GSI devoted to impact parameter sensitive studies of inner shell atomic processes for bare and He-like xenon ions (Xe54+, Xe52+) colliding with neutral xenon gas atoms. The projectile and target x-rays have been measured at different observation angles for all impact parameters as well as for the impact parameter range of ∼35 - 70 fm.
NeuLAND (New Large-Area Neutron Detector) is the next-generation neutron detector for the R3B (Reactions with Relativistic Radioactive Beams) experiment at FAIR (Facility for Antiproton and Ion Research). NeuLAND detects neutrons with energies from 100 to 1000 MeV, featuring a high detection efficiency, a high spatial and time resolution, and a large multi-neutron reconstruction efficiency. This is achieved by a highly granular design of organic scintillators: 3000 individual submodules with a size of 5 × 5 × 250 cm3 are arranged in 30 double planes with 100 submodules each, providing an active area of 250 × 250 cm2 and a total depth of 3 m. The spatial resolution due to the granularity together with a time resolution of 150 ps ensures high-resolution capabilities. In conjunction with calorimetric properties, a multi-neutron reconstruction efficiency of 50% to 70% for four-neutron events will be achieved, depending on both the emission scenario and the boundary conditions allowed for the reconstruction method. We present in this paper the final design of the detector as well as results from test measurements and simulations on which this design is based.
We discuss the electron-optical properties of a toroidal magnetic sector spectrometer and its suitablilty for electron-positron pair spectroscopy in relativistic ion-atom collisions in the future HESR storage ring at FAIR. With the simultaneous mapping of electrons and positrons and geometric invariants in the lepton trajectorties this instrument offers a very high efficiency for studies of vector momentum correlation in free-free pair production.
New experimental data for dissociation of relativistic 17Ne projectiles incident on targets of lead, carbon, and polyethylene targets at GSI are presented. Special attention is paid to the excitation and decay of narrow resonant states in 17Ne. Distributions of internal energy in the three-body system have been determined together with angular and partial-energy correlations between the decay products in different energy regions. The analysis was done using existing experimental data on 17Ne and its mirror nucleus 17N. The isobaric multiplet mass equation is used for assignment of observed resonances and their spins and parities. A combination of data from the heavy and light targets yielded cross sections and transition probabilities for the Coulomb excitations of the narrow resonant states. The resulting transition probabilities provide information relevant for a better understanding of the 17Ne structure.
Cross sections for neutron-induced reactions of short-lived nuclei are essential for nuclear astrophysics since these reactions in the stars are responsible for the production of most heavy elements in the universe. These reactions are also key in applied domains like energy production and medicine. Nevertheless, neutron-induced cross-section measurements can be extremely challenging or even impossible to perform due to the radioactivity of the targets involved. Indirect measurements through the surrogate-reaction method can help to overcome these difficulties.
The surrogate-reaction method relies on the use of an alternative reaction that will lead to the formation of the same excited nucleus as in the neutron-induced reaction of interest. The decay probabilities (for fission, neutron and gamma-ray emission) of the nucleus produced via the surrogate reaction allow one to constrain models and the prediction of the desired neutron cross sections.
We propose to perform surrogate reaction measurements in inverse kinematics at heavy-ion storage rings, in particular at the CRYRING@ESR of the GSI/FAIR facility. We present the conceptual idea of the most promising setup to measure for the first time simultaneously the fission, neutron and gamma-ray emission probabilities. The results of the first simulations considering the 238U(d,d') reaction are shown, as well as new technical developments that are being carried out towards this set-up.
Quasi-universal behavior of the threshold mass in unequal-mass, spinning binary neutron star mergers
(2021)
The lifetime of the remnant produced by the merger of two neutron stars can provide a wealth of information on the equation of state of nuclear matter and on the processes leading to the electromagnetic counterpart. Hence, it is essential to determine when this lifetime is the shortest, corresponding to when the remnant has a mass equal to the threshold mass, Mth, to prompt collapse to a black hole. We report on the results of more than 360 simulations of merging neutron-star binaries covering 40 different configurations differing in mass ratio and spin of the primary. Using this data, we have derived a quasi-universal relation for Mth and expressed its dependence on the mass ratio and spin of the binary. The new expression recovers the results of Koeppel et al. for equal-mass, irrotational binaries and reveals that Mth can increase (decrease) by 5% (10%) for binaries that have spins aligned (antialigned) with the orbital angular momentum and provides evidence for a nonmonotonic dependence of Mth on the mass asymmetry in the system. Finally, we extend to unequal masses and spinning binaries the lower limits that can be set on the stellar radii once a neutron star binary is detected, illustrating how the merger of an unequal-mass, rapidly spinning binary can significantly constrain the allowed values of the stellar radii.
Neutron-induced cross sections of short-lived nuclei are highly relevant in many domains such as fundamental nuclear physics, astrophysics and applications in nuclear technology. In particular, these cross sections are essential for understanding the synthesis of elements via the s- and r stellar processes. However, the measurement of such cross sections with current techniques is very difficult or even impossible, because of the difficulties to produce and handle the necessary amounts of radioactive nuclei. Reaching the nuclei of interest is only possible by inverting the reaction kinematics with radioactive beams.
In this contribution we present a project for indirectly determining neutron cross sections via the surrogate-reaction method. This project is based on the measurement of transfer- or inelastic-scattering-induced decay probabilities in inverse kinematics at storage rings. The measured probabilities are then used to tune nuclear-reaction models that will provide much more accurate predictions of the desired neutron cross sections. We also discuss a very ambitious, long-term project to directly measure neutron cross sections in inverse kinematics. It consists in the combination of a radioactive beam facility, an ion storage ring and a spallation neutron source.
Im Zentrum dieser Arbeit steht die Diagnostik eines Wasserstoff-Theta-Pinch-Plasmas hinsichtlich der integrierten Elektronen- und Neutralgasdichte mittels Zweifarben Interferometrie. Die integrierte Elektronen- und Neutralgasdichte sind essenzielle Größen, aus welchen sich die Ratenkoeffizienten der Ionisation und Rekombination bei einer Plasma-Ionenstrahl-Wechselwirkung bestimmen lassen.
Ein Theta-Pinch-Plasma ist ein induktiv gezündetes Plasma, wobei das zur Zündung notwendige elektrische Feld durch ein magnetisches Wechselfeld generiert wird. Das induzierte, azimutale elektrische Feld beschleunigt freie Elektronen im Arbeitsgas, welches durch Stoßionisation in den Plasmazustand gebracht wird. Der azimutale Plasmastrom erzeugt einen radialen magnetischen Druckgradienten, der das Plasma komprimiert. Da in axialer Richtung keine Kompressionskraft wirkt, weicht das Plasma einer weiteren Kompression aus, wodurch es zu einer axialen Expansion des Plasmas kommt. Die Expansion erzeugt eine Ionisationswelle im kalten Restgas und es wird eine lange, hoch ionisierte Plasmasäule gebildet.
Dieser hochdynamische Prozess ist mit einem Mach-Zehnder-Interferometer bei der Verwendung von zwei verschiedenen Versionen des Theta-Pinchs zeitaufgelöst untersucht worden. Der Unterschied dieser Versionen liegt in der Geometrie und Induktivität der Spulen, wobei zum einen eine zylindrische und zum anderen eine sphärische Spule eingesetzt worden ist. Das grundlegende Messprinzip beruht darauf, dass das Plasma einen Brechungsindex besitzt, welcher von den Dichten der im Plasma enthaltenen Teilchenspezies abhängt. In einem Wasserstoffplasmas sind dies der Beitrag der freien Elektronen und der des Neutralgases, wodurch ein Zweifarben-Interferometer eingesetzt wird. Um eine von den Laserintensitäten unabhängige Messung zu ermöglichen, wird das heterodyne Verfahren benutzt, bei dem die Referenzstrahlen beider Wellenlängen jeweils mit einem akusto-optischen Modulator frequenzverschoben werden. Durch einen Vergleich mit einem stationären Referenzsignal mittels eines I/Q-Demodulators wird die interferometrische Phasenverschiebung aus dem Messsignal extrahiert.
Mit diesem diagnostischen Verfahren ist die integrierte Elektronen- und Neutralgasdichte des Theta-Pinch-Plasmas bei Variation des Arbeitsdrucks und der Ladespannung der Kondensatorbank untersucht worden. Mit der zylindrischen Experimentversion ist eine optimale Kombination aus integrierter Elektronendichte und effektivem Ionisationsgrad η von (1,45 ± 0,04) · 1018 cm−2 bei η = (0,826 ± 0,022) bei einem Arbeitsdruck von 20 Pa und einer Ladespannung von 16 kV ermittelt worden. Dagegen beträgt die optimale Kombination bei einem Arbeitsdruck von 20 Pa und einer Ladespannung von 18 kV bei Verwendung der sphärischen Experimentversion lediglich (1,23 ± 0,03) · 1018 cm−2 bei η = (0,699 ± 0,019).
Des Weiteren ist bei beiden Experimentversionen nachgewiesen worden, dass die integrierte Elektronendichte dem oszillierenden Strom folgend periodische lokale Maxima zeigt, welche zeitlich mit signifikanten Einbrüchen in der integrierten Neutralgasdichte zusammenfallen. Diese Einbrüche werden durch die axiale Expansion des Plasmas und der damit verbundenen Ionisationswelle im Restgas erzeugt. Neben diesem zentralen Teil dieser Arbeit ist eine lasergestützte polarimetrische Diagnostik durchgeführt worden, mit der die longitudinale Komponente der magnetischen Flussdichte der Theta-Pinch-Spulen zeit- und ortsaufgelöst bestimmt worden ist. Als Messprinzip ist der Faraday-Effekt eines magneto-optischen TGGKristalls verwendet worden.
Vor der polarimetrischen Diagnostik ist der TGG-Kristall bezüglich seiner Verdet- Konstante kalibriert worden, wobei ein Wert von V = (−149,7 ± 6,4) rad/Tm gemessen worden ist. Die ortsaufgelöste polarimetrische Diagnostik ist durch einen Seilzug ermöglicht worden, mit dem der TGG-Kristall auf einem Schlitten an unterschiedliche Positionen entlang der Spulenachse gefahren werden konnte. An den jeweiligen Messpunkten ist für beide Experimentversionen die magnetische Flussdichte für verschiedene Ladespannungen zeitaufgelöst bestimmt worden. Als Messverfahren ist dabei das Δ/Σ-Verfahren eingesetzt worden, mit dem sich eine intensitätsunabhängige Messung erzielen ließ.
Die ortsaufgelösten Messergebnisse fallen gegenüber Simulationen allerdings zu niedrig aus. Bei der zylindrischen Spule betragen die Abweichungen im Spulenzentrum circa 14 - 16% und bei der sphärischen Spule in etwa 16 - 18%. Bei einer Normierung der Messwerte und der simulierten Werte auf den jeweiligen Wert im Zentrum ist dagegen innerhalb der Fehler eine völlige Übereinstimmung zwischen den Messwerten und der Simulation für die zylindrische Spule erzielt worden. Als Ursache der negativen Abweichungen wird die Hysterese des TGG-Kristalls diskutiert. Es zeigt sich insbesondere zu Beginn der Entladung eine zeitliche Verzögerung der gemessenen magnetischen Flussdichte gegenüber dem Strom, die in der Umgebung des Stromnulldurchgangs besonders stark ausgeprägt ist.
The aim of this thesis is to provide a complete and consistent derivation of second-order dissipative relativistic spin hydrodynamics from quantum field theory. We will proceed in two main steps. The first one is the formulation of spin kinetic theory from quantum field theory using the Wigner-function formalism and performing an expansion in powers of the Planck constant. The essential ingredient here is the nonlocal collision term. We will find that the nonlocality of the collision term arises at first order in the Planck constant and is responsible for the spin alignment with vorticity, as it allows for conversion between spin and orbital angular momentum.
In the second step, this kinetic theory is used as the starting point to derive hydrodynamics including spin degrees of freedom. The so-called canonical form of the conserved currents follows from Noether’s theorem.
Applying an HW pseudo-gauge transformation, we obtain a spin tensor and energy-momentum tensor with obvious physical interpretation. Promoting all components of the HW tensors to be dynamical, we derive
second-order dissipative spin hydrodynamics. The additional equations of motion for the dissipative currents are obtained from kinetic theory generalizing the method of moments to include spin degrees of freedom.
Die vorgelegte Dissertation behandelt den Einfluss homöostatischer Adaption auf die Informationsverarbeitung und Lenrprozesse in neuronalen Systemen. Der Begriff Homöostase bezeichnet die Fähigkeit eines dynamischen Systems, bestimmte interne Variablen durch Regelmechanismen in einem dynamischen Gleichgewicht zu halten. Ein klassisches Beispiel neuronaler Homöostase ist die dynamische Skalierung synaptischer Gewichte, wodurch die Aktivität bzw. Feuerrate einzelner Neuronen im zeitlichen Mittel konstant bleibt. Bei den von uns betrachteten Modellen handelt es sich um eine duale Form der neuronalen Homöostase. Das bedeutet, dass für jedes Neuron zwei interne Parameter an eine intrinsische Variable wie die bereits erwähnte mittlere Aktivität oder das Membranpotential gekoppelt werden. Eine Besonderheit dieser dualen Adaption ist die Tatsache, dass dadurch nicht nur das zeitliche Mittel einer dynamischen Variable, sondern auch die zeitliche Varianz, also die stärke der Fluktuation um den Mittelwert, kontrolliert werden kann. In dieser Arbeit werden zwei neuronale Systeme betrachtet, in der dieser Aspekt zum Tragen kommt.
Das erste behandelte System ist ein sogennantes Echo State Netzwerk, welches unter die Kategorie der rekurrenten Netzwerke fällt. Rekurrente neuronale Netzwerke haben im Allgemeinen die Eigenschaft, dass eine Population von Neuronen synaptische Verbindungen besitzt, die auf die Population selbst projizieren, also rückkoppeln. Rekurrente Netzwerke können somit als autonome (falls keinerlei zusätzliche externe synaptische Verbindungen existieren) oder nicht-autonome dynamische Systeme betrachtet werden, die durch die genannte Rückkopplung komplexe dynamische Eigenschaften besitzen. Abhängig von der Struktur der rekurrenten synaptischen Verbindungen kann beispielsweise Information aus externem Input über einen längeren Zeitraum gespeichert werden. Ebenso können dynamische Fixpunkte oder auch periodische bzw. chaotische Aktivitätsmuster entstehen. Diese dynamische Vielseitigkeit findet sich auch in den im Gehirn omnipräsenten rekurrenten Netzwerken und dient hier z.B. der Verarbeitung sensorischer Information oder der Ausführung von motorischen Bewegungsmustern. Das von uns betrachtete Echo State Netzwerk zeichnet sich dadurch aus, dass rekurrente synaptische Verbindungen zufällig generiert werden und keiner synaptischen Plastizität unterliegen. Verändert werden im Zuge eines Lernprozesses nur Verbindungen, die von diesem sogenannten dynamischen Reservoir auf Output-Neuronen projizieren. Trotz der Tatsache, dass dies den Lernvorgang stark vereinfacht, ist die Fähigkeit des Reservoirs zur Verarbeitung zeitabhängiger Inputs stark von der statistischen Verteilung abhängig, die für die Generierung der rekurrenten Verbindungen verwendet wird. Insbesondere die Varianz bzw. die Skalierung der Gewichte ist hierbei von großer Bedeutung. Ein Maß für diese Skalierung ist der Spektralradius der rekurrenten Gewichtsmatrix.
In vorangegangenen theoretischen Arbeiten wurde gezeigt, dass für das betrachtete System ein Spektralradius nahe unterhalb des kritischen Wertes von 1 zu einer guten Performance führt. Oberhalb dieses Wertes kommt es im autonomen Fall zu chaotischem dynamischen Verhalten, welches sich negativ auf die Informationsverarbeitung auswirkt. Der von uns eingeführte und als Flow Control bezeichnete duale Adaptionsmechanismus zielt nun darauf ab, über eine Skalierung der synaptischen Gewichte den Spektralradius auf den gewünschten Zielwert zu regulieren. Essentiell ist hierbei, dass die verwendete Adaptionsdynamik im Sinne der biologischen Plausibilität nur auf lokale Größen zurückgreift. Dies geschieht im Falle von Flow Control über eine Regulation der im Membranpotential der Zelle auftretenden Fluktuationen. Bei der Evaluierung der Effektivität von Flow Control zeigte sich, dass der Spektralradius sehr präzise kontrolliert werden kann, falls die Aktivitäten der Neuronen in der rekurrenten Population nur schwach korreliert sind. Korrelationen können beispielsweise durch einen zwischen den Neuronen stark synchronisierten externen Input induziert werden, der sich dementsprechend negativ auf die Präzision des Adaptionsmechanismus auswirkt.
Beim Testen des Netzwerks in einem Lernszenario wirkte sich dieser Effekt aber nicht negativ auf die Performance aus: Die optimale Performance wurde unabhängig von der stärke des korrelierten Inputs für einen Spektralradius erreicht, der leicht unter dem kritischen Wert von 1 lag. Dies führt uns zu der Schlussfolgerung, dass Flow Control unabhängig von der Stärke der externen Stimulation in der Lage ist, rekurrente Netze in einen für die Informationsverarbeitung optimalen Arbeitsbereich einzuregeln.
Bei dem zweiten betrachteten Modell handelt es sich um ein Neuronenmodell mit zwei Kompartimenten, welche der spezifischen Anatomie von Pyramidenneuronen im Kortex nachempfunden ist. Während ein basales Kompartiment synaptischen Input zusammenfasst, der in Dendriten nahe des Zellkerns auftritt, repräsentiert das zweite apikale Kompartiment die im Kortex anzutreffende komplexe dendritische Baumstruktur. In früheren Experimenten konnte gezeigt werden, dass eine zeitlich korrelierte Stimulation sowohl im basalen als auch apikalen Kompartiment eine deutlich höhere neuronale Aktivität hervorrufen kann als durch Stimulation nur einer der beiden Kompartimente möglich ist. In unserem Modell können wir zeigen, dass dieser Effekt der Koinzidenz-Detektion es erlaubt, den Input im apikalen Kompartiment als Lernsignal für synaptische Plastizität im basalen Kompartiment zu nutzen. Duale Homöostase kommt auch hier zum Tragen, da diese in beiden Kompartimenten sicherstellt, dass sich der synaptische Input hinsichtlich des zeitlichen Mittels und der Varianz in einem für den Lernprozess benötigten Bereich befindet. Anhand eines Lernszenarios, das aus einer linearen binären Klassifikation besteht, können wir zeigen, dass sich das beschriebene Framework für biologisch plausibles überwachtes Lernen eignet.
Die beiden betrachteten Modelle zeigen beispielhaft die Relevanz dualer Homöostase im Hinblick auf zwei Aspekte. Das ist zum einen die Regulation rekurrenter neuronaler Netze in einen dynamischen Zustand, der für Informationsverarbeitung optimal ist. Der Effekt der Adaption zeigt sich hier also im Verhalten des Netzwerks als Ganzes. Zum anderen kann duale Homöostase, wie im zweiten Modell gezeigt, auch für Plastizitäts- und Lernprozesse auf der Ebene einzelner Neuronen von Bedeutung sein. Während neuronale Homöostase im klassischen Sinn darauf beschränkt ist, Teile des Systems möglichst präzise auf einen gewünschten Mittelwert zu regulieren, konnten wir Anhand der diskutierten Modelle also darlegen, dass eine Kontrolle des Ausmaßes von Fluktuationen ebenfalls Einfluss auf die Funktionalität neuronaler Systeme haben kann.
In der vorliegenden Arbeit wurde die Ionisation von Atomen und Molekülen in starken Laserfeldern experimentell untersucht. Hierbei kam die COLTRIMS-Technik zur koinzidenten Messung der Impulse aller aus einem Ionisationsereignis stammender Ionen und Elektronen zum Einsatz. Unter Mitwirkung des Autors wurde ein COLTRIMS-Reaktionsmikroskop umgebaut und mit einem neuen Spektrometer sowie einer atomaren Wasserstoffquelle ausgestattet. Des Weiteren entstand ein interferometrischer Aufbau zur Erzeugung von Zwei-Farben-Feldern. Aus jedem der vorgestellten Experimente konnten Informationen über die elektronische Wellenfunktion an der Grenze zum klassisch verbotenen Bereich gewonnen werden. Dies geschah sowohl im Hinblick auf die Amplitude, als auch auf die Phase der Wellenfunktion. Mit dem Wasserstoffatom (Kapitel 9), dem Wasserstoffmolekül (Kapitel 10) und dem Argondimer (Kapitel 11) wurden drei Systeme unterschiedlicher Komplexität gewählt.
Folgend auf den ersten Realisierungen von Bose-Einstein Kondensaten erschienen weitere innovative Experimente, die sich in den optischen Gittern gefangenen Quantengasen widmeten. In diesen zahlreichen, wissenschaftlichen Untersuchungen konnten die Eigenschaften von Bose-Einstein Kondensaten besser verstanden werden. Das Prinzip von Vielteilchensystemen, gefangen in einem periodischen Potential, bot eine Plattform zur Untersuchung weiterer Quantenphasen.
Eine konzeptionell einfache Modifikation von solchen Systemen erhält man durch die Kopplung der Grundzustände der gefangenen Teilchen an hoch angeregten Zuständen mithilfe einer externen Lichtquelle. Im Falle dessen, dass diese Zustände nahe der Ionisationsgrenze des Atoms liegen, spricht man von Rydberg-Zuständen und Atome, welche zu diesen Zuständen angeregt werden, bezeichnet man als Rydberg-Atome. Eines der vielen charakteristischen Eigenschaften von Rydberg-Atomen ist die Fähigkeit über große Entfernungen jenseits der atomaren Längenskalen zu wechselwirken. Im Rahmen von Vielteilchensystemen wurden dementsprechend Kristallstrukturen aus gefangenen Rydberg-Atomen experimentell beobachtet.
Nun stellt sich die Frage, was mit einem gefangenen Bose-Einstein Kondensat passiert, dessen Teilchen an langreichweitig wechselwirkenden Zuständen gekoppelt sind. Gibt es ein Parameterregime, in dem sowohl Kristallstruktur als auch Suprafluidität in solchen Systemen koexistieren können? Dies ist die zentrale Frage dieser Arbeit, die sich mit der Theorie von gefangenen Quantengasen gekoppelt an Rydberg-Zuständen auseinandersetzt.
When a very strong light field is applied to a molecule an electron can be ejected by tunneling. In order to quantify the time-resolved dynamics of this ionization process, the concept of the Wigner time delay can be used. The properties of this process can depend on the tunneling direction relative to the molecular axis. Here, we show experimental and theoretical data on the Wigner time delay for tunnel ionization of H2 molecules and demonstrate its dependence on the emission direction of the electron with respect to the molecular axis. We find, that the observed changes in the Wigner time delay can be quantitatively explained by elongated/shortened travel paths of the emitted electrons, which occur due to spatial shifts of the electrons’ birth positions after tunneling. Our work provides therefore an intuitive perspective towards the Wigner time delay in strong-field ionization.