Refine
Year of publication
Document Type
- Preprint (451)
- Article (219)
- Conference Proceeding (1)
Language
- English (671)
Has Fulltext
- yes (671)
Is part of the Bibliography
- no (671)
Keywords
Institute
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.
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.
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.
The total charm-quark production cross section per unit of rapidity dσ(cc)/dy, and the fragmentation fractions of charm quarks to different charm-hadron species f(c → hc), are measured for the first time in p–Pb collisions at √sNN = 5.02 TeV at midrapidity (−0.96 < y < 0.04 in the centre-ofmass frame) using data collected by ALICE at the CERN LHC. The results are obtained based on all the available measurements of prompt production of ground-state charm-hadron species: D0, D+,D+s, and J/ψ mesons, and Λ+cand Ξ0cbaryons. The resulting cross section is dσ(cc)/dy = 219.6±6.3 (stat.)+10.5−11.8(syst.)+7.6−2.9(extr.)±5.4 (BR)±4.6 (lumi.)±19.5 (rapidity shape) +15.0 (Ω0c) mb, which is consistent with a binary scaling of pQCD calculations from pp ollisions. The measured fragmentation fractions are compatible with those measured in pp collisions at √s = 5.02 and 13 TeV, showing an increase in the relative production rates of charm baryons with respect to charm mesons in pp and p–Pb collisions compared with e+e − and e−p collisions. The pT-integrated nuclear modification factor of charm quarks, RpPb(cc) = 0.91±0.04 (stat.) +0.08 −0.09 (syst.) +0.04 −0.03 (extr.)±0.03 (lumi.), is found to be consistent with unity and with theoretical predictions including nuclear modifications of the parton distribution functions.
This work aims to differentiate strangeness produced from hard processes (jet-like) and softer processes (underlying event) by measuring the angular correlation between a high-momentum trigger hadron (h) acting as a jet-proxy and a produced strange hadron (φ(1020) meson). Measuring h–φ correlations at midrapidity in p–Pb collisions at √sNN = 5.02 TeV as a function of event multiplicity provides insight into the microscopic origin of strangeness enhancement in small collision systems. The jet-like and the underlying-event-like strangeness production are investigated as a function of event multiplicity. They are also compared between a lower and higher momentum region. The evolution of the per-trigger yields within the near-side (aligned with the trigger hadron) and away-side (in the opposite direction of the trigger hadron) jet is studied separately, allowing for the characterization of two distinct jet-like production regimes. Furthermore, the h–φ correlations within the underlying event give access to a production regime dominated by soft production processes, which can be compared directly to the in-jet production. Comparisons between h–φ and dihadron correlations show that the observed strangeness enhancement is largely driven by the underlying event, where the φ/h ratio is significantly larger than within the jet regions. As multiplicity increases, the fraction of the total φ(1020) yield coming from jets decreases compared to the underlying event production, leading to high-multiplicity events being dominated by the increased strangeness production from the underlying event
Measurements of the production of electrons from heavy-flavour hadron decays in pp collisions at s√=13 TeV at midrapidity with the ALICE detector are presented down to a transverse momentum (pT) of 0.2 GeV/c and up to pT=35 GeV/c, which is the largest momentum range probed for inclusive electron measurements in ALICE. In p−Pb collisions, the production cross section and the nuclear modification factor of electrons from heavy-flavour hadron decays are measured in the pT range 0.5<pT<26 GeV/c at sNN−−−√=8.16 TeV. The nuclear modification factor is found to be consistent with unity within the statistical and systematic uncertainties. In both collision systems, first measurements of the yields of electrons from heavy-flavour hadron decays in different multiplicity intervals normalised to the multiplicity-integrated yield (self-normalised yield) at midrapidity are reported as a function of the self-normalised charged-particle multiplicity estimated at midrapidity. The self-normalised yields in pp and p−Pb collisions grow faster than linear with the self-normalised multiplicity. A strong pT dependence is observed in pp collisions, where the yield of high-pT electrons increases faster as a function of multiplicity than the one of low-pT electrons. The measurement in p−Pb collisions shows no pT dependence within uncertainties. The self-normalised yields in pp and p−Pb collisions are compared with measurements of other heavy-flavour, light-flavour, and strange particles, and with Monte Carlo simulations.
The production cross section of inclusive isolated photons has been measured by the ALICE experiment at the CERN LHC in pp collisions at centre-of-momentum energy of s√=13 TeV collected during the LHC Run 2 data-taking period. The measurement is performed by combining the measurements of the electromagnetic calorimeter EMCal and the central tracking detectors ITS and TPC, covering a pseudorapidity range of |ηγ|<0.67 and a transverse momentum range of 7<pγT<200 GeV/c. The result extends to lower pγT and xγT=2pγT/s√ ranges, the lowest xγT of any isolated photon measurements to date, extending significantly those measured by the ATLAS and CMS experiments towards lower pγT at the same collision energy with a small overlap between the measurements. The measurement is compared with next-to-leading order perturbative QCD calculations and the results from the ATLAS and CMS experiments as well as with measurements at other collision energies. The measurement and theory prediction are in agreement with each other within the experimental and theoretical uncertainties.
In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD)1. These partons subsequently emit further partons in a process that can be described as a parton shower2, which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass mQ and energy E, within a cone of angular size mQ/E around the emitter3. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques4,5 to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.
At particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD). The vacuum is not transparent to the partons and induces gluon radiation and quark pair production in a process that can be described as a parton shower. Studying the pattern of the parton shower is one of the key experimental tools in understanding the properties of QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass m and energy E, within a cone of angular size m/E around the emitter. A direct observation of the dead-cone effect in QCD has not been possible until now, due to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible bound hadronic states. We report the first direct observation of the QCD dead-cone by using new iterative declustering techniques to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD, which is derived more generally from its origin as a gauge quantum field theory. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.
In particle collider experiments, elementary particle interactions with large momentum transfer produce quarks and gluons (known as partons) whose evolution is governed by the strong force, as described by the theory of quantum chromodynamics (QCD). These partons subsequently emit further partons in a process that can be described as a parton shower which culminates in the formation of detectable hadrons. Studying the pattern of the parton shower is one of the key experimental tools for testing QCD. This pattern is expected to depend on the mass of the initiating parton, through a phenomenon known as the dead-cone effect, which predicts a suppression of the gluon spectrum emitted by a heavy quark of mass mQ and energy E, within a cone of angular size mQ/E around the emitter. Previously, a direct observation of the dead-cone effect in QCD had not been possible, owing to the challenge of reconstructing the cascading quarks and gluons from the experimentally accessible hadrons. We report the direct observation of the QCD dead cone by using new iterative declustering techniques to reconstruct the parton shower of charm quarks. This result confirms a fundamental feature of QCD. Furthermore, the measurement of a dead-cone angle constitutes a direct experimental observation of the non-zero mass of the charm quark, which is a fundamental constant in the standard model of particle physics.
Microscopic calculations of central collisions between heavy nuclei are used to study fragment production and the creation of collective flow. It is shown that the final phase space distributions are compatible with the expectations from a thermally equilibrated source, which in addition exhibits a collective transverse expansion. However, the microscopic analyses of the transient states in the reaction stages of highest density and during the expansion show that the system does not reach global equilibrium. Even if a considerable amount of equilibration is assumed, the connection of the measurable final state to the macroscopic parameters, e.g. the temperature, of the transient "equilibrium" state remains ambiguous.
Equilibrium properties of infinite relativistic hadron matter are investigated using the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) model. The simulations are performed in a box with periodic boundary conditions. Equilibration times depend critically on energy and baryon densities. Energy spectra of various hadronic species are shown to be isotropic and consistent with a single temperature in equilibrium. The variation of energy density versus temperature shows a Hagedorn-like behavior with a limiting temperature of 130 +/- 10 MeV. Comparison of abundances of different particle species to ideal hadron gas model predictions show good agreement only if detailed balance is implemented for all channels. At low energy densities, high mass resonances are not relevant; however, their importance raises with increasing energy density. The relevance of these different conceptual frameworks for any interpretation of experimental data is questioned.
Ratios of hadronic abundances are analyzed for pp and nucleus-nucleus collisions at sqrt(s)=20 GeV using the microscopic transport model UrQMD. Secondary interactions significantly change the primordial hadronic cocktail of the system. A comparison to data shows a strong dependence on rapidity. Without assuming thermal and chemical equilibrium, predicted hadron yields and ratios agree with many of the data, the few observed discrepancies are discussed.
Quantum Molecular Dynamics (QMD) calculations of central collisions between heavy nuclei are used to study fragment production and the creation of collective flow. It is shown that the final phase space distributions are compatible with the expectations from a thermally equilibrated source, which in addition exhibits a collective transverse expansion. However, the microscopic analyses of the transient states in the intermediate reaction stages show that the event shapes are more complex and that equilibrium is reached only in very special cases but not in event samples which cover a wide range of impact parameters as it is the case in experiments. The basic features of a new molecular dynamics model (UQMD) for heavy ion collisions from the Fermi energy regime up to the highest presently available energies are outlined.
In this paper, the concepts of microscopic transport theory are introduced and the features and shortcomings of the most commonly used ansatzes are discussed. In particular, the Ultrarelativistic Quantum Molecular Dynamics (UrQMD) transport model is described in great detail. Based on the same principles as QMD and RQMD, it incorporates a vastly extended collision term with full baryon-antibaryon symmetry, 55 baryon and 32 meson species. Isospin is explicitly treated for all hadrons. The range of applicability stretches from E lab < 100$ MeV/nucleon up to E lab> 200$ GeV/nucleon, allowing for a consistent calculation of excitation functions from the intermediate energy domain up to ultrarelativistic energies. The main physics topics under discussion are stopping, particle production and collective flow.
Basic problems of the semiclassical microscopic modelling of strongly interacting systems are discussed within the framework of Quantum Molecular Dynamics (QMD). This model allows to study the influence of several types of nucleonic interactions on a large variety of observables and phenomena occur- ring in heavy ion collisions at relativistic energies. It is shown that the same predictions can be obtained with several numerically completely di erent and independently written programs as far as the same model parameters are employed and the same basic approximations are made. Many observ- ables are robust against variations of the details of the model assumptions used. Some of the physical results, however, depend also on rather technical parameters like the preparation of the initial configuration in phase space. This crucial problem is connected with the description of the ground state of single nuclei, which di ers among the various approaches. An outlook to an improved molecular dynamics scheme for heavy ion collisions is given.
We present simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the Arctic winter 2002/2003. We integrated a Lagrangian denitrification scheme into the three-dimensional version of CLaMS that calculates the growth and sedimentation of nitric acid trihydrate (NAT) particles along individual particle trajectories. From those, we derive the HNO3 downward flux resulting from different particle nucleation assumptions. The simulation results show a clear vertical redistribution of total inorganic nitrogen (NOy), with a maximum vortex average permanent NOy removal of over 5 ppb in late December between 500 and 550 K and a corresponding increase of NOy of over 2 ppb below about 450 K. The simulated vertical redistribution of NOy is compared with balloon observations by MkIV and in-situ observations from the high altitude aircraft Geophysica. Assuming a globally uniform NAT particle nucleation rate of 3.4·10−6 cm−3 h−1 in the model, the observed denitrification is well reproduced. In the investigated winter 2002/2003, the denitrification has only moderate impact (<=10%) on the simulated vortex average ozone loss of about 1.1 ppm near the 460 K level. At higher altitudes, above 600 K potential temperature, the simulations show significant ozone depletion through NOx-catalytic cycles due to the unusual early exposure of vortex air to sunlight.
We present simulations with the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the Arctic winter 2002/2003. We integrated a Lagrangian denitrification scheme into the three-dimensional version of CLaMS that calculates the growth and sedimentation of nitric acid trihydrate (NAT) particles along individual particle trajectories. From those, we derive the HNO3 downward flux resulting from different particle nucleation assumptions. The simulation results show a clear vertical redistribution of total inorganic nitrogen ( ), with a maximum vortex average permanent removal of over 5ppb in late December between 500 and 550K and a corresponding increase of of over 2ppb below about 450K. The simulated vertical redistribution of is compared with balloon observations by MkIV and in-situ observations from the high altitude aircraft Geophysica. Assuming a globally uniform NAT particle nucleation rate of 7.8x10-6cm-3h-1 in the model, the observed denitrification is well reproduced.
In the investigated winter 2002/2003, the denitrification has only moderate impact (≤14%) on the simulated vortex average ozone loss of about 1.1ppm near the 460K level. At higher altitudes, above 600K potential temperature, the simulations show significant ozone depletion through -catalytic cycles due to the unusual early exposure of vortex air to sunlight.
The Match method for the quantification of polar chemical ozone loss is investigated mainly with respect to the impact of the transport of air masses across the vortex edge. For the winter 2002/03, we show that significant transport across the vortex edge occurred and was simulated by the Chemical Lagrangian Model of the Stratosphere. In-situ observations of inert tracers and ozone from HAGAR on the Geophysica aircraft and balloon-borne sondes, and remote observations from MIPAS on the ENVISAT satellite were reproduced well by CLaMS. The model even reproduced a small vortex remnant that remained a distinct feature until June 2003 and was also observed in-situ by a balloon-borne whole air sampler. We use this CLaMS simulation to quantify the impact of transport across the vortex edge on ozone loss estimates from the Match method. We show that a time integration of the determined vortex average ozone loss rates, as performed in Match, results in a larger ozone loss than the polar vortex average ozone loss in CLaMS. The determination of the Match ozone loss rates is also influenced by the transport of air across the vortex edge. We use the model to investigate how the sampling of the ozone sondes on which Match is based represents the vortex average ozone loss rate. Both the time integration of ozone loss and the determination of ozone loss rates for Match are evaluated using the winter 2002/2003 CLaMS simulation. These impacts can explain the majority of the differences between CLaMS and Match column ozone loss. While the investigated effects somewhat reduce the apparent discrepancy in January ozone loss rates reported earlier, a distinct discrepancy between simulations and Match remains. However, its contribution to the accumulated ozone loss over the winter is not large.
The Match method for quantification of polar chemical ozone loss is investigated mainly with respect to the impact of mixing across the vortex edge onto this estimate. We show for the winter 2002/03 that significant mixing across the vortex edge occurred and was accurately modeled by the Chemical Lagrangian Model of the Stratosphere. Observations of inert tracers and ozone in-situ from HAGAR on the Geophysica aircraft and sondes and also remote from MIPAS on ENVISAT were reproduced well. The model even reproduced a small vortex remnant that was isolated until June 2003 and was observed in-situ by a balloon-borne whole air sampler. We use this CLaMS simulation to quantify the impact of cross vortex edge mixing on the results of the Match method. It is shown that a time integration of the determined vortex average ozone loss rates as performed in Match results in larger ozone loss than the polar vortex average ozone loss in CLaMS. Also, the determination of the Match ozone loss rates can be influenced by mixing. This is especially important below 430 K, where ozone outside the vortex is lower than inside and the vortex boundary is not a strong transport barrier. This effect and further sampling effects cause an offset between vortex average ozone loss rates derived from Match and deduced from CLaMS with an even sampling for the entire vortex. Both, the time-integration of ozone loss and the determination of ozone loss rates for Match are evaluated using the winter 2002/03 CLaMS simulation. These impacts can explain the differences between CLaMS and Match column ozone loss. While the investigated effects somewhat reduce the apparent discrepancy in January ozone loss rates, a discrepancy between simulations and Match remains. However, its contribution to the accumulated ozone loss over the winter is not large.
The behavior of hadronic matter at high baryon densities is studied within Ultrarelativistic Quantum Molecular Dynamics (URQMD). Baryonic stopping is observed for Au+Au collisions from SIS up to SPS energies. The excitation function of flow shows strong sensitivities to the underlying equation of state (EOS), allowing for systematic studies of the EOS. Dilepton spectra are calculated with and without shifting the rho pole. Except for S+Au collisions our calculations reproduce the CERES data.
The behavior of hadronic matter at high baryon densities is studied within Ultrarelativistic Quantum Molecular Dynamics (URQMD). Baryonic stopping is observed for Au+Au collisions from SIS up to SPS energies. The excitation function of flow shows strong sensitivities to the underlying equation of state (EOS), allowing for systematic studies of the EOS. Effects of a density dependent pole of the rho-meson propagator on dilepton spectra are studied for different systems and centralities at CERN energies.
We study the thermodynamic properties of infinite nuclear matter with the Ultrarelativistic Quantum Molecular Dynamics (URQMD), a semiclassical transport model, running in a box with periodic boundary conditions. It appears that the energy density rises faster than T4 at high temperatures of T approx. 200 - 300 MeV. This indicates an increase in the number of degrees of freedom. Moreover, We have calculated direct photon production in Pb+Pb collisions at 160 GeV/u within this model. The direct photon slope from the microscopic calculation equals that from a hydrodynamical calculation without a phase transition in the equation of state of the photon source.
Nonequilibrium models (three-fluid hydrodynamics, UrQMD, and quark molecular dynamics) are used to discuss the uniqueness of often proposed experimental signatures for quark matter formation in relativistic heavy ion collisions from the SPS via RHIC to LHC. It is demonstrated that these models - although they do treat the most interesting early phase of the collisions quite differently (thermalizing QGP vs. coherent color fields with virtual particles) -- all yield a reasonable agreement with a large variety of the available heavy ion data. Hadron/hyperon yields, including J/Psi meson production/suppression, strange matter formation, dileptons, and directed flow (bounce-off and squeeze-out) are investigated. Observations of interesting phenomena in dense matter are reported. However, we emphasize the need for systematic future measurements to search for simultaneous irregularities in the excitation functions of several observables in order to come close to pinning the properties of hot, dense QCD matter from data. The role of future experiments with the STAR and ALICE detectors is pointed out.
The microscopic phasespace approach URQMD is used to investigate the stopping power and particle production in heavy systems at SPS and RHIC energies. We find no gap in the baryon rapidity distribution even at RHIC. For CERN energies URQMD shows a pile up of baryons and a supression of multi-nucleon clusters at midrapidity.
The quantum statistical model (QSM) is used to calculate nuclear fragment distributions in chemical equilibrium. Several observable isotopic effects are predicted for intermediate energy heavy ion collisions. It is demonstrated that particle ratios for different systemsdo not depend on the breakup density-the only free parameter in our model.The importance of entropy measurements is discussed. Specific particle ratios for the system Au-Au are predicted, which can be used to determine the chemical potentials of the hot midrapidity fragment source in nearly central heavy ion collisions. Pacs-Nr. 25.70 Pq
Spectra of various particle species have been calculated with the Quantum Molecular Dynamics (QMD) model for very central collisions of Au+Au. They are compatible with the idea of a fully stopped thermal source which exhibits a transversal expansion besides the thermal distribution of an ideal gas. How- ever, the microscopic analyses of the local flow velocities and temperatures indicate much lower temperatures at densities associated with the freeze-out. The results express the overall impossibility of a model-independent determi- nation of nuclear temperatures from heavy ion spectral data, also at other energies (e.g. CERN) or for other species (i.e. pions, kaons, hyperons)
A quasiclassical Pauli potential is used to simulate the Fermi motion of nucleons in a molecular dynamical simulation of heavy ion collisions. The thermostatic properties of a Fermi gas with and without interactions are presented. The inclusion of this Pauli potential into the quantum molecular dynamics (QMD) approach yields a model with well defined fermionic ground states, which is therefore also able to give the excitation energies of the emitted fragments. The deexcitation mechanisms (particle evaporation and multifragmentation) of the new model are investigated. The dynamics of the QMD with Pauli potential is tested by a wide range of comparisons of calculated and experimental double-differential cross sections for inclusive p-induced reactions at incident energies of 80 to 160 MeV. Results at 256 and 800 MeV incident proton energy are presented as predictions for completed experiments which are as yet unpublished.
This paper investigates the global stratospheric Brewer–Dobson circulation (BDC) in the ERA5 meteorological reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF). The analysis is based on simulations of stratospheric mean age of air, including the full age spectrum, with the Lagrangian transport model CLaMS (Chemical Lagrangian Model of the Stratosphere), driven by reanalysis winds and total diabatic heating rates. ERA5-based results are compared to results based on the preceding ERA-Interim reanalysis. Our results show a significantly slower BDC for ERA5 than for ERA-Interim, manifesting in weaker diabatic heating rates and higher age of air. In the tropical lower stratosphere, heating rates are 30 %–40 % weaker in ERA5, likely correcting a bias in ERA-Interim. At 20 km and in the Northern Hemisphere (NH) stratosphere, ERA5 age values are around the upper margin of the uncertainty range from historical tracer observations, indicating a somewhat slow–biased BDC. The age trend in ERA5 over the 1989–2018 period is negative throughout the stratosphere, as climate models predict in response to global warming. However, the age decrease is not linear but steplike, potentially caused by multi-annual variability or changes in the observations included in the assimilation. During the 2002–2012 period, the ERA5 age shows a similar hemispheric dipole trend pattern as ERA-Interim, with age increasing in the NH and decreasing in the Southern Hemisphere (SH). Shifts in the age spectrum peak and residual circulation transit times indicate that reanalysis differences in age are likely caused by differences in the residual circulation. In particular, the shallow BDC branch accelerates in both reanalyses, whereas the deep branch accelerates in ERA5 and decelerates in ERA-Interim.