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
- Physik (724)
- Frankfurt Institute for Advanced Studies (FIAS) (75)
- MPI für Biophysik (55)
- ELEMENTS (17)
- Präsidium (13)
- Biochemie, Chemie und Pharmazie (11)
- Buchmann Institut für Molekulare Lebenswissenschaften (BMLS) (9)
- Biochemie und Chemie (8)
- Biowissenschaften (7)
- Informatik und Mathematik (5)
- Medizin (5)
- Helmholtz International Center for FAIR (3)
- MPI für Hirnforschung (2)
- Zentrum für Biomolekulare Magnetische Resonanz (BMRZ) (2)
- Ernst Strüngmann Institut (1)
- Exzellenzcluster Makromolekulare Komplexe (1)
- Geowissenschaften / Geographie (1)
- Informatik (1)
- Psychologie und Sportwissenschaften (1)
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.
The new heavy ion superconducting continuous wave HElmholtz LInear ACcelerator (HELIAC) is under construction at GSI. A normal conducting injector, comprising an ECR ion source, an RFQ and a DTL, is recently in development. The new Interdigital H-mode DTL, presented in this paper, accelerates the heavy ion beam from 300 to 1400 keV/u, applying an Alternating Phase Focusing (APF) beam dynamics scheme. This APF section, consisting of two separately controlled tanks, has to provide for stable routine operation with assistance of dedicated beam diagnostics devices in the Intertank section. The installed quadrupole lenses and beam steerers installed there ensure full transmission in a wide range of input beam parameters.
The article presents the results of numerical and experimental investigations of guided wave propagation in aluminum plates with variable thickness. The shapes of plate surfaces have been specially designed and manufactured using a CNC milling machine. The shapes of the plates were defined by sinusoidal functions varying in phase shift, which forced the changes in thickness variability alongside the propagation path. The main aim of the study is to analyze the wave propagation characteristics caused by non-uniform thickness. In the first step, the influence of thickness variability on the time course of propagating waves has been analyzed theoretically. The study proves that the wave propagation signals can be determined based on knowledge about the statistical description of the specimen geometry. The histograms of thickness distribution together with the a priori knowledge of the dispersion curves were used to develop an iterative procedure assuming that the signal from the previous step becomes the excitation in the next step. Such an approach allowed for taking into account the complex geometry of the plate and rejecting the assumption about the constant average thickness alongside the propagation path. In consequence, it was possible to predict correctly the signal time course, as well as the time of flight and number of propagating wave modes in specimens with variable thickness. It is demonstrated that theoretical signals predicted in this way coincide well with numerical and experimental results. Moreover, the novel procedure allowed for the correct prediction of the occurrence of higher-order modes.
The appearance of strangeness in the form of hyperons within the inner core of neutron stars is expected to affect its detectable properties, such as its global structure or gravitational wave emission. This work explores the parameter space of hyperonic stars within the framework of the Relativistic Mean Field model allowed by the present uncertainties in the state-of-the-art nuclear and hypernuclear experimental data. We impose multi-physics constraints at different density regimes to restrict the parameter space: Chiral effective field theory, heavy-ion collision data, and multi-messenger astrophysical observations of neutron stars. We investigate possible correlations between empirical nuclear and hypernuclear parameters, particularly the symmetry energy and its slope, with observable properties of neutron stars. We do not find a correlation for the hyperon parameters and the astrophysical data. However, the inclusion of hyperons generates a tension between the astrophysical and heavy-ion data constraining considerably the available parameter space.
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.
Bounded rationality is one crucial component in human behaviours. It plays a key role in the typical collective behaviour of evacuation, in which heterogeneous information can lead to deviations from optimal choices. In this study, we propose a framework of deep learning to extract a key dynamical parameter that drives crowd evacuation behaviour in a cellular automaton (CA) model. On simulation data sets of a replica dynamic CA model, trained deep convolution neural networks (CNNs) can accurately predict dynamics from multiple frames of images. The dynamical parameter could be regarded as a factor describing the optimality of path-choosing decisions in evacuation behaviour. In addition, it should be noted that the performance of this method is robust to incomplete images, in which the information loss caused by cutting images does not hinder the feasibility of the method. Moreover, this framework provides us with a platform to quantitatively measure the optimal strategy in evacuation, and this approach can be extended to other well-designed crowd behaviour experiments.
Highlights
• Sampling the large conformational space of disordered proteins requires extensive molecular dynamics (MD) simulations.
• Fragment assembly complements MD simulations to produce extensive ensembles of disordered proteins with atomic detail.
• Hierarchical chain growth (HCG) ensembles capture key experimental descriptors “out of the box”.
• HCG has revealed local structural characteristics associated with protein dysfunction in neurodegeneration.
Abstract
Disordered proteins and nucleic acids play key roles in cellular function and disease. Here, we review recent advances in the computational exploration of the conformational dynamics of flexible biomolecules. While atomistic molecular dynamics (MD) simulation has seen a lot of improvement in recent years, large-scale computing resources and careful validation are required to simulate full-length disordered biopolymers in solution. As a computationally efficient alternative, hierarchical chain growth (HCG) combines pre-sampled chain fragments in a statistically reproducible manner into ensembles of full-length atomically detailed biomolecular structures. Experimental data can be integrated during and after chain assembly. Applications to the neurodegeneration-linked proteins α-synuclein, tau, and TDP-43, including as condensate, illustrate the use of HCG. We conclude by highlighting the emerging connections to AI-based structural modeling including AlphaFold2.
In this work, the phase diagram of the 2+1-dimensional Gross-Neveu model is investigated with baryon chemical potential as well as chiral chemical potential in the mean-field approximation. We study the theory using two lattice discretizations, which are both based on naive fermions. An inhomogeneous chiral phase is observed only for one of the two discretizations. Our results suggest that this phase disappears in the continuum limit.
This Ph. D. thesis with the title "Characterisation of laser-driven radiation beams: Gamma-ray dosimetry and Monte Carlo simulations of optimised target geometry for record-breaking efficiency of MeV gamma-sources" is dedicated to the study of the acceleration of electrons by intense sub-picosecond laser pulses propagating in a sub-millimeter plasma with near-critical electron density (NCD) and resulting generation of the gamma bremsstrahlung and positrons in the targets of different materials and thickness.
Laser-driven particle acceleration is an area of increasing scientific interest since the recent development of short pulse, high-intensity laser systems. The interaction of intense high-energy, short-pulse lasers with solid targets leads to the production of high-energy electrons in the relativistic laser intensity regime of more than 1018 W /cm2. These electrons play the leading role in the first stage of the interaction of laser with matter, which leads to the creation of laser sources of particles and radiation. Therefore, the optimisation of the electron beam parameters in the direction of increasing the effective temperature and beam charge, together with a slight divergence, plays a decisive role, especially for further detection and characterisation of laser-driven photon and positron beams.
In the context of this work, experiments were carried out at the PHELIX laser system (Petawatt High-Energy Laser for Heavy Ion eXperiments) at GSI Helmholtz Center for Heavy-Ion Research GmbH in Darmstadt, Germany. This thesis presents a thermoluminescence dosimetry (TLD) based method for the measurement of bremsstrahlung spectra in the energy range from 30 keV to 100 MeV. The results of the TLD measurements reinforced the observed tendency towards the strong increase of the mean electron energy and number of super-ponderomotive electrons. In the case of laser interaction with long-scale NCD-plasmas, the dose caused by the gamma-radiation measured in the direction of the laser pulse propagation showed a 1000-fold increase compared to the high contrast shots onto plane foils and doses measured perpendicular to the laser propagation direction for all used combinations of targets and laser parameters.
In this thesis I present novel characterisation method using a combination of TLD measurements and Monte Carlo FLUKA simulations applicable to laser-driven beams. The thermoluminescence detector-based spectrometry method for simultaneous detection of electrons and photons from relativistic laser-induced plasmas initially developed by Behrens et al. (Behrens et al., 2003) and further applied in experiments at PHELIX laser (Horst et al., 2015) delivered good spectral information from keV energies up to some MeV, but as it was presented in (Horst et al., 2015) this method was not really suitable to resolve the content of photon spectra above 10 MeV because of the dominant presence of electrons. Therefore, I created new evaluation method of the incident electron spectra from the readings of TLDs. For this purpose, by means of MatLab programming language an unfolding algorithm was written. It was based on a sequential enumeration of matching data series of the dose values measured by the dosimeters and calculated with of FLUKA-simulations. The significant advantage of this method is the ability to obtain the spectrum of incident electrons in the low energy range from 1 keV, which is very difficult to measure reliably using traditional electron spectrometers.
The results of the evaluation of the effective temperature of super-ponderomotive electrons retrieved from the measured TLD-doses by means of the Monte-Carlo simulations demonstrated, that application of low density polymer foam layers irradiated by the relativistic sub-ps laser pulse provided a strong increase of the electron effective temperature from 1.5 - 2 MeV in the case of the relativistic laser interaction with a metallic foil up to 13 MeV for the laser shots onto the pre-ionized foam and more than 10 times higher charge carried by relativistic electrons.
The progressive simulation method of whole electron spectra described with two -temperatures Maxwellian distribution function has been developed and the results of dose simulations were compared with the acquired experimental data. The advanced feature of this method, which distinguishes it from the results of the simulation of the photon spectrum using the interaction with the target of mono-energetic electron beams (Nilgün Demir, 2013; Nilgün Demir, 2019) or the initial electron spectrum expressed as a function of one electron temperature (Fiorini, 2012), is the ability to simulate the initial electron spectrum described by the Maxwellian distribution function with two temperatures.
The important objective of this thesis was dedicated to the study and characterisation of laser-driven photon beams. In addition to this, the positron beams were evaluated. The investigation of bremsstrahlung photons and positrons spectra from high Z targets by varying the target thickness from 10 µm to 4 mm in simulated models of the interactions of electron spectra with Maxwellian distribution functions allowed to define an optimal thickness when the fluences of photons and positrons are maximal. Furthermore based on the results of FLUKA simulations the gold material was found to be the most suitable for the future experiments as e − γ target because of its highest bremsstrahlung yield.
Additionally Monte Carlo simulations were performed applying the obtained electron beam parameters from the electron acceleration process in laser-plasma interactions simulated with particle-in-cell (PIC) code for two laser energies of 20 J and 200 J. The corresponding electron spectra were imported into a Monte Carlo code FLUKA to simulate the production process of bremsstrahlung photons and positrons in Au converter. FLUKA simulations showed the record conversion of efficiency in MeV gammas can reach 10%, which reinforces the generation of positrons. The obtained results demonstrate the advantages of long-scale plasmas of near critical density (NCD) to increase the parameters of MeV particles and photon beams generated in relativistic laser-plasma interaction. The efficiency of the laser-driven generation of MeV electrons and photons by application of low-density polymer foams is essentially enhanced.
We study the μ-μ45-T phase diagram of the 2+1-dimensional Gross-Neveu model, where μ denotes the ordinary chemical potential, μ45 the chiral chemical potential and T the temperature. We use the mean-field approximation and two different lattice regularizations with naive chiral fermions. An inhomogeneous phase at finite lattice spacing is found for one of the two regularizations. Our results suggest that there is no inhomogeneous phase in the continuum limit. We show that a chiral chemical potential is equivalent to an isospin chemical potential. Thus, all results presented in this work can also be interpreted in the context of isospin imbalance.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √sNN = 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconstructed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[72±63 (stat.)±35 (syst.)] keV are compatible with predictions from effective field theories and conclusively confirm that the 3ΛH is a weakly-bound system.
The most precise measurements to date of the 3ΛH lifetime τ and Λ separation energy BΛ are obtained using the data sample of Pb-Pb collisions at √sNN = 5.02 TeV collected by ALICE at the LHC. The 3ΛH is reconstructed via its charged two-body mesonic decay channel (3ΛH→ 3He + π− and the charge-conjugate process). The measured values τ=[253±11 (stat.)±6 (syst.)] ps and BΛ=[72±63 (stat.)±36 (syst.)] keV are compatible with predictions from effective field theories and conclusively confirm that the 3ΛH is a weakly-bound system.
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.
The gas-phase reaction of O + H₃⁺ has two exothermic product channels: OH+ + H2 and H2O+ + H. In the present study, we analyze experimental data from a merged-beams measurement to derive thermal rate coefficients resolved by product channel for the temperature range from 10 to 1000 K. Published astrochemical models either ignore the second product channel or apply a temperature-independent branching ratio of 70% versus 30% for the formation of OH+ + H2 versus H2O+ + H, respectively, which originates from a single experimental data point measured at 295 K. Our results are consistent with this data point, but show a branching ratio that varies with temperature reaching 58% versus 42% at 10 K. We provide recommended rate coefficients for the two product channels for two cases, one where the initial fine-structure population of the O(3PJ) reactant is in its J = 2 ground state and the other one where it is in thermal equilibrium.
Alternating acquisition of background and sample spectra is often employed in conventional Fourier-transform infrared spectroscopy or ultraviolet–visible spectroscopy for accurate background subtraction. For example, for solvent background correction, typically a spectrum of a cuvette with solvent is measured and subtracted from a spectrum of a cuvette with solvent and solute. Ultrafast spectroscopies, though, come with many peculiarities that make the collection of well-matched, subtractable background and sample spectra challenging. Here, we present a demountable split-sample cell in combination with a modified Lissajous scanner to overcome these challenges. It allows for quasi-simultaneous measurements of background and sample spectra, mitigating the effects of drifts of the setup and maintaining the beam and sample geometry when swapping between background and sample measurements. The cell is moving between subsequent laser shots to refresh the excited sample volume. With less than 45 μl of solution for 150 μm optical thickness, sample usage is economical. Cell assembly is a key step and covered in an illustrated protocol.
Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of ATP. Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt-bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.
Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of ATP. Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt-bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.
Cryo-electron tomography (CryoET) resolves individual macromolecules inside living cells. However, the complex composition and high density of cells challenge the faithful identification of features in tomograms. Here, we capitalize on recent advances in electron tomography and demonstrate that 3D template matching (TM) localizes a wide range of structures inside crowded eukaryotic cells with confidence 10 to 100-fold above the noise level. We establish a TM pipeline with systematically tuned parameters for automated, objective and comprehensive feature identification. High-fidelity and high-confidence localizations of nuclear pore complexes, vaults, ribosomes, proteasomes, lipid membranes and microtubules, and individual subunits, demonstrate that TM is generic. We resolve ~100-kDa proteins, connect the functional states of complexes to their cellular localization, and capture vaults carrying ribosomal cargo in situ. By capturing individual molecular events inside living cells with defined statistical confidence, high-confidence TM greatly speeds up the CryoET workflow and sets the stage for visual proteomics.