Refine
Year of publication
- 2023 (327) (remove)
Document Type
- Preprint (170)
- Article (102)
- Doctoral Thesis (30)
- Conference Proceeding (12)
- Contribution to a Periodical (9)
- Book (3)
- Bachelor Thesis (1)
Has Fulltext
- yes (327)
Is part of the Bibliography
- no (327)
Keywords
Institute
- Physik (327) (remove)
Using a sample of (10.09 ± 0.04) × 109 J/ψ decays collected with the BESIII detector, partial wave analyses of the decay J/ψ → γK0SK0Sπ0 are performed within the K0SK0Sπ0 invariant mass region below 1.6 GeV/c2. The covariant tensor amplitude method is used in both mass independent and mass dependent approaches. Both analysis approaches exhibit dominant pseudoscalar and axial vector components, and show good consistency for the other individual components. Furthermore, the mass dependent analysis reveals that the K0SK0 Sπ0 invariant mass spectrum for the pseudoscalar component can be well described with two isoscalar resonant states using relativistic Breit-Wigner model, i.e., the η(1405) with a mass of 1391.7±0.7+11.3 −0.3 MeV/c 2 and a width of 60.8±1.2+5.5 −12.0 MeV, and the η(1475) with a mass of 1507.6±1.6+15.5−32.2 MeV/c2 and a width of 115.8±2.4 +14.8 −10.9 MeV. The first and second uncertainties are statistical and systematic, respectively. Alternate models for the pseudoscalar component are also tested, but the description of the K0SK0Sπ0invariant mass spectrum deteriorates significantly.
We present our recent results on antiheavy-antiheavy-light-light tetraquark systems using lattice QCD. Our study of the b¯b¯us four-quark system with quantum numbers JP=1+ and the b¯c¯ud four-quark systems with I(JP)=0(0+) and I(JP)=0(1+) utilizes scattering operators at the sink to improve the extraction of the low-lying energy levels. We found a bound state for b¯b¯us with Ebind,b¯b¯us=(−86±22±10)MeV, but no indication for a bound state in both b¯c¯ud channels. Moreover, we show preliminary results for b¯b¯ud with I(JP)=0(1+), where we used scattering operators both at the sink and the source. We found a bound state and determined its infinite-volume binding energy with a scattering analysis, resulting in Ebind,b¯b¯ud=(−103±8)MeV.
This thesis is concerned with the investigation of static and dynamic properties of quantum Heisenberg paramagnets in the absence of a magnetic field and therefore for vanishing magnetization. For this purpose a new formulation of the spin functional renormalization group (SFRG) is employed. The first manifestations of the SFRG were developed by Krieg and Kopietz, motivated by the FRG approach to ordinary field theories and the older works of Vaks, Larkin and Pikin on diagrammatic methods for spin operators.
The main idea is to study quantum spin systems by considering the evolution of correlation functions under a continuous deformation of the interaction between magnetic moments, starting from a solvable limit. This leads to nonperturbative results for quantities like the spin-spin correlation function. After a basic introduction to the phenomena and concomitant problems discussed in this thesis, a detailed description of the SFRG method in its initial formulation is given in the second chapter. We start with the generating functional of connected imaginary-time spin-correlation functions GΛ [h], for which an exact flow equation is derived. A particular issue, already pointed out by Krieg and Kopietz, arises here, namely the singular non-interacting limit of its subtracted Legendre transform ΓΛ [m]. As a consequence the initial condition of that functional does not have a proper series expansion in powers of m. This prevents us from working directly within a pure one-particle irreducible (1-PI) parametrization of the correlation functions, as is often done in the context of field theories. Thus motivated, we develop a workaround explicitly tailored to paramagnets, which provides us with a functional that has a well-behaved Legendre transform. The new approach is based on a different treatment of fluctuations at zero and finite frequencies, analogous to a previous hybrid formulation for the symmetry-broken phase. Certain properties, considered to be highly relevant for isotropic paramagnets, as well as previous observations, already made in the study of simpler spin systems like the Ising model, serve as additional justifications for choosing this construction.
In the third chapter our new method is assessed by calculating the dynamic susceptibility G(k, iω) and thus the dynamic structure factor S(k, ω) in the symmetric phase. For this purpose an approximate integral equation for the dynamic polarization function Π̃(k, iω) was derived. This equation results from a truncation of the hierarchy of flow equations and contains static quantities, that are assumed to be known from another source. Our first application is the high-temperature limit T → ∞ in d ≤ 3 dimensions. Salient features, believed to be part of the spin dynamics in isotropic Heisenberg magnets are also exhibited by our solution, like (anomalous) diffusion in a suitable hydrodynamic limit. Moreover we obtain the same order of magnitude for the diffusion coefficient D as in experiments and other theoretical calculations. Other aspects do not entirely agree with previous approaches.
Afterwards we continue by investigating systems close to the critical point Tc. Dynamic scaling forms for Π̃(k, iω) and S(k, ω), which, like spin diffusion, are postulated on the basis of quite general physical arguments, are reproduced. Agreement of the line-shapes 2with neutron scattering experiments at T = Tc is found to be satisfying, with deviations for ω → 0, that may be attributed to the simplicity of the approximation, like at infinite temperature.
Finally, we focus our attention on the thermodynamic properties of isotropic Heisenberg paramagnets by calculating the static susceptibility G(k). For this purpose we employ simple truncation schemes of the flow equations for the static self-energy ΣΛ (k) and four-spin vertex ΓΛ , together with a basic ansatz for the dynamic polarization Π̃(k, iω) in quantum systems. As a result we obtain transition temperatures Tc of three-dimensional nonfrustrated magnets within an accuracy of 5 percent compared to established benchmark values from Quantum Monte Carlo and high temperature expansion series. We conclude this chapter by giving an outlook on the application of our method to frustrated systems, which may require a combined non-trivial calculation of static and dynamic properties.
First study of reaction Ξ⁰n → Ξ⁻ p using Ξ⁰-nucleus scattering at an electron-positron collider
(2023)
Using ð1.0087 0.0044Þ × 1010 J=ψ events collected with the BESIII detector at the BEPCII storage ring, the process Ξ0n → Ξ−p is studied, where the Ξ0 baryon is produced in the process J=ψ → Ξ0Ξ¯ 0 and the neutron is a component of the 9 Be, 12C, and 197Au nuclei in the beam pipe. A clear signal is observed with a statistical significance of 7.1σ. The cross section of the reaction Ξ0 þ 9 Be → Ξ− þ p þ 8 Be is determined to be σðΞ0 þ 9 Be → Ξ− þ p þ 8 BeÞ¼ð22.1 5.3stat 4.5sysÞ mb at the Ξ0 momentum of 0.818 GeV=c, where the first uncertainty is statistical and the second is systematic. No significant H-dibaryon signal is observed in the Ξ−p final state. This is the first study of hyperon-nucleon interactions in electron-positron collisions and opens up a new direction for such research.
The quantum entangled J=ψ → ΣþΣ¯ − pairs from ð1.0087 0.0044Þ × 1010 J=ψ events taken by the BESIII detector are used to study the nonleptonic two-body weak decays Σþ → nπþ and Σ¯ − → n¯π−. The CP-odd weak decay parameters of the decays Σþ → nπþ (αþ) and Σ¯ − → n¯π− (α¯−) are determined to be 0.0481 0.0031stat 0.0019syst and −0.0565 0.0047stat 0.0022syst, respectively. The decay parameter α¯− is measured for the first time, and the accuracy of αþ is improved by a factor of 4 compared to the previous results. The simultaneously determined decay parameters allow the first precision CP symmetry test for any hyperon decay with a neutron in the final state with the measurement of ACP ¼ ðαþ þ α¯−Þ=ðαþ − α¯−Þ ¼ −0.080 0.052stat 0.028syst. Assuming CP conservation, the average decay parameter is determined as hαþi¼ðαþ − α¯−Þ=2 ¼ −0.0506 0.0026stat 0.0019syst, while the ratios αþ=α0 and α¯−=α¯ 0 are −0.0490 0.0032stat 0.0021syst and −0.0571 0.0053stat 0.0032syst, where α0 and α¯ 0 are the decay parameters of the decays Σþ → pπ0 and Σ¯ − → p¯ π0, respectively.
The quantum entangled J/ψ→Σ+Σ¯− pairs from (1.0087±0.0044)×1010 J/ψ events taken by the BESIII detector are used to study the non-leptonic two-body weak decays Σ+→nπ+ and Σ¯−→n¯π−. The CP-odd weak decay parameters of the decays Σ+→nπ+ (α+) and Σ¯−→n¯π− (α¯−) are determined to be −0.0565±0.0047stat±0.0022syst and 0.0481±0.0031stat±0.0019syst, respectively. The decay parameter α¯− is measured for the first time, and the accuracy of α+ is improved by a factor of four compared to the previous results. The simultaneously determined decay parameters allow the first precision CP symmetry test for any hyperon decay with a neutron in the final state with the measurement of ACP=(α++α¯−)/(α+−α¯−) = −0.080±0.052stat±0.028syst. Assuming CP conservation, the average decay parameter is determined as ⟨α+⟩=(α+−α¯−)/2 = −0.0506±0.0026stat±0.0019syst, while the ratios α+/α0 and α¯−/α¯0 are −0.0490±0.0032stat±0.0021syst and −0.0571±0.0053stat±0.0032syst, where α0 and α¯0 are the decay parameters of the decays Σ+→pπ0 and Σ¯−→p¯π0, respectively.
Precision measurements of the semileptonic decays 𝐷+𝑠→𝜂𝑒+𝜈𝑒 and 𝐷+𝑠→𝜂′𝑒+𝜈𝑒 are performed with 7.33 fb−1 of 𝑒+𝑒− collision data collected at center-of-mass energies between 4.128 and 4.226 GeV with the BESIII detector. The branching fractions obtained are ℬ(𝐷+𝑠→𝜂𝑒+𝜈𝑒) = (2.255±0.039stat±0.051syst)% and ℬ(𝐷+𝑠→𝜂′𝑒+𝜈𝑒)=(0.810±0.038stat±0.024syst)%. Combining these results with the ℬ(𝐷+→𝜂𝑒+𝜈𝑒) and ℬ(𝐷+→𝜂′𝑒+𝜈𝑒) obtained from previous BESIII measurements, the 𝜂−𝜂′ mixing angle in the quark flavor basis is determined to be 𝜙P=(40.0±2.0stat±0.6syst)°. Moreover, from the fits to the partial decay rates of 𝐷+𝑠→𝜂𝑒+𝜈𝑒 and 𝐷+𝑠→𝜂′𝑒+𝜈𝑒, the products of the hadronic transition form factors 𝑓𝜂(′)+(0) and the modulus of the 𝑐→𝑠 Cabibbo-Kobayashi-Maskawa matrix element |𝑉𝑐𝑠| are determined by using different hadronic transition form factor parametrizations. Based on the two-parameter series expansion, the products 𝑓𝜂+(0)|𝑉𝑐𝑠| = 0.4519±0.0071stat±0.0065syst and 𝑓𝜂′+(0)|𝑉𝑐𝑠| = 0.525±0.024stat±0.009syst are extracted. All results determined in this work supersede those measured in the previous BESIII analyses based on the 3.19 fb−1 subsample of data at 4.178 GeV.
Quantum-correlated 𝐷¯𝐷 pairs collected by the BESIII experiment at the 𝜓(3770) resonance corresponding to an integrated luminosity of 2.93 fb−1 are used to study the 𝐷0→𝐾0𝑆𝜋+𝜋−𝜋0 decay mode. The 𝐶𝑃-even fraction of 𝐷0→𝐾0𝑆𝜋+𝜋−𝜋0 decays is determined to be 0.235±0.010±0.002, where the first uncertainty is statistical and the second is systematic.
We present the first observation of the singly Cabibbo-suppressed decay Λ+c→ΛK+π0 with a significance of 5.7σ and the first evidence of Λ+c→ΛK+π+π− decay with a significance of 3.1σ, based on e+e− annihilation data recorded by the BESIII detector at the BEPCII collider. The data correspond to an integrated luminosity of 6.4 fb−1, in the center-of-mass energy range from 4.600 GeV to 4.950 GeV. We determine the branching fractions of Λ+c→ΛK+π0 and Λ+c→ΛK+π+π− relative to their Cabibbo-favored counterparts to be B(Λ+c→ΛK+π0)B(Λ+c→Λπ+π0)=(2.09±0.39stat.±0.07syst.)×10−2 and B(Λ+c→ΛK+π+π−)B(Λ+c→Λπ+π+π−)=(1.13±0.41stat.±0.06syst.)×10−2, respectively. Moreover, by combining our measured result with the world average of B(Λ+c→Λπ+π0), we obtain the branching fraction B(Λ+c→ΛK+π0)=(1.49±0.27stat.±0.05syst.±0.08ref.)×10−3. This result significantly departs from theoretical predictions based on quark SU(3) flavor symmetry, which is underpinned by the presumption of meson pair S-wave amplitude dominance.
A light scalar X0 or vector X1 particles have been introduced as a possible explanation for the (g−2)μ anomaly and dark matter phenomena.
Using (8.998±0.039)×109 $\jpsi$ events collected by the BESIII detector, we search for a light muon philic scalar X0 or vector X1 in the processes J/ψ→μ+μ−X0,1 with X0,1 invisible decays. No obvious signal is found, and the upper limits on the coupling g′0,1 between the muon and the X0,1 particles are set to be between 1.1×10−3 and 1.0×10−2 for the X0,1 mass in the range of 1<M(X0,1)<1000~MeV/c2 at 90% confidence level.
The Cabbibo-favored decay Λ+c→Ξ0K+π0 is studied for the first time using 6.1 fb−1 of e+e− collision data at center-of-mass energies between 4.600 and 4.840 GeV, collected with the BESIII detector at the BEPCII collider. With a double-tag method, the branching fraction of the three-body decay Λ+c→Ξ0K+π0 is measured to be (7.79±1.46±0.71)×10−3, where the first and second uncertainties are statistical and systematic, respectively. The branching fraction of the two-body decay Λ+c→Ξ(1530)0K+ is (5.99±1.04±0.29)×10−3, which is consistent with the previous result of (5.02±0.99±0.31)×10−3. In addition, the upper limit on the branching fraction of the doubly Cabbibo-suppressed decay Λ+c→nK+π0 is 7.1×10−4 at the 90% confidence level. The upper limits on the branching fractions of Λ+c→Σ0K+π0 and ΛK+π0 are also determined to be 1.8×10−3 and 2.0×10−3, respectively.
The Cabbibo-favored decay Λ+c→Ξ0K+π0 is studied for the first time using 6.1 fb−1 of e+e− collision data at center-of-mass energies between 4.600 and 4.840 GeV, collected with the BESIII detector at the BEPCII collider. With a double-tag method, the branching fraction of the three-body decay Λ+c→Ξ0K+π0 is measured to be (7.79±1.46±0.71)×10−3, where the first and second uncertainties are statistical and systematic, respectively. The branching fraction of the two-body decay Λ+c→Ξ(1530)0K+ is (5.99±1.04±0.29)×10−3, which is consistent with the previous result of (5.02±0.99±0.31)×10−3. In addition, the upper limit on the branching fraction of the doubly Cabbibo-suppressed decay Λ+c→nK+π0 is 7.1×10−4 at the 90% confidence level. The upper limits on the branching fractions of Λ+c→Σ0K+π0 and ΛK+π0 are also determined to be 1.8×10−3 and 2.0×10−3, respectively.
The family of cubic noncentrosymmetric 3-4-3 compounds has become a fertile ground for the discovery of novel correlated metallic and insulating phases. Here, we report the synthesis of a new heavy fermion compound, Ce3Bi4Ni3. It is an isoelectronic analog of the prototypical Kondo insulator Ce3Bi4Pt3 and of the recently discovered Weyl-Kondo semimetal Ce3Bi4Pd3. In contrast to the volume-preserving Pt-Pd substitution, structural and chemical analyses reveal a positive chemical pressure effect in Ce3Bi4Ni3 relative to its heavier counterparts. Based on the results of electrical resistivity, Hall effect, magnetic susceptibility, and specific heat measurements, we identify an energy gap of 65-70 meV, about 8 times larger than that in Ce3Bi4Pt3 and about 45 times larger than that of the Kondo-insulating background hosting the Weyl nodes in Ce3Bi4Pd3. We show that this gap as well as other physical properties do not evolve monotonically with increasing atomic number, i.e., in the sequence Ce3Bi4Ni3-Ce3Bi4Pd3-Ce3Bi4Pt3, but instead with increasing partial electronic density of states of the d orbitals at the Fermi energy. To understand under which condition topological states form in these materials is a topic for future studies.
Inverse modeling of circular lattices via orbit response measurements in the presence of degeneracy
(2023)
The number and location of beam position monitors (BPMs) and steerers with respect to the quadrupoles in a circular lattice can lead to degeneracy in the context of fitting linear optics and extracting lattice information from measured closed orbits. Furthermore, the measurement uncertainties due to the imperfection of BPMs and steerers can be propagated by the fitting process in ways that prohibit the successful extraction of discrepancies between lattice elements in the real machine and their description in the corresponding model. We systematically studied the influence of the placement of BPMs and steerers on the reconstruction of linear optics and corresponding lattice information. The derivative of orbit response coefficients with respect to the quadrupole strengths, the Jacobian, is derived as an analytical formula. This analytical version of the Jacobian is used to further derive the theoretical limitations of fitting linear optics from closed orbits in terms of the placement of BPMs and steerers. It is further demonstrated that when evaluating the Jacobian during the fitting procedure, the analytical version can be used in place of the conventional finite-difference computation. This allows for greatly improved efficiency when computing the Jacobian during each iteration of the fitting procedure. The approach is tested with large-scale simulations and the findings are verified by measurement data taken on SIS18 synchrotron at GSI Helmholtz Centre for Heavy Ion Research. The presented methods are of general nature and can be applied to other accelerator lattices as well. The fitting procedure by using the analytical Jacobian is tested in conjunction with various methods for mitigating quasidegeneracy and the results agree with those obtained by using the conventional Jacobian via finite-difference approximation.
We report a search for a dark photon using 14.9~fb−1 of e+e− annihilation data taken at center-of-mass energies from 4.13 to 4.60~GeV with the BESIII detector operated at the BEPCII storage ring. The dark photon is assumed to be produced in the radiative annihilation process of e+e− and to predominantly decay into light dark matter particles, which escape from the detector undetected. The mass range from 1.5 to 2.9~GeV is scanned for the dark photon candidate, and no significant signal is observed. The mass dependent upper limits at the 90% confidence level on the coupling strength parameter ϵ for a dark photon coupling with an ordinary photon vary between 1.6×10−3 and 5.7×10−3.
Using a data sample corresponding to an integrated luminosity of 11.3 fb−1 collected at center-of-mass energies from 4.23 to 4.70 GeV with the BESIII detector, we observe the process e+e− → π0π0ψ2(3823) for the first time with a statistical significance of 6.0 standard deviations. The ratio of average cross sections for e+e− → π0π0ψ2(3823) and π+π−ψ2(3823) is determined to be R = σ[e+e− → π0π0ψ2(3823)] σ[e+e−→π+π−ψ2(3823)] = 0.57 ± 0.14 ± 0.05, which is consistent with expectations from isospin symmetry. Here and below, the first uncertainties are statistical and the second are systematic. The mass of the ψ2(3823) is measured to be M[ψ2(3823)] = 3824.5±2.4±1.0 MeV/c2. Due to the limited data sample, an upper limit of 18.8 MeV at 90% confidence level is set on the intrinsic width of ψ2(3823).
We report the measurement of the cross sections for e+e−→hadrons at center-of-mass (c.m.) energies from 3.645 to 3.871 GeV. We observe a new resonance R(3810) in the cross sections for the first time, and observe the R(3760) resonance with high significance in the cross sections. The R(3810) has a mass of (3804.5±0.9±0.9) ~MeV/c2, a total width of (5.4±3.5±3.2)~MeV, and an electronic partial width of (19.4±7.4±12.1)~eV. Its significance is 7.7σ. The R(3810) could be interpreted as a hadro-charmonium resonance predicted by Quantum Chromodynamics (QCD). In addition, we measure the mass (3751.9±3.8±2.8) ~MeV/c2, the total width (32.8±5.8±8.7)~MeV, and the electronic partial width (184±75±86)~eV with improved precision for the R(3760). Furthermore, for the R(3780) we measure the mass (3778.7±0.5±0.3) ~MeV/c2 and total width (20.3±0.8±1.7)~MeV with improved precision, and the electronic partial width (265±69±83)~eV. The R(3780) can be interpreted as the 13D1 state of charmonium. Its mass and total width differ significantly from the corresponding fitted values given by the Particle Data Group in 2022 by 7.1 and 3.2 times the uncertainties for ψ(3770), respectively. ψ(3770) has been interpreted as the 13D1 state for 45 years.
Using 2.93 fb−1 of e+e− collision data collected with the BESIII detector at the center-of-mass energy 3.773 GeV, we perform the first amplitude analysis of the decay D+ → π+π0π0 and determine the relative magnitudes and phases of different intermediate processes. The absolute branching fraction of D+ → π+π0π0 is measured to be (2.888 ± 0.058stat. ± 0.069syst.)%. The dominant intermediate processes are D+ → a1(1260)+(→ ρ+π0) and D+ → *0ρ+, with branching fractions of (8.66 ± 1.04stat. ± 1.39syst.) × 10−3 and (9.70 ± 0.81stat. ± 0.53syst.) × 10−3, respectively.
We report a search for a heavier partner of the recently observed Zcs(3985)− state, denoted as Z′−cs, in the process e+e−→K+D∗−sD∗0+c.c., based on e+e− collision data collected at the center-of-mass energies of s√=4.661, 4.682 and 4.699 GeV with the BESIII detector. The Z′−cs is of interest as it is expected to be a candidate for a hidden-charm and open-strange tetraquark. A partial-reconstruction technique is used to isolate K+ recoil-mass spectra, which are probed for a potential contribution from Z′−cs→D∗−sD∗0 (c.c.). We find an excess of Z′−cs→D∗−sD∗0 (c.c.) candidates with a significance of 2.1σ, after considering systematic uncertainties, at a mass of (4123.5±0.7stat.±4.7syst.) MeV/c2. As the data set is limited in size, the upper limits are evaluated at the 90\% confidence level on the product of the Born cross sections (σBorn) and the branching fraction (B) of Z′−cs→D∗−sD∗0, under different assumptions of the Z′−cs mass from 4.120 to 4.140 MeV and of the width from 10 to 50 MeV at the three center-of-mass energies. The upper limits of σBorn⋅B are found to be at the level of O(1) pb at each energy. Larger data samples are needed to confirm the Z′−cs state and clarify its nature in the coming years.
We search for an axion-like particle (ALP) a through the process ψ(3686)→π+π−J/ψ, J/ψ→γa, a→γγ in a data sample of (2.71±0.01)×109 ψ(3686) events collected by the BESIII detector. No significant ALP signal is observed over the expected background, and the upper limits on the branching fraction of the decay J/ψ→γa and the ALP-photon coupling constant gaγγ are set at 95% confidence level in the mass range of 0.165≤ma≤2.84GeV/c2. The limits on B(J/ψ→γa) range from 8.3×10−8 to 1.8×10−6 over the search region, and the constraints on the ALP-photon coupling are the most stringent to date for 0.165 ≤ ma ≤ 1.468GeV/c2.
We search for an axion-like particle (ALP) a through the process ψ(3686)→π+π−J/ψ, J/ψ→γa, a→γγ in a data sample of (2.71±0.01)×109 ψ(3686) events collected by the BESIII detector. No significant ALP signal is observed over the expected background, and the upper limits on the branching fraction of the decay J/ψ→γa and the ALP-photon coupling constant gaγγ are set at 95% confidence level in the mass range of 0.165≤ma≤2.84GeV/c2. The limits on B(J/ψ→γa) range from 8.3×10−8 to 1.8×10−6 over the search region, and the constraints on the ALP-photon coupling are the most stringent to date for 0.165≤ma≤1.468GeV/c2.
The cross sections of the 𝑒+𝑒−→𝜙𝜂′ process at center-of-mass energies from 3.508 to 4.951 GeV are measured with high precision using 26.1 fb−1 data collected with the BESIII detector operating at the BEPCII storage ring. The cross sections are of the order of a few picobarn and decrease as the center-of-mass energy increases as 𝑠−𝑛/2 with 𝑛=4.35±0.14. This result is in agreement with the Nambu-Jona-Lasinio model prediction of 𝑛=3.5±0.9. In addition, the charmless decay 𝜓(3770)→𝜙𝜂′ is searched for by fitting the measured cross sections, yet no significant signal is observed. The upper limit of ℬ(𝜓(3770)→𝜙𝜂′) at the 90% confidence level is determined to be 2.3×10−5.
With data samples collected with the BESIII detector at seven energy points at √s = 3.68 − 3.71 GeV, corresponding to an integrated luminosity of 333 pb−1, we present a study of the Λ transverse polarization in the e+e− → ΛΛ¯ reaction. The signifcance of polarization by combining the seven energy points is found to be 2.6σ including the systematic uncertainty, which implies a non-zero phase between the transition amplitudes of the ΛΛ¯ helicity states. The modulus ratio and the relative phase of EM-psionic form factors combined with all energy points are measured to be RΨ = 0.71+0.10−0.10 ± 0.03 and ∆ΦΨ = 23+8.8−8.0 ± 1.6◦, where the frst uncertainties are statistical and the second systematic.
Efficient modeling and mitigation of quadrupole errors in synchrotrons and their beam transfer lines
(2023)
This thesis investigates the problem of estimating quadrupole errors on synchrotrons as well as how to minimize the influence of quadrupole errors for beam transfer lines (beamlines). It emphasizes the importance to treat possible error sources in all parts of an accelerator in order to provide constantly high beam quality to the experimental stations. While the presented methods have been investigated by using the example of the SIS18 synchrotron and the HEST beamlines at GSI Helmholtz Centre for Heavy Ion Research, they are equally relevant for the future synchrotrons and beamlines of the Facility for Antiproton and Ion Research in Europe (FAIR).
Part 1 discusses the problem of estimating quadrupole errors via orbit response measurements at synchrotrons. An emphasis is put on investigating the influence of the availability of steerer magnets and beam position monitors (BPMs) on the solvability of the inverse problem as well as on the propagation of measurement uncertainty for the estimation of quadrupole errors. The problem is approached via analytical considerations as well as via dedicated simulation studies. By developing an analytical expression for the Jacobian matrix, the theoretical boundaries for the solvability of the inverse problem are derived. Moreover, it is shown that the analytical expressions for the Jacobian matrix can be used during the fitting procedure to achieve a significant improvement in the computational efficiency by a factor $N_{steerers} \times N_{quadrupoles}$, where $N$ denotes the number of lattice elements of the respective type. The presented results are tested via dedicated measurements at the SIS18 synchrotron.
Part 2 discusses – complementary to part 1 – the influence of quadrupole errors in beam transfer lines with respect to the beam quality requirements given by the experimental stations. A preventive approach is presented which allows to minimize the influence of possible quadrupole errors on the degradation of beam quality. By identifying and selecting robust quadrupole configurations, a stable operation of the beamline can be enabled and the time needed by operators to readjust the beamline parameters can be reduced. The concept of beamline robustness is developed and is studied with the help of dedicated simulations. The simulation results are used to identify certain properties that distinguish robust from nonrobust quadrupole configurations. Also, various methods for improving the computational process of identifying robust quadrupole configurations are presented. The methods and results are tested via dedicated measurements at two different beamlines at GSI Helmholtz Centre for Heavy Ion Research and at Forschungszentrum Jülich.
Using 15.6 fb−1 of e+e− collision data collected at twenty-four center-of-mass energies from 4.0 to 4.6 GeV with the BESIII detector, the helicity amplitudes of the process e+e− → π+π−ω are analyzed for the first time. Born cross section measurements of two-body intermediate resonance states with statistical significance greater than 5σ are presented, such as f0(500), f0(980), f2(1270), f0(1370), b1(1235)±, and ρ(1450)±. In addition, evidence of a resonance state in e+e− → π+π−ω production is found. The mass of this state obtained by line shape fitting is about 4.2 GeV/c2, which is consistent with the production of ψ(4160) or Y(4220).
Using 15.6 fb−1 of e+e− collision data collected at twenty-four center-of-mass energies from 4.0 to 4.6 GeV with the BESIII detector, the helicity amplitudes of the process e+e−→π+π−ω are analyzed for the first time. Born cross section measurements of two-body intermediate resonance states with statistical significance greater than 5σ are presented, such as f0(500), f0(980), f2(1270), f0(1370), b1(1235)±, and ρ(1450)±. In addition, evidence of a resonance state in e+e−→π+π−ω production is found. The mass of this state obtained by line shape fitting is about 4.2 GeV/c2, which is consistent with the production of ψ(4160) or Y(4220).
Using an 𝑒+𝑒− collision data sample of (27.08±0.14)×108 𝜓(3686) events collected by the BESIII detector, we report the first observation of 𝜒𝑐𝐽→Ω−¯Ω+ (𝐽=0, 1, 2) decays with significances of 5.6𝜎, 6.4𝜎, and 18𝜎, respectively, where the 𝜒𝑐𝐽 mesons are produced in the radiative 𝜓(3686) decays. The branching fractions are determined to be ℬ(𝜒𝑐0→Ω−¯Ω+) = (3.51±0.54±0.29)×10−5, ℬ(𝜒𝑐1→Ω−¯Ω+)=(1.49±0.23±0.10)×10−5, and ℬ(𝜒𝑐2→Ω−¯Ω+)=(4.52±0.24±0.18)×10−5, where the first and second uncertainties are statistical and systematic, respectively.
Using (1.0087±0.0044)×1010 𝐽/𝜓 events collected by the BESIII detector at the BEPCII collider, we report the first search for the baryon and lepton number violating decays Ξ0→𝐾−𝑒+ with Δ(𝐵−𝐿)=0 and Ξ0→𝐾+𝑒− with |Δ(𝐵−𝐿)|=2, where 𝐵 (𝐿) is the baryon (lepton) number. While no signal is observed, the upper limits on the branching fractions of these two decays are set to ℬ(Ξ0→𝐾−𝑒+)<3.6×10−6 and ℬ(Ξ0→𝐾+𝑒−)<1.9×10−6 at the 90% confidence level, respectively. These results offer a direct probe of baryon number violating interactions involving a strange quark.
The Ξ0 asymmetry parameters are measured using entangled quantum Ξ0 − Ξ¯ 0 pairs from a sample of ð448.1 2.9Þ × 106 ψð3686Þ events collected with the BESIII detector at BEPCII. The relative phase between the transition amplitudes of the Ξ0Ξ¯ 0 helicity states is measured to be ΔΦ ¼ −0.050 0.150 0.020 rad, which implies that there is no obvious polarization at the current level of statistics. The decay parameters of the Ξ0 hyperon ðαΞ0 ; αΞ¯ 0 ; ϕΞ0 ; ϕΞ¯ 0 Þ and the angular distribution parameter ½αψð3686Þ and ΔΦ are measured simultaneously for the first time. In addition, the CP asymmetry observables are determined to be AΞ0 CP ¼ ðαΞ0 þ αΞ¯ 0 Þ=ðαΞ0 − αΞ¯ 0 Þ ¼ −0.007 0.082 0.025 and ΔϕΞ0 CP ¼ ðϕΞ0 þ ϕΞ¯ 0 Þ=2 ¼ −0.079 0.082 0.010 rad, which are consistent with CP conservation.
We report on the first search for ¯Λ−Λ oscillations in the decay 𝐽/𝜓→𝑝𝐾−¯Λ+c.c. by analyzing 1.31×109 𝐽/𝜓 events accumulated with the BESIII detector at the BEPCII collider. The 𝐽/𝜓 events are produced using 𝑒+𝑒− collisions at a center of mass energy √𝑠=3.097 GeV. No evidence for hyperon oscillations is observed. The upper limit for the oscillation rate of ¯Λ to Λ hyperons is determined to be 𝒫(Λ)=[ℬ(𝐽/𝜓→𝑝𝐾−Λ+c.c.)/ℬ(𝐽/𝜓→𝑝𝐾−¯Λ+c.c.)]<4.4×10−6 corresponding to an oscillation parameter 𝛿𝑚Λ¯Λ of less than 3.8×10−18 GeV at the 90% confidence level.
The theoretical and experimental investigation of exotic hadrons like tetraquarks is an important branch of modern elementary particle physics. In this thesis I investigate different four-quark systems using lattice QCD and search for evidence of stable tetraquark states or resonances.
Lattice QCD as a non-perturbative approach to QCD allows an accurate and reliable determination of the masses of strongly bound hadrons.
However, most tetraquarks appear as weakly bound states or resonances, which makes a theoretical investigation using lattice QCD difficult due to the finite spatial volume. A rigorous treatment of such systems is feasible using the so-called Lüscher method. This allows to calculate the scattering amplitude based on the finite-volume energy spectrum determined in a lattice QCD calculation. Similarly to the analysis of experimental data, this scattering amplitude can be used to determine the binding energies of bound states or the masses and decay widths of resonances in the infinite volume.
In my work I calculate the low-energy energy spectra of different four-quark systems and use - if necessary - the Lüscher method to determine the masses of potential tetraquark states.
I focus on systems consisting of two heavy antiquarks and two light quarks, where at least one of the heavy antiquarks is a bottom quark.
Even though such tetraquarks have not yet been experimentally detected, they are considered promising candidates for particles that are stable with respect to the strong interaction.
A decisive step for successfully calculating low-lying energy levels for such four-quark systems is a carefully chosen set of creation operators, which represent the physical states most accurately. In addition to operators that generate a local structure where all four quarks are located at the same space-time point, I also use so-called scattering operators that resemble two spatially separated mesons. These scattering operators turned out to be relevant for successfully determining the lowest energy levels and are therefore essential, especially if a Lüscher analysis is carried out.
In my work, I considered two different lattice setups to study the four-quark systems $\bar{b}\bar{b}ud$ with $I(J^P)=0(1^+) $, $\bar{b}\bar{b}us$ with $J^P=1^+ $ and $\bar{b}\bar{c}ud$ with $I(J^P)=0(0^+) $ and $I(J^P)=0(1^+) $ and to predict potential tetraquark states. In both setups, I considered scattering operators. While in the first setup I used them only as annihilation operators, in the second setup they were included both as creation and annihilation operators. Additionally, in the second lattice setup, I performed a simplified investigation of the $\bar{b}\bar{b}ud$ system with $I(J^P)=0(1^-) $, which is a potential candidate for a tetraquark resonance. The results of the investigation of the mentioned four-quark systems can be summarized as follows:
For the $ \bar{b}\bar{b}ud $ four-quark system with $ I(J^P)=0(1^+) $ I found a deeply bound ground state slightly more than $ 100\,\textrm{MeV} $ below the lowest meson-meson threshold. The existence of a corresponding $\bar{b}\bar{b}ud$ tetraquark in the infinite volume was confirmed using a Lüscher analysis and possible systematic errors due to the use of lattice QCD were taken into account.
Similar results were obtained for the $ \bar{b}\bar{b}us $ four-quark system with $ J^P=1^+ $. Again, I found a ground state well below the lowest meson-meson threshold, but slightly weaker bound than for the $ \bar{b}\bar{b}ud $ system. Effects due to the finite volume turned out to be negligible for this system, as already predicted for the $ \bar{b}\bar{b}ud $ system. \item For the $ \bar{b}\bar{c}ud $ four-quark systems with $ (J^P)=0(0^+) $ and $ (J^P)=0(1^+) $ I was able to rule out the existence of a deeply bound tetraquark states based on the energy spectrum in the finite volume. However, by means of a scattering analysis using the Lüscher method, I found evidence a broad resonance for both channels.
In the case of the $ \bar{b}\bar{b}ud $ four-quark system with $ I(J^P)=0(1^-) $, I could neither confirm the existence of a resonance, nor rule out its existence with certainty.
In particular, my investigations showed that the results of the two different lattice simulations are consistent. The theoretical prediction of the bound tetraquark states $\bar{b}\bar{b}ud$ and $\bar{b}\bar{b}us$ as well as the tetraquark resonances in the $\bar{b}\bar{c}ud$ system in this work represent an important contribution to the future experimental search for exotic hadrons and can support the discovery of previously unobserved particles.
We search for an axion-like particle (ALP) a through the process ψ(3686)→π+π−J/ψ, J/ψ→γa, a→γγ in a data sample of (2.71±0.01)×109 ψ(3686) events collected by the BESIII detector. No significant ALP signal is observed over the expected background, and the upper limits on the branching fraction of the decay J/ψ→γa and the ALP-photon coupling constant gaγγ are set at 95% confidence level in the mass range of 0.165≤ma≤2.84GeV/c2. The limits on B(J/ψ→γa) range from 8.3×10−8 to 1.8×10−6 over the search region, and the constraints on the ALP-photon coupling are the most stringent to date for 0.165≤ma≤1.468GeV/c2.
Using data samples with an integrated luminosity of 22.42 fb−1 collected by the BESIII detector operating at the BEPCII storage ring, we measure the cross sections of the $e^{+}e^{-}\rightarrow\etaJ/\psi$ process at center-of-mass energies from 3.808 to 4.951 GeV. Three structures are observed in the line shape of the measured cross sections. A maximum-likelihood fit with ψ(4040), two additional resonances, and a non-resonant component is performed. The mass and width of the first additional state are (4219.7±2.5±4.5)MeV/c2 and (80.7±4.4±1.4)MeV, respectively, consistent with the ψ(4230). For the second state, the mass and width are (4386±13±17)MeV/c2 and (177±32±13)MeV, respectively, consistent with the ψ(4360). The first uncertainties are statistical and the second ones are systematic. The statistical significance of ψ(4040) is 8.0σ and those for ψ(4230) and ψ(4360) are more than 10.0σ.
Using a sample of (10087±44)×106 J/ψ events, which is about fifty times larger than that was previously analyzed, a further investigation on the J/ψ→γ3(π+π−) decay is performed. A significant distortion at 1.84 GeV/c2 in the line-shape of the 3(π+π−) invariant mass spectrum is observed for the first time, which is analogous to the behavior of X(1835) and could be resolved by two overlapping resonant structures, X(1840) and X(1880). The new state X(1880) is observed with a statistical significance of 14.7σ. The mass and width of X(1880) are determined to be 1882.1±1.7±0.7 MeV/c2 and 30.7±5.5±2.4 MeV, respectively, which indicates the existence of a pp¯ bound state.
Based on (10.09±0.04)×109 J/ψ events collected with the BESIII detector operating at the BEPCII collider, a partial wave analysis of the decay J/ψ→ϕπ0η is performed. We observe for the first time two new structures on the ϕη invariant mass distribution, with statistical significances of 24.0σ and 16.9σ; the first with JPC = 1+−, mass M = (1911 ± 6 (stat.) ± 14 (sys.))~MeV/c2, and width Γ= (149 ± 12 (stat.) ± 23 (sys.))~MeV, the second with JPC = 1−−, mass M = (1996 ± 11 (stat.) ± 30 (sys.))~MeV/c2, and width Γ = (148 ± 16 (stat.) ± 66 (sys.))~MeV. These measurements provide important input for the strangeonium spectrum. In addition, the f0(980)−a0(980)0 mixing signal in J/ψ→ϕf0(980)→ϕa0(980)0 and the corresponding electromagnetic decay J/ψ→ϕa0(980)0 are measured with improved precision, providing crucial information to understand the nature of a0(980)0 and f0(980).
Based on (10.09±0.04)×109 J/ψ events collected with the BESIII detector operating at the BEPCII collider, a partial wave analysis of the decay J/ψ→ϕπ0η is performed. We observe for the first time two new structures on the ϕη invariant mass distribution, with statistical significances of 24.0σ and 16.9σ; the first with JPC = 1+−, mass M = (1911 ± 6 (stat.) ± 14 (sys.))~MeV/c2, and width Γ= (149 ± 12 (stat.) ± 23 (sys.))~MeV, the second with JPC = 1−−, mass M = (1996 ± 11 (stat.) ± 30 (sys.))~MeV/c2, and width Γ = (148 ± 16 (stat.) ± 66 (sys.))~MeV. These measurements provide important input for the strangeonium spectrum. In addition, the f0(980)−a0(980)0 mixing signal in J/ψ→ϕf0(980)→ϕa0(980)0 and the corresponding electromagnetic decay J/ψ→ϕa0(980)0 are measured with improved precision, providing crucial information to understand the nature of a0(980)0 and f0(980).
Using 7.33 fb−1 of e+e− collision data collected by the BESIII detector at center-of-mass energies between 4.128 and 4.226~GeV, we observe for the first time the decay D±s→ωπ±η with a statistical significance of 7.6σ. The measured branching fraction of this decay is (0.54±0.12±0.04)%, where the first uncertainty is statistical and the second is systematic.
Evidence for the singly Cabibbo suppressed decay Λ+c→pπ0 is reported for the first time with a statistical significance of 3.7σ based on 6.0 fb−1 of e+e− collision data collected at center-of-mass energies between 4.600 and 4.843 GeV with the BESIII detector at the BEPCII collider. The absolute branching fraction of Λ+c→pπ0 is measured to be (1.56+0.72−0.58±0.20)×10−4, which distinctly exceeds the upper limit measured by Belle experiment. Combining with the branching fraction of Λ+c→nπ+, (6.6±1.3)×10−4, the ratio of the branching fractions of Λ+c→nπ+ and Λ+c→pπ0 is calculated to be 3.2+2.2−1.2. As an important input for the theoretical models describing the decay mechanisms of charmed baryons, our result indicates that the non-factorizable contributions play an essential role and their interference with the factorizable contributions should not be significant. In addition, the absolute branching fraction of Λ+c→pη is measured to be (1.63±0.31stat±0.11syst)×10−3, which is consistent with previous measurements.
Evidence for the singly Cabibbo suppressed decay Λ+c→pπ0 is reported for the first time with a statistical significance of 3.7σ based on 6.0 fb−1 of e+e− collision data collected at center-of-mass energies between 4.600 and 4.843 GeV with the BESIII detector at the BEPCII collider. The absolute branching fraction of Λ+c→pπ0 is measured to be (1.56+0.72−0.58±0.20)×10−4. Combining with the branching fraction of Λ+c→nπ+, (6.6±1.2±0.4)×10−4, the ratio of the branching fractions Λ+c→nπ+ and Λ+c→pπ0 is calculated to be 4.2+2.2−1.9; this is an important input for the understanding of the decay mechanisms of charmed baryons. In addition, the absolute branching fraction of Λ+c→pη is measured to be (1.63±0.31stat±0.11syst)×10−3, which is consistent with previous measurements.
Based on 4.5 fb−1 data taken at seven center-of-mass energies ranging from 4.600 to 4.699 GeV with the BESIII detector at the BEPCII collider, we measure the branching fractions of Λ + c → Σ + + hadrons relative to Λ + c → Σ +π +π −. Combining with the world average branching fraction of Λ + c → Σ +π +π −, their branching fractions are measured to be (0.377 ± 0.042 ± 0.020 ± 0.021)% for Λ + c → Σ +K+K−, (0.200 ± 0.023 ± 0.011 ± 0.011)% for Λ + c → Σ+K+π−, (0.414 ± 0.080 ± 0.030 ± 0.023)% for Λ + c → Σ +φ and (0.197 ± 0.036 ± 0.009 ± 0.011)% for Λ + c → Σ +K+K−(non-φ). In all the above results, the first uncertainties are statistical, the second are systematic and the third are from external input of the branching fraction of Λ + c → Σ +π +π −. Since no signal for Λ + c → Σ +K+π−π 0 is observed, the upper limit of its branching fraction is determined to be 0.13% at the 90% confidence level.
We search for the semi-leptonic decays Λ + c → Λπ+π−e+νe and Λ + c → pK0 Sπ−e+νe in a sample of 4.5 fb−1 of e+e− annihilation data collected in the center-of-mass energy region between 4.600 GeV and 4.699 GeV by the BESIII detector at the BEPCII. No significant signals are observed, and the upper limits on the decay branching fractions are set to be B(Λ+c → Λπ+π−e+νe ) < 3.9 × 10−4 and B(Λ + c → pK0Sπ−e+νe ) < 3.3 × 10−4 at the 90% confidence level, respectively.
We present cross sections for the reaction e+e−→K0SK0L at center-of-mass energies ranging from 3.51 GeV to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5 fb−1. The ratio of neutral-to-charged kaon form factors at large momentum transfers (12 GeV2<Q2<25 GeV2) is determined to be 0.21±0.01, which indicates a small but significant effect of flavor-SU(3) breaking in the kaon wave function, and consequently excludes the possibility that flavor-SU(3) breaking is the primary reason for the strong experimental violation of the pQCD prediction |F(π±)|/|F(K±)|=f2π/f2K, where F(π±) and F(K±) are the form factors, and fπ and fK are the decay constants of charged pions and kaons, respectively. We also observe a significant signal for the charmless decay ψ(3770)→K0SK0L for the first time. Within a 1σ contour of the likelihood value, the the branching fraction for ψ(3770)→K0SK0L is determined to be B=(2.63+1.40−1.59)×10−5, and the relative phase between the continuum and ψ(3770) amplitudes is ϕ=(−0.39+0.05−0.10)π. The branching fraction is in good agreement with the S- and D-wave charmonia mixing scheme proposed in the interpretation of the "ρπ puzzle" between J/ψ and ψ(3686) decays.
A search has been performed for the semileptonic decays D0→K0SK−e+νe, D+→K0SK0Se+νe and D+→K+K−e+νe, using 7.9 fb−1 of e+e− annihilation data collected at the center-of-mass energy s√=3.773 GeV by the BESIII detector operating at the BEPCII collider. No significant signals are observed, and upper limits are set at the 90\% confidence level of 2.13×10−5, 1.54×10−5 and 2.10×10−5 for the branching fractions of D0→K0SK−e+νe, D+→K0SK0Se+νe and D+→K+K−e+νe, respectively.
The Heidelberg Ion-Beam Therapy Centre (HIT) provides proton, helium, and carbon-ion beams with different energies and intensities for cancer treatment and oxygen-ion beams for experiments. For several experiments and possible future applications, such as helium ion beam radiography, a low-intensity ion beam monitor integrated into the dose delivery feedback system for the accelerator control is a necessary pre-requisite. The updated 2D prototype for this purpose consists of scintillating fibres with enhanced radiation hardness, silicon photomultipliers (SiPMs) to amplify the emitted light, and a dedicated front-end readout system (FERS) to process and record the generated signals. This setup was tested successfully on monitoring ion-beam position and profile horizontally and vertically, as well as the beam intensity, for all four ion types with energies from 50 to 430 MeV/u and intensities from 1E2 to 1E7 ions/s. Additionally, time-of-arrival (ToA) measurements on single ions have been successfully performed for a limited intensity range, allowing for ion tracking in a further update. This will reduce noise, and will also improve the accuracy and usability of ion radiography.
The free energy of TAP-solutions for the SK-model of mean field spin glasses can be expressed as a nonlinear functional of local terms: we exploit this feature in order to contrive abstract REM-like models which we then solve by a classical large deviations treatment. This allows to identify the origin of the physically unsettling quadratic (in the inverse of temperature) correction to the Parisi free energy for the SK-model, and formalizes the true cavity dynamics which acts on TAP-space, i.e. on the space of TAP-solutions. From a non-spin glass point of view, this work is the first in a series of refinements which addresses the stability of hierarchical structures in models of evolving populations.
Off-central heavy-ion collisions are known to feature magnetic fields with magnitudes and characteristic gradients corresponding to the scale of the strong interactions. In this work, we employ equilibrium lattice simulations of the underlying theory, QCD, involving similar inhomogeneous magnetic field profiles to achieve a better understanding of this system. We simulate three flavors of dynamical staggered quarks with physical masses at a range of magnetic fields and temperatures, and extrapolate the results to the continuum limit. Analyzing the impact of the field on the quark condensate and the Polyakov loop, we find non-trivial spatial features that render the QCD medium qualitatively different as in the homogeneous setup, especially at temperatures around the transition. In addition, we construct leading-order chiral perturbation theory for the inhomogeneous background and compare its prediction to our lattice results at low temperature. Our findings will be useful to benchmark effective theories and low-energy models of QCD for a better description of peripheral heavy-ion collisions.
Off-central heavy-ion collisions are known to feature magnetic fields with magnitudes and characteristic gradients corresponding to the scale of the strong interactions. In this work, we employ equilibrium lattice simulations of the underlying theory, QCD, involving similar inhomogeneous magnetic field profiles to achieve a better understanding of this system. We simulate three flavors of dynamical staggered quarks with physical masses at a range of magnetic fields and temperatures, and extrapolate the results to the continuum limit. Analyzing the impact of the field on the quark condensate and the Polyakov loop, we find non-trivial spatial features that render the QCD medium qualitatively different as in the homogeneous setup, especially at temperatures around the transition. In addition, we construct leading-order chiral perturbation theory for the inhomogeneous background and compare its prediction to our lattice results at low temperature. Our findings will be useful to benchmark effective theories and low-energy models of QCD for a better description of peripheral heavy-ion collisions.
Off-central heavy-ion collisions are known to feature magnetic fields with magnitudes and characteristic gradients corresponding to the scale of the strong interactions. In this work, we employ equilibrium lattice simulations of the underlying theory, QCD, involving similar inhomogeneous magnetic field profiles to achieve a better understanding of this system. We simulate three flavors of dynamical staggered quarks with physical masses at a range of magnetic fields and temperatures, and extrapolate the results to the continuum limit. Analyzing the impact of the field on the quark condensate and the Polyakov loop, we find non-trivial spatial features that render the QCD medium qualitatively different as in the homogeneous setup, especially at temperatures around the transition. In addition, we construct leading-order chiral perturbation theory for the inhomogeneous background and compare its prediction to our lattice results at low temperature. Our findings will be useful to benchmark effective theories and low-energy models of QCD for a better description of peripheral heavy-ion collisions.
Human feline leukemia virus subgroup C receptor-related proteins 1 and 2 (FLVCR1 and 2) are members of the major facilitator superfamily1. Their dysfunction is linked to several clinical disorders, including PCARP, HSAN, and Fowler syndrome2–7. Earlier studies concluded that FLVCR1 may function as a putative heme exporter8–12, while FLVCR2 was suggested to act as a heme importer13, yet conclusive biochemical and detailed molecular evidence remained elusive for the function of both transporters14–17. Here, we show that FLVCR1 and FLVCR2 facilitate the transport of choline and ethanolamine across human plasma membranes, utilizing a concentration-driven substrate translocation process. Through structural and computational analyses, we have identified distinct conformational states of FLVCRs and unraveled the coordination chemistry underlying their substrate interactions. Within the binding pocket of both transporters, we identify fully conserved tryptophan and tyrosine residues holding a central role in the formation of cation-π interactions, essential for choline and ethanolamine selectivity. Our findings not only clarify the mechanisms of choline and ethanolamine transport by FLVCR1 and FLVCR2, enhancing our comprehension of disease-associated mutations that interfere with these vital processes, but also shed light on the conformational dynamics of these MFS-type proteins during the transport cycle.
The polarization of Λ and Λ¯ hyperons along the beam direction has been measured relative to the second and third harmonic event planes in isobar Ru+Ru and Zr+Zr collisions at √sNN = 200 GeV. This is the first experimental evidence of the hyperon polarization by the triangular flow originating from the initial density fluctuations. The amplitudes of the sine modulation for the second and third harmonic results are comparable in magnitude, increase from central to peripheral collisions, and show a mild pT dependence. The azimuthal angle dependence of the polarization follows the vorticity pattern expected due to elliptic and triangular anisotropic flow, and qualitatively disagree with most hydrodynamic model calculations based on thermal vorticity and shear induced contributions. The model results based on one of existing implementations of the shear contribution lead to a correct azimuthal angle dependence, but predict centrality and pT dependence that still disagree with experimental measurements. Thus, our results provide stringent constraints on the thermal vorticity and shear-induced contributions to hyperon polarization. Comparison to previous measurements at RHIC and the LHC for the second-order harmonic results shows little dependence on the collision system size and collision energy.
The production of K∗(892)± meson resonance is measured at midrapidity (|y|<0.5) in Pb-Pb collisions at sNN−−−√=5.02 TeV using the ALICE detector at the LHC. The resonance is reconstructed via its hadronic decay channel K∗(892)±→K0Sπ±. The transverse momentum distributions are obtained for various centrality intervals in the pT range of 0.4-16 GeV/c. The reported measurements of integrated yields, mean transverse momenta, and particle yield ratios are consistent with previous ALICE measurements for K∗(892)0. The pT-integrated yield ratio 2K∗(892)±/(K++K−) in central Pb-Pb collisions shows a significant suppression (9.3σ) relative to pp collisions. Thermal model calculations overpredict the particle yield ratio. Although both simulations consider the hadronic phase, only HRG-PCE accurately represents the measurements, whereas MUSIC+SMASH tends to overpredict them. These observations, along with the kinetic freeze-out temperatures extracted from the yields of light-flavored hadrons using the HRG-PCE model, indicate a finite hadronic phase lifetime, which increases towards central collisions. The pT-differential yield ratios 2K∗(892)±/(K++K−) and 2K∗(892)±/(π++π−) are suppressed by up to a factor of five at pT<2 GeV/c in central Pb-Pb collisions compared to pp collisions at s√= 5.02 TeV. Both particle ratios and are qualitatively consistent with expectations for rescattering effects in the hadronic phase. The nuclear modification factor shows a smooth evolution with centrality and is below unity at pT>8 GeV/c, consistent with measurements for other light-flavored hadrons. The smallest values are observed in most central collisions, indicating larger energy loss of partons traversing the dense medium.
A new, more precise measurement of the Λ hyperon lifetime is performed using a large data sample of Pb–Pb collisions at √sNN p ¼ 5.02 TeV with ALICE. The Λ and Λ¯ hyperons are reconstructed at midrapidity using their two-body weak decay channel Λ → p þ π− and Λ¯ → p¯ þ πþ. The measured value of the Λ lifetime is τΛ ¼ ½261.07 0.37ðstat:Þ 0.72ðsyst:Þ ps. The relative difference between the lifetime of Λ and Λ¯ , which represents an important test of CPT invariance in the strangeness sector, is also measured. The obtained value ðτΛ − τΛ¯Þ=τΛ ¼ 0.0013 0.0028ðstat:Þ 0.0021ðsyst:Þ is consistent with zero within the uncertainties. Both measurements of the Λ hyperon lifetime and of the relative difference between τΛ and τΛ¯ are in agreement with the corresponding world averages of the Particle Data Group and about a factor of three more precise.
First measurement of Λ+c production down to pT = 0 in pp and p-Pb collisions at √𝑠NN = 5.02 TeV
(2023)
The production of prompt +c baryons has been measured at midrapidity in the transverse momentum interval 0 < pT < 1 GeV/c for the first time, in pp and p–Pb collisions at a center-of-mass energy per nucleon-nucleon collision √sNN = 5.02 TeV. The measurement was performed in the decay channel +c → pK0S by applying new decay reconstruction techniques using a Kalman-Filter vertexing algorithm and adopting a machine-learning approach for the candidate selection. The pT -integrated +c production cross sections in both collision systems were determined and used along with the measured yields in Pb–Pb collisions to compute the pT -integrated nuclear modification factors RpPb and RAA of +c baryons, which are compared to model calculations that consider nuclear modification of the parton distribution functions. The +c /D0 baryon-to-meson yield ratio is reported for pp and p–Pb collisions. Comparisons with models that include modified hadronization processes are presented, and the implications of the results on the understanding of charm hadronization in hadronic collisions are discussed. A significant (3.7σ) modification of the mean transverse momentum of + c baryons is seen in p–Pb collisions with respect to pp collisions, while the pT -integrated +c /D0 yield ratio was found to be consistent between the two collision systems within the uncertainties.
The measurement of the production of deuterons, tritons and 3He and their antiparticles in Pb-Pb collisions at √sNN = 5.02 TeV is presented in this article. The measurements are carried out at midrapidity (y|< 0.5) as a function of collision centrality using the ALICE detector. The pT-integrated yields, the coalescence parameters and the ratios to protons and antiprotons are reported and compared with nucleosynthesis models. The comparison of these results in different collision systems at different center-of-mass collision energies reveals a suppression of nucleus production in small systems. In the Statistical Hadronisation Model framework, this can be explained by a small correlation volume where the baryon number is conserved, as already shown in previous fluctuation analyses. However, a different size of the correlation volume is required to describe the proton yields in the same data sets. The coalescence model can describe this suppression by the fact that the wave functions of the nuclei are large and the fireball size starts to become comparable and even much smaller than the actual nucleus at low multiplicities.
The knowledge of the material budget with a high precision is fundamental for measurements of direct photon production using the photon conversion method due to its direct impact on the total systematic uncertainty. Moreover, it influences many aspects of the charged-particle reconstruction performance. In this article, two procedures to determine data-driven corrections to the material-budget description in ALICE simulation software are developed. One is based on the precise knowledge of the gas composition in the Time Projection Chamber. The other is based on the robustness of the ratio between the produced number of photons and charged particles, to a large extent due to the approximate isospin symmetry in the number of produced neutral and charged pions. Both methods are applied to ALICE data allowing for a reduction of the overall material budget systematic uncertainty from 4.5% down to 2.5%. Using these methods, a locally correct material budget is also achieved. The two proposed methods are generic and can be applied to any experiment in a similar fashion.
We study equilibrium as well as out-of-equilibrium properties of the strongly interacting QGP medium under extreme conditions of high temperature T and high baryon densities or baryon chemical potentials μB within a kinetic approach. We present the thermodynamic and transport properties of the QGP close to equilibrium in the framework of effective models with Nf=3 active quark flavours such as the Polyakov extended Nambu-Jona Lasinio (PNJL) and dynamical quasiparticle model with the CEP (DQPM-CP). Considering the transport coefficients and the EoS of the QGP phase, we compare our results with various results from the literature. Furthermore, out-of equilibrium properties of the QGP medium and in particular, the effect of a μB- dependence of thermodynamic and transport properties of the QGP are studied within the Parton-Hadron-String-Dynamics (PHSD) transport approach, which covers the full evolution of the system during HICs. We find that bulk observables and flow coefficients for strange hadrons as well as for antiprotons are more sensitive to the properties of the QGP, in particular to the μB - dependence of the QGP interactions.
Presolar grain isotopic ratios as constraints to nuclear physics inputs for s-process calculations
(2023)
The isotopic abundances in presolar SiC grains of AGB origin provide important and precise constraints to those star nucleosynthesis models. By comparing the values of the s-element abundances resulting from calculations with the ones measured in these dust grains, it turns out that new measurements of weak-interaction rates in ionized plasmas, as well as of neutron-capture cross sections, are needed, especially in the region near the neutron magic numbers 50 and 82.
Results on proton and Λ flow, calculated with the UrQMD model that incorporates different realistic density dependent equations of state, are presented. It is shown that the proton and hyperon flow shows sensitivity to the equation of state and especially to the appearance of a phase transition at densities below 4n0. Even though qualitatively hyperons and protons exhibit the same beam energy dependence of the flow, the quantitative results are different. In this context it is suggested that the hyperon measurements can be used to study the density dependence of the hyperon interaction in high density QCD matter.
We introduce a novel technique that utilizes a physics-driven deep learning method to reconstruct the dense matter equation of state from neutron star observables, particularly the masses and radii. The proposed framework involves two neural networks: one to optimize the EoS using Automatic Differentiation in the unsupervised learning scheme; and a pre-trained network to solve the Tolman–Oppenheimer–Volkoff (TOV) equations. The gradient-based optimization process incorporates a Bayesian picture into the proposed framework. The reconstructed EoS is proven to be consistent with the results from conventional methods. Furthermore, the resulting tidal deformation is in agreement with the limits obtained from the gravitational wave event, GW170817.
This thesis aims to investigate the properties of hadronic matter by analyzing fluctuations of conserved charges. A transport model (SMASH) is used for these studies to achieve this. The first part of this thesis focuses on examining transport coefficients, specifically the diffusion coefficients of conserved charges and the shear viscosity. The second part investigates equal-time correlations of particle numbers in the form of cumulants. The last chapter studies different aspects of the isobar collision systems Ru and Zr.
As a first step, the hadronic medium and interactions between its constituents are introduced, and simultaneously, their impact on transport coefficients is investigated. The methodology is verified by comparing the results of SMASH with Chapman-Enskog calculations, followed by examining 3-to-1 multi-particle reactions, revealing their influence on shear viscosity and electrical diffusion. The analysis of the full hadron gas considers angle-dependent cross-sections and additional elastic cross-sections via the AQM description, showing significant impacts on transport coefficients. The dependency on the number of degrees of freedom is explored, with noticeable effects on diffusion coefficients but a smaller influence on the shear viscosity. At non-zero baryon chemical potential, the diffusion coefficients are strongly influenced, while the shear viscosity remains unaffected. Overall, the study underscores the importance of individual cross-sections and the modeling of interactions on transport coefficients.
The following chapter explores fluctuations of conserved charges, crucial for understanding phase transitions in heavy-ion collision from the quark-gluon plasma to the hadronic phase. Using SMASH, the impact of global charge conservation on particle number cumulants in subvolumes of boxes simulating infinite matter is studied. Comparisons with simpler systems highlights the influence of hadronic interactions on cumulants, especially via charge annihilation processes and the results from SMASH shows agreement with analytical calculations. Calculations at finite baryon chemical potential reveals a transition from a Poisson to Skellam distribution within the net proton cumulants. It is shown that an unfolding procedure to obtain the net baryon fluctuations from the net proton ones deviates from the actual net baryon result, particularly in larger volumes. Finally, net proton correlations at vanishing baryon chemical potential align with ALICE measurements and the net proton cumulants are unaffected by deuteron formation.
In the next step, the goal is to investigate critical fluctuations in the hadronic medium. Therefore, the hadronic system is initialized with critical equilibrium fluctuations by coupling the hadron resonance gas with the 3D Ising model. The single-particle probability distributions are derived from the principle of maximum entropy. Evolving these distributions in SMASH, their development in an expanding sphere adjusted to experimental conditions can be analyzed. It reveals resonance decay and formations as the primary source that affects the particle cumulants. Because of isospin randomization processes, critical fluctuations are better preserved in net nucleon numbers. However, for the strongest coupling investigated in this work, correlations of the critical field are still present in the final state of the evolution in the net proton fluctuations. Examining cumulant dependence on rapidity windows shows a non-monotonic trend.
In the third part, collisions involving the isobars Ru and Zr are studied at a center-of-mass energy of 200 GeV. Initially, SMASH is used to study the initial conditions to hydrodynamical simulations, emphasizing the importance of the nuclear structure of isobars on the geometry of the collision area. It is found that the deformation parameters notably influence the initial state. Correlations between nucleon-nucleon pairs on eccentricity fluctuations yield no significant effect. Subsequently, the hydrodynamic model vHLLE evolves the previously explored initial conditions and for the transition between the hydrodynamic and kinetic descriptions, the Cooper-Frye formula is used. Usage of the canonical ensemble ensures the exact conservation of the conserved charges B, Q, and S. The neutron skin effect, which changes the charge distribution within Ru nuclei, is additionally considered. Fluctuations are assessed, revealing suppression in large rapidity windows due to global charge conservation. The hadronic phase modifies fluctuations of net pions, net kaons, and net protons via annihilation processes, yet fluctuations remain unaffected by the neutron skin effect.
The equation of state (EoS) of matter at extremely high temperatures and densities is currently not fully understood, and remains a major challenge in the field of nuclear physics. Neutron stars harbor such extreme conditions and therefore serve as celestial laboratories for constraining the dense matter EoS. In this thesis, we present a novel algorithm that utilizes the idea of Bayesian analysis and the computational efficiency of neural networks to reconstruct the dense matter equation of state from mass-radius observations of neutron stars. We show that the results are compatible with those from earlier works based on conventional methods, and are in agreement with the limits on tidal deformabilities obtained from the gravitational wave event, GW170817. We also observe that the resulting squared speed of sound from the reconstructed EoS features a peak, indicating a likely convergence to the conformal limit at asymptotic densities, as expected from quantum chromodynamics. The novel algorithm can also be applied across various fields faced with computational challenges in solving inverse problems. We further examine the efficiency of deep learning methods for analyzing gravitational waves from compact binary coalescences in this thesis. In particular, we develop a deep learning classifier to segregate simulated gravitational wave data into three classes: signals from binary black hole mergers, signals from binary neutron star mergers, or white noise without any signals. A second deep learning algorithm allows for the regression of chirp mass and combined tidal deformability from simulated binary neutron star mergers. An accurate estimation of these parameters is crucial to constrain the underlying EoS. Lastly, we explore the effects of finite temperatures on the binary neutron star merger remnant from GW170817. Isentropic EoSs are used to infer the frequencies of the rigidly rotating remnant and are noted to be significantly lower compared to previous estimates from zero temperature EoSs. Overall, this thesis presents novel deep learning methods to constrain the neutron star EoS, which will prove useful in future, as more observational data is expected in the upcoming years.
We have investigated the systematic differences introduced when performing a Bayesian-inference analysis of the equation of state (EOS) of neutron stars employing either variable- or constant-likelihood functions. The former has the advantage of retaining the full information on the distributions of the measurements, making exhaustive usage of the data. The latter, on the other hand, has the advantage of a much simpler implementation and reduced computational costs. In both approaches, the EOSs have identical priors and have been built using the sound speed parameterization method so as to satisfy the constraints from X-ray and gravitational waves observations, as well as those from chiral effective theory and perturbative quantum chromodynamics. In all cases, the two approaches lead to very similar results and the 90% confidence levels essentially overlap. Some differences do appear, but in regions where the probability density is extremely small and are mostly due to the sharp cutoff on the binary tidal deformability L˜ 720 set in the constant-likelihood approach. Our analysis has also produced two additional results. First, an inverse correlation between the normalized central number density, nc,TOV/ns, and the radius of a maximally massive star, RTOV. Second, and most importantly, it has confirmed the relation between the chirp mass and the binary tidal deformability. The importance of this result is that it relates chirp, which is measured very accurately, and L˜ , which contains important information on the EOS. Hence, when chirp is measured in future detections, our relation can be used to set tight constraints on L˜ .
The amplification of magnetic fields plays an important role in explaining numerous astrophysical phenomena associated with binary neutron star mergers, such as mass ejection and the powering of short gamma-ray bursts. Magnetic fields in isolated neutron stars are often assumed to be confined to a small region near the stellar surface, while they are normally taken to fill the whole star in numerical modeling of mergers. By performing high-resolution, global, and high-order general-relativistic magnetohydrodynamic simulations, we investigate the impact of a purely crustal magnetic field and contrast it with the standard configuration consisting of a dipolar magnetic field with the same magnetic energy but filling the whole star. While the crust configurations are very effective in generating strong magnetic fields during the Kelvin–Helmholtz-instability stage, they fail to achieve the same level of magnetic-field amplification of the full-star configurations. This is due to the lack of magnetized material in the neutron-star interiors to be used for further turbulent amplification and to the surface losses of highly magnetized matter in the crust configurations. Hence, the final magnetic energies in the two configurations differ by more than 1 order of magnitude. We briefly discuss the impact of these results on astrophysical observables and how they can be employed to deduce the magnetic topology in merging binaries.
Post-merger gravitational-wave signal from neutron-star binaries: a new look at an old problem
(2023)
The spectral properties of the post-merger gravitational-wave signal from a binary of neutron stars encodes a variety of information about the features of the system and of the equation of state describing matter around and above nuclear saturation density. Characterizing the properties of such a signal is an “old” problem, which first emerged when a number of frequencies were shown to be related to the properties of the binary through “quasiuniversal” relations. Here we take a new look at this old problem by computing the properties of the signal in terms of the Weyl scalar ψ4. In this way, and using a database of more than 100 simulations, we provide the first evidence for a new instantaneous frequency, y f0 4, associated with the instant of quasi-time-symmetry in the dynamics, and which also follows a quasi-universal relation. We also derive a new quasi-universal relation for the merger frequency f h mer, which provides a description of the data that is 4 times more accurate than previous expressions while requiring fewer fitting coefficients. Finally, consistent with the findings of numerous studies before ours, and using an enlarged ensemble of binary systems, we point out that the ℓ = 2, m = 1 gravitational-wave mode could become comparable with the traditional ℓ = 2, m = 2 mode on sufficiently long timescales, with strain amplitudes in a ratio |h21|/|h22| ∼ 0.1–1 under generic orientations of the binary, which could be measured by present detectors for signals with a large signal-to-noise ratio or by third-generation detectors for generic signals should no collapse occur.
Post-merger gravitational-wave signal from neutron-star binaries: a new look at an old problem
(2023)
The spectral properties of the post-merger gravitational-wave signal from a binary of neutron stars encodes a variety of information about the features of the system and of the equation of state describing matter around and above nuclear saturation density. Characterising the properties of such a signal is an “old” problem, which first emerged when a number of frequencies were shown to be related to the properties of the binary through “quasi-universal” relations. Here we take a new look at this old problem by computing the properties of the signal in terms of the Weyl scalar ψ4. In this way, and using a database of more than 100 simulations, we provide the first evidence for a new instantaneous frequency, f ψ4 0, associated with the instant of quasi timesymmetry in the postmerger dynamics, and which also follows a quasi-universal relation. We also derive a new quasi-universal relation for the merger frequency f h mer, which provides a description of the data that is four times more accurate than previous expressions while requiring fewer fitting coefficients. Finally, consistently with the findings of numerous studies before ours, and using an enlarged ensamble of binary systems we point out that the ℓ = 2, m = 1 gravitational-wave mode could become comparable with the traditional ℓ = 2, m = 2 mode on sufficiently long timescales, with strain amplitudes in a ratio |h 21|/|h 22| ∼ 0.1 − 1 under generic orientations of the binary, which could be measured by present detectors for signals with large signal-to-noise ratio or by third-generation detectors for generic signals should no collapse occur.
A considerable effort has been dedicated recently to the construction of generic equations of state (EOSs) for matter in neutron stars. The advantage of these approaches is that they can provide model-independent information on the interior structure and global properties of neutron stars. Making use of more than 106 generic EOSs, we asses the validity of quasi-universal relations of neutron star properties for a broad range of rotation rates, from slow-rotation up to the mass-shedding limit. In this way, we are able to determine with unprecedented accuracy the quasi-universal maximum-mass ratio between rotating and nonrotating stars and reveal the existence of a new relation for the surface oblateness, i.e., the ratio between the polar and equatorial proper radii. We discuss the impact that our findings have on the imminent detection of new binary neutron-star mergers and how they can be used to set new and more stringent limits on the maximum mass of nonrotating neutron stars, as well as to improve the modelling of the X-ray emission from the surface of rotating stars.
According to the inflationary theory of cosmology, most elementary particles in the current Universe were created during a period of reheating after inflation. In this Letter, we self-consistently couple the Einstein-inflaton equations to a strongly coupled quantum field theory as described by holography. We show that this leads to an inflating universe, a reheating phase, and finally a universe dominated by the quantum field theory in thermal equilibrium.
According to the inflationary theory of cosmology, most elementary particles in the current universe were created during a period of reheating after inflation. In this work we self-consistently couple the Einstein-inflaton equations to a strongly coupled quantum field theory (QFT) as described by holography. We show that this leads to an inflating universe, a reheating phase and finally a universe dominated by the QFT in thermal equilibrium.
The core of this work is represented by the investigation of the chiral phase transition, using Monte Carlo simulations and unimproved staggered fermions, both in the weak and strong coupling regimes of Quantum Chromodynamics. Based on recent results from Monte Carlo simulations, both using unimproved staggered fermions and Wilson fermions, the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf (number of flavours) in range [2:7], at zero baryon chemical potential. This achievement relies on the analytic continuation of Nf to non-integer values on the lattice, which allows to make use of extrapolation techniques to the chiral limit, where simulations are not possible. Furthermore, these results provide a resolution to the ambiguous scenario for Nf = 2 in the chiral limt. The first part of this thesis is devoted to the investigation of the chiral phase transition when a non-zero imaginary baryon chemical potential is involved, whose value corresponds to the 81% of the Roberge-Weiss one. Using the same extrapolation techniques aforementioned, the order of the chiral phase transition in the continuum and chiral limit shows compatibility with a second-order phase transition for Nf in range [2:6], highlighting a lack of dependence of the order of the chiral phase transition on the imaginary baryon chemical potential value. The second part of this thesis is about the study of the extension of the first-order chiral region in the strong coupling regime, at zero baryon chemical potential. Using Monte Carlo techniques, this can be done by investigating the Z2 boundary on a coarse lattice, whose temporal extent reads Nt = 2, and simulations are realised for Nf = 4, 8. The results in the weak coupling regime show, for $Nt = 8, 6, 4 and fixed Nf value, an inflating first-order chiral region. As in the strong coupling limit a second-order chiral phase transition is expected, the first-order chiral region has to shrink as the strong coupling regime is approached, resulting in a non-monotonic behaviour of the Z2 boundary. For Nf = 8, a critical mass on the Z2 boundary has been obtained, confirming the expected non-monotonic behaviour. For Nf = 4 the results do not provide a unique conclusion: Either a Z2 boundary at extremely low bare quark mass or a second-order chiral phase transition in the O(2) universality class in the chiral limit can take place. In addition to the two main topics, the performances of the second-order minimum norm integrator (2MN) and the fourth-order minimum norm integrator (4MN) have been compared, after implementing the 4MN one in the CL2QCD code used to realise our simulations. The 2MN integrator had already been implemented in the code since the first version was released. The two integrators belong to the class of symplectic integrators and represent an essential component of the RHMC algorithm, involved in our investigation. This step is extremely important, in order to guarantee the best quality when collecting data from simulations, and the results of the comparison suggested to favor the 2MN integrator, for both the topics.
In the framework of the LHC Injectors Upgrade Project (LIU), the CERN Proton Synchrotron Booster (PSB) went through major upgrades resulting in new effects to study, challenges to overcome and new parameter regimes to explore. To assess the achievable beam brightness limit of the machine, a series of experimental and computational studies in the transverse planes were performed. In particular, the new injection scheme induces optics perturbations that are strongly enhanced near the half-integer resonance. In this thesis, methods for dynamically measuring and correcting these perturbations and their impact on the beam performance will be presented. Additionally, the quality of the transverse beam distributions and strategies for improvement will be addressed. Finally, the space charge effects when dynamically crossing the half-integer resonance will be characterized. The results of these studies and their broader significance beyond the PSB will be discussed.
This thesis provides a detailed derivation of dissipative spin hydrodynamics from quantum field theory for systems composed of spin-0, spin-1/2, or spin-1 particles.
The Wigner function formalism is introduced for quantum fields in the respective representations of the Poincaré group, and the conserved currents, i.e., the energy-momentum tensor and the total angular momentum tensor, in various so-called pseudogauges are derived. An expansion around the semiclassical limit in powers of the Planck constant is performed.
Subsequently, kinetic equations are obtained for binary elastic scattering, using both the de Groot-van Leeuwen-van Weert and Kadanoff-Baym method, with the latter retaining the effect of quantum statistics. The resulting collision term features both local and nonlocal contributions, with the latter providing a relaxation mechanism for the spin degrees of freedom of the quasiparticles. The local-equilibrium distribution function is derived from the requirement that the local part of the collision term vanishes.
From quantum kinetic theory, dissipative spin hydrodynamics is then constructed via the method of moments, extended to particles with spin. The system of moment equations is closed via the Inverse-Reynolds Dominance (IReD) approach, resulting in a set of equations of motion describing the evolution of both ideal and dissipative degrees of freedom. The application to polarization phenomena relevant to heavy-ion collisions is discussed.
Substantial progress in the field of neuroscience has been made from anaesthetized preparations. Ketamine is one of the most used drugs in electrophysiology studies, but how ketamine affects neuronal responses is poorly understood. Here, we used in vivo electrophysiology and computational modelling to study how the auditory cortex of bats responds to vocalisations under anaesthesia and in wakefulness. In wakefulness, acoustic context increases neuronal discrimination of natural sounds. Neuron models predicted that ketamine affects the contextual discrimination of sounds regardless of the type of context heard by the animals (echolocation or communication sounds). However, empirical evidence showed that the predicted effect of ketamine occurs only if the acoustic context consists of low-pitched sounds (e.g., communication calls in bats). Using the empirical data, we updated the naïve models to show that differential effects of ketamine on cortical responses can be mediated by unbalanced changes in the firing rate of feedforward inputs to cortex, and changes in the depression of thalamo-cortical synaptic receptors. Combined, our findings obtained in vivo and in silico reveal the effects and mechanisms by which ketamine affects cortical responses to vocalisations.
Hadron lists based on experimental studies summarized by the Particle Data Group (PDG) are a crucial input for the equation of state and thermal models used in the study of strongly-interacting matter produced in heavy-ion collisions. Modeling of these strongly-interacting systems is carried out via hydrodynamical simulations, which are followed by hadronic transport codes that also require a hadronic list as input. To remain consistent throughout the different stages of modeling of a heavy-ion collision, the same hadron list with its corresponding decays must be used at each step. It has been shown that even the most uncertain states listed in the PDG from 2016 are required to reproduce partial pressures and susceptibilities from Lattice Quantum Chromodynamics with the hadronic list known as the PDG2016+. Here, we update the hadronic list for use in heavy-ion collision modeling by including the latest experimental information for all states listed in the Particle Data Booklet in 2021. We then compare our new list, called PDG2021+, to Lattice Quantum Chromodynamics results and find that it achieves even better agreement with the first principles calculations than the PDG2016+ list. Furthermore, we develop a novel scheme based on intermediate decay channels that allows for only binary decays, such that PDG2021+ will be compatible with the hadronic transport framework SMASH. Finally, we use these results to make comparisons to experimental data and discuss the impact on particle yields and spectra.
DNA binding redistributes activation domain ensemble and accessibility in pioneer factor Sox2
(2023)
More than 1600 human transcription factors orchestrate the transcriptional machinery to control gene expression and cell fate. Their function is conveyed through intrinsically disordered regions (IDRs) containing activation or repression domains but lacking quantitative structural ensemble models prevents their mechanistic decoding. Here we integrate single-molecule FRET and NMR spectroscopy with molecular simulations showing that DNA binding can lead to complex changes in the IDR ensemble and accessibility. The C-terminal IDR of pioneer factor Sox2 is highly disordered but its conformational dynamics are guided by weak and dynamic charge interactions with the folded DNA binding domain. Both DNA and nucleosome binding induce major rearrangements in the IDR ensemble without affecting DNA binding affinity. Remarkably, interdomain interactions are redistributed in complex with DNA leading to variable exposure of two activation domains critical for transcription. Charged intramolecular interactions allowing for dynamic redistributions may be common in transcription factors and necessary for sensitive tuning of structural ensembles.
We show examples of the impact of the Maxwellian averaged capture cross sections determined at n_TOF over the past 20 years on AGB stellar nucleosynthesis models. In particular, we developed an automated procedure to derive MACSs from evaluated data libraries, which are subsequently used as input to stellar models computed by means of the FuNS code. In this contribution, we present a number of s-process abundances obtained using different data libraries as input to stellar models, with a focus on the role of n_TOF data.
Subensemble is a type of statistical ensemble which is the generalization of grand canonical and canonical ensembles. The subensemble acceptance method (SAM) provides general formulas to correct the cumulants of distributions in heavy-ion collisions for the global conservation of all QCD charges. The method is applicable for an arbitrary equation of state and sufficiently large systems, such as those created in central collisions of heavy ions. The new fluctuation measures insensitive to global conservation effects are presented. The main results are illustrated in the hadron resonance gas and van der Waals fluid frameworks.
Present nuclear reaction network computations for astrophysical simulations involve many different types of rates, including neutron-capture reactions of interest for the modeling of heavy-element nucleosynthesis. While for many of them we still have to rely on theoretical calculations, an increasing number of experimentally-determined cross sections have now become available. In this contribution, we present “ASTrophysical Rate and rAw data Library” (ASTRAL), a new online database for neutron-capture cross sections based on experimental results, mainly obtained through activation and timeof-flight measurements. For the evaluation process, cross sections were re-calculated starting from raw data and by considering recent changes in physical properties of the involved isotopes (e.g., half-life and γ-ray intensities). We show the current status of the database, the techniques adopted to derive the recommended Maxwellian-averaged cross sections, and future developments.
Prediction for hyper nuclei multiplicities from GSI to LHC energies from the Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model combined with a final state coalescence approach is presented and compared to the thermal model. The influence of the coalescence radius on the collision energy and centrality dependence of the Λ3H/ΛΛ3H/Λ ratio is discussed.
This article summarizes some of the current theoretical developments and the experimental status of hypernuclei in relativistic heavy-ion collisions and elementary collisions. In particular, the most recent results of hyperhydrogen of mass A = 3 and 4 are discussed. The highlight at SQM2022 in this perspective was the discovery of the anti-hyperhydrogen-4 by the STAR Collaboration, in a large data set consisting of different collision systems. Furthermore, the production yields of hyperhydrogen-4 and hyperhelium-4 from the STAR Collaboration can be described nicely by the thermal model when the excited states of these hypernuclei are taken into account. In contrast, the production measurements in small systems (pp and p–Pb) from the ALICE Collaboration tends to favour the coalescence model over the thermal description. New measurements from STAR, ALICE and HADES Collaborations of the properties, e.g. lifetime, of A = 3 and 4 hypernuclei give similar results of these properties. Also the anti-hyperhydrogen-4 lifetime is in rather good agreement with previous measurements. Interestingly, the new STAR measurement on the R3 value, that is connected to the branching ratio, points to a Λ separation energy that is below 100 keV but definitely consistent with the value of 130 keV assumed since the 70s.
Asymptotic giant branch (AGB) stars are responsible for the production of the main component of the solar s-process distribution. Despite enormous progress in the theoretical modeling of these objects over the last few decades, many uncertainties remain. The still-unknown mechanism leading to the production of 13C neutron source is one example. The nucleosynthetic signature of AGB stars can be examined in a number of stellar sources, from spectroscopic observations of intrinsic and extrinsic stars to the heavy-element isotopic composition of presolar grains found in meteorites. The wealth of available observational data allows for constraining the processes occurring in AGB interiors. In this view, we discuss recent results from new AGB models including the effects of mixing triggered by magnetic fields, and show comparisons of the related s-process nucleosynthesis with available observations.
Lattice QCD and functional methods are making significant progress in constraining the QCD phase diagram. As an important milestone, the chiral phase transition with massless u, d-quarks at zero density is now understood to be of second order for all strange quark masses, and a smooth crossover as soon as mu,d, ≠ 0. Together with information on fluctuations and refined reweighted simulations, this bounds a possible critical point to be at µB/T ≲3. On the other hand, an approximately chiral-spin symmetric temperature window has been discovered above the chiral crossover, Tch<T ≳3Tch, with distinct correlator multiplet patterns and a pion spectral function suggesting resonance-like degrees of freedom, which dissolve graduallly with temperature.
In this work we study the 3+1-dimensional Nambu-Jona-Lasinio (NJL) model in the mean field-approximation. We carry out calculations using five different regularization schemes (two continuum and three lattice regularization schemes) with particular focus on inhomogeneous phases and condensates. The regularization schemes lead to drastically different inhomogeneous regions. We provide evidence that inhomogeneous condensates appear for all regularization schemes almost exclusively at values of the chemical potential and with wave numbers, which are of the order of or even larger than the corresponding regulators. This can be interpreted as indication that inhomogeneous phases in the 3+1-dimensional NJL model are rather artifacts of the regularization and not a consequence of the NJL Lagrangian and its symmetries.
In this work, inhomogeneous chiral phases are studied in a variety of Four-Fermion and Yukawa models in 2+1 dimensions at zero and non-zero temperature and chemical potentials. Employing the mean-field approximation, we do not find indications for an inhomogeneous phase in any of the studied models. We show that the homogeneous phases are stable against inhomogeneous perturbations. At zero temperature, full analytic results are presented.
We continue previous investigations of the (inhomogeneous) phase structure of the Gross-Neveu model in a noninteger number of spatial dimensions (1≤d<3) in the limit of an infinite number of fermion species (N→∞) at (non)zero chemical potential μ. In this work, we extend the analysis from zero to nonzero temperature T.
The phase diagram of the Gross-Neveu model in 1≤d<3 spatial dimensions is well known under the assumption of spatially homogeneous condensation with both a symmetry broken and a symmetric phase present for all spatial dimensions. In d=1 one additionally finds an inhomogeneous phase, where the order parameter, the condensate, is varying in space. Similarly, phases of spatially varying condensates are also found in the Gross-Neveu model in d=2 and d=3, as long as the theory is not fully renormalized, i.e., in the presence of a regulator. For d=2, one observes that the inhomogeneous phase vanishes, when the regulator is properly removed (which is not possible for d=3 without introducing additional parameters).
In the present work, we use the stability analysis of the symmetric phase to study the presence (for 1≤d<2) and absence (for 2≤d<3) of these inhomogeneous phases and the related moat regimes in the fully renormalized Gross-Neveu model in the μ,T-plane. We also discuss the relation between "the number of spatial dimensions" and "studying the model with a finite regulator" as well as the possible consequences for the limit d→3.
Inhomogeneous condensation in the Gross-Neveu model in noninteger spatial dimensions 1 ≤ d < 3
(2023)
The Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various 2+1-dimensional four-fermion and Yukawa models. All models are studied at non-zero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially non-uniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave function renormalization is negative, in these models is ruled out.
Inhomogeneous condensation in the Gross-Neveu model in non-integer spatial dimensions 1 ≤ d < 3
(2023)
he Gross-Neveu model in the N→∞ approximation in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to non-integer spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to non-integer spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability towards an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
We show the absence of an instability of homogeneous (chiral) condensates against spatially inhomogeneous perturbations for various (2+1)-dimensional four-fermion and Yukawa models. All models are studied at nonzero baryon chemical potential, while some of them are also subjected to chiral and isospin chemical potential. The considered theories contain up to 16 Lorentz-(pseudo)scalar fermionic interaction channels. We prove the stability of homogeneous condensates by analyzing the bosonic two-point function, which can be expressed in a purely analytical form at zero temperature. Our analysis is presented in a general manner for all of the different discussed models. We argue that the absence of an inhomogeneous chiral phase (where the chiral condensate is spatially nonuniform) follows from this lack of instability. Furthermore, the existence of a moat regime, where the bosonic wave-function renormalization is negative, in these models is ruled out.
Inhomogeneous condensation in the Gross-Neveu model in noninteger spatial dimensions 1 ≤ d < 3
(2023)
The Gross-Neveu model in the N→∞ limit in d=1 spatial dimensions exhibits a chiral inhomogeneous phase (IP), where the chiral condensate has a spatial dependence that spontaneously breaks translational invariance and the Z2 chiral symmetry. This phase is absent in d=2, while in d=3 its existence and extent strongly depends on the regularization and the value of the finite regulator. This work connects these three results smoothly by extending the analysis to noninteger spatial dimensions 1≤d<3, where the model is fully renormalizable. To this end, we adapt the stability analysis, which probes the stability of the homogeneous ground state under inhomogeneous perturbations, to noninteger spatial dimensions. We find that the IP is present for all d<2 and vanishes exactly at d=2. Moreover, we find no instability toward an IP for 2≤d<3, which suggests that the IP in d=3 is solely generated by the presence of a regulator.
In this work I investigate two different systems - spin systems and charge-density-waves. The same theoretical method is used to investigate both types of system. My investigations are motivated by experimental investigations and the goal is to describe the experimental results theoretically. For this purpose I formulate kinetic equations starting from the microscopical dynamics of the systems.
First of all, a method is formulated to derive the kinetic equations diagrammatically. Within this method an expansion in equal-time connected correlation functions is carried out. The generating functional of connected correlations is employed to derive the method.
The first system to be investigated is a thin stripe of the magnetic insulator yttrium-iron-garnet (YIG). Magnons are pumped parametrically with an external microwave field. The motivation of my theoretical investigations is to explain the experimental observations. In a small parameter range close to the confluence field strength where confluence processes of two parametrically pumped magnons with the same wave vector becomes kinematically possible the efficiency of the pumping is reduced or enhanced depending on the pumping field strength. Because it is expected that that confluence and splitting processes of magnons are essential for the experimental observations I go beyond the kinetic theories that are conventionally applied in the context of parametric excitations in YIG and investigate the influence of cubic vertices on the parametric instability of magnons in YIG.
Furthermore, the influence of phonons is investigated. Usually in the literature these are taken into account as heat bath. Here, I want to explain experiments where an accumulation of magnetoelastic bosons - magnon-phonon-quasi-particles - has been observed. I employ the method of kinetic equations to investigate this phenomenon theoretically. The kinetic theory is able to reproduce the experimental observations and it is shown that the accumulation of magnetoelastic bosons is purely incoherent.
Finally, charge-density waves (CDW) in quasi-one-dimensional materials will be investigated. Charge-density waves emerge from a Peierls-instability and are a prime example for spontaneous symmetry breaking in solids. Again, the motivation for my theoretical investigations are an experiment where the spectrum of amplitude and phase phonon modes has been measured. Starting from the Fröhlich-Hamiltonian I derive kinetic equations and from these kinetic equations the equations of motion for the CDW order parameter can be derived. The frequencies and damping rates of amplitude and phase phonon modes will be derived from the linearized equations of motion. I compare my theory with existing methods. Furthermore, I also investigate the influence of Coulomb interaction.
We use the topological heavy fermion (THF) model and its Kondo Lattice (KL) formulation to study the symmetric Kondo state in twisted bilayer graphene. Via a large-N approximation, we find a symmetric Kondo (SK) state in KL mode at fillings ν=0,±1,±2. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at ν=0,±1,±2,±3, we also find a similar symmetric state. We study the stability of the symmetric state by comparing its energy with the ordered states and find the ordered states to have lower energy. However, moving away from integer fillings by doping holes to the light bands, we find the energy difference is reduced, which suggests the loss of ordering and a tendency towards Kondo screening. In order to include many-body effects beyond the mean-field approximation, we perform dynamical mean-field theory (DMFT) calculations on the THF model. We find the spin susceptibility follows a Curie behavior at ν=0,±1,±2 down to ∼2K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At non-integer filling ν=±0.5,±0.8,±1.2 DMFT shows deviations from a 1/T-susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Finally, we study the effect of a C3z-rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C3z symmetry) can be stabilized at sufficiently large strain at ν=0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at non-integer fillings or by applying strain.
The existence of bound states induced by local impurities coupled to an insulating host depends decisively on the global topological properties of the host's electronic structure. In this context, we consider magnetic impurities modelled as classical unit-length spins that are exchange-coupled to the spinful Haldane model on the honeycomb lattice. We investigate the spectral flow of bound states with the coupling strength J in both the topologically trivial and Chern-insulating phases. In addition to conventional k-space topology, an additional, spatially local topological feature is available, based on the space of impurity-spin configurations forming, in case of R impurities, an R-fold direct product of two-dimensional spheres. Global k-space and local S-space topology are represented by different topological invariants, the first (k-space) Chern number and the R-th (S-space) spin-Chern number. We demonstrate that there is a local S-space topological transition as a function of J associated with a change in the spin Chern number and work out the implications of this for the J-dependent local electronic structure close to the impurities and, in particular, for in-gap bound states. The critical exchange couplings' dependence on the parameters of the Haldane model, and thus on the k-space topological state, is obtained numerically to construct local topological phase diagrams for systems with R=1 and R=2 impurity spins.
Topological semimetal antiferromagnets provide a rich source of exotic topological states which can be controlled by manipulating the orientation of the Néel vector, or by modulating the lattice parameters through strain. We investigate via ab initio density functional theory calculations, the effects of shear strain on the bulk and surface states n two antiferromagnetic EuCd2As2 phases with out-of-plane and in-plane spin configurations. When magnetic moments are along the c-axis, a 3% longitudinal or diagonal shear strain can tune the Dirac semimetal phase to an axion insulator phase, characterized by the parity-based invariant η4I=2. For an in-plane magnetic order, the axion insulator phase remains robust under all shear strains. We further find that for both magnetic orders, the bulk gap increases and a surface gap opens on the (001) surface up to 16 meV. Because of a nonzero η4I index and gapped states on the (001) surface, hinge modes are expected to happen on the side surface states between those gapped surface states. This result can provide a valuable insight in the realization of the long-sought axion states.
Strontium ruthenate Sr2RuO4 is an unconventional superconductor whose pairing symmetry has not been fully clarified, despite more than two decades of intensive research. Recent NMR Knight shift experiments have rekindled the Sr2RuO4 pairing debate by giving strong evidence against all odd-parity pairing states, including chiral p-wave pairing that was for a long time the leading pairing candidate. Here, we exclude additional pairing states by analyzing recent elastocaloric measurements [YS. Li et al., Nature 607, 276--280 (2022)]. To be able to explain the elastocaloric experiment, we find that unconventional even-parity pairings must include either large dx2−y2-wave or large {dxz∣dyz}-wave admixtures, where the latter possibility arises because of the body-centered point group symmetry. These {dxz∣dyz}-wave admixtures take the form of distinctively body-centered-periodic harmonics that have horizontal line nodes. Hence gxy(x2−y2)-wave and dxy-wave pairings are excluded as possible dominant even pairing states.
Strontium ruthenate Sr2RuO4 is an unconventional superconductor whose pairing symmetry has not been fully clarified, despite more than two decades of intensive research. Recent NMR Knight shift experiments have rekindled the Sr2RuO4 pairing debate by giving strong evidence against all odd-parity pairing states, including chiral p-wave pairing that was for a long time the leading pairing candidate. Here, we exclude additional pairing states by analyzing recent elastocaloric measurements [YS. Li et al., Nature 607, 276--280 (2022)]. To be able to explain the elastocaloric experiment, we find that unconventional even-parity pairings must include either large dx2−y2-wave or large {dxz∣dyz}-wave admixtures, where the latter possibility arises because of the body-centered point group symmetry. These {dxz∣dyz}-wave admixtures take the form of distinctively body-centered-periodic harmonics that have horizontal line nodes. Hence gxy(x2−y2)-wave and dxy-wave pairings are excluded as possible dominant even pairing states.
In magic angle twisted bilayer graphene, transport, thermodynamic and spectroscopic experiments pinpoint at a competition between distinct low-energy states with and without electronic order, as well as a competition between localized and delocalized charge carriers. In this study, we utilize Dynamical Mean Field Theory (DMFT) on the topological heavy Fermion (THF) model of twisted bilayer graphene to investigate the emergence of electronic correlations and long-range order in the absence of strain. We explain the nature of emergent insulating and correlated metallic states, as well as transitions between them driven by three central phenomena: (i) the formation of local spin and valley isospin moments around 100K, (ii) the ordering of the local isospin moments around 10K, and (iii) a cascadic redistribution of charge between localized and delocalized electronic states upon doping. At integer fillings, we find that low energy spectral weight is depleted in the symmetric phase, while we find insulating states with gaps enhanced by exchange coupling in the zero-strain ordered phases. Doping away from integer filling results in distinct metallic states: a "bad metal" above the ordering temperature, where coherence of the low-energy electronic excitations is suppressed by scattering off the disordered local moments, and a "good metal" in the ordered states with coherence of quasiparticles facilitated by isospin order. Upon doping, there is charge transfer between the localized and delocalized orbitals of the THF model such that they get periodically filled and emptied in between integer fillings. This charge reshuffling manifests itself in cascades of doping-induced Lifshitz transitions, local spectral weight redistributions and periodic variations of the electronic compressibility ranging from nearly incompressible to negative.
Twisted heterostructures of van der Waals materials have received much attention for their many remarkable properties. Here, we present a comprehensive theory of the long-range ordered magnetic phases of twisted bilayer α-RuCl3 via a combination of first-principles calculations and atomistic simulations. While a monolayer exhibits zigzag antiferromagnetic order with three possible ordering wave vectors, a rich phase diagram is obtained for moiré superlattices as a function of interlayer exchange and twist angle. For large twist angles, each layer spontaneously picks a single zigzag ordering wave vector, whereas, for small twist angles, the ground state involves a combination of all three wave vectors in a complex hexagonal domain structure. This multi-domain order minimizes the interlayer energy while enduring the energy cost due to the domain wall formation. Our results indicate that magnetic frustration due to stacking-dependent interlayer exchange in moiré superlattices can be used to tune the magnetic ground state and enhance quantum fluctuations in α-RuCl3.