Refine
Year of publication
Document Type
- Article (23)
- In Proceedings (12)
- ZIB-Report (2)
Language
- English (37)
Is part of the Bibliography
- no (37)
The coma of comet 67P/Churyumov-Gerasimenko has been probed by the Rosetta spacecraft and shows a variety of different molecules. The ROSINA COmet Pressure Sensor and the Double Focusing Mass Spectrometer provide in-situ densities for many volatile compounds including the 14 gas species H2O, CO2, CO, H2S, O2, C2H6, CH3OH, H2CO, CH4, NH3, HCN, C2H5OH, OCS, and CS2. We fit the observed densities during the entire comet mission between August 2014 and September 2016 to an inverse coma model. We retrieve surface emissions on a cometary shape with 3996 triangular elements for 50 separated time intervals. For each gas we derive systematic error bounds and report the temporal evolution of the production, peak production, and the time-integrated total production. We discuss the production for the two lobes of the nucleus and for the northern and southern hemispheres. Moreover we provide a comparison of the gas production with the seasonal illumination.
Computing the Hierarchical Equations of Motion (HEOM) is by itself a challenging problem, and so is writing portable production code that runs efficiently on a variety of architectures while scaling from PCs to supercomputers. We combined both challenges to push the boundaries of simulating quantum systems, and to evaluate and improve methodologies for scientific software engineering.
Our contributions are threefold: We present the first distributed memory implementation of the HEOM method (DM-HEOM), we describe an interdisciplinary development workflow, and we provide guidelines and experiences for designing distributed, performance-portable HPC applications with MPI-3, OpenCL and other state-of-the-art programming models. We evaluated the resulting code on multi- and many-core CPUs as well as GPUs, and demonstrate scalability on a Cray XC40 supercomputer for the PS I molecular light harvesting complex.
Context. The change of the rotation period and the orientation of the rotation axis of comet 67P/Churyumov-Gerasimenko (67P/C-G) is deducible from images taken by the scientific imaging instruments on-board the Rosetta mission with high precision. Non gravitational forces are a natural explanation for these data.
Aims. We describe observed changes for the orientation of the rotation axis and the rotation period of 67P/C-G. For these changes we give an explanation based on a sublimation model with a best-fit for the surface active fraction (model P). Torque effects of periodically changing gas emissions on the surface are considered.
Methods. We solve the equation of state for the angular momentum in the inertial and the body- fixed frames and provide an analytic theory of the rotation changes in terms of Fourier coefficients, generally applicable to periodically forced rigid body dynamics.
Results. The torque induced changes of the rotation state constrain the physical properties of the surface, the sublimation rate and the local active fraction of the surface.
Conclusions. We determine a distribution of the local surface active fraction in agreement with the rotation properties, period and orientation, of 67P/C-G. The torque movement confirms that the sublimation increases faster than the insolation towards perihelion. The derived relatively uniform activity pattern is discussed in terms of related surface features.
Cometary activity affects the orbital motion and rotation state due to sublimation induced forces. The availability of precise rotation-axis orientation and position data from the Rosetta mission allows one to accurately determine the outgassing of comet Churyumov-Gerasimenko/67P (67P).
We derive the observed non-gravitational acceleration of 67P directly from the Rosetta spacecraft trajectory. From the non-gravitational acceleration we recover the diurnal outgassing variations and study a possible delay of the sublimation response with respect to the peak solar illumination. This allows us to compare the non-gravitational acceleration of 67P with expectations based on empirical models and common assumptions about the sublimation process.
We use an iterative orbit refinement and Fourier decomposition of the diurnal activity to derive the outgassing induced non-gravitational acceleration. The uncertainties of the data reduction are established by a sensitivity analysis of an ensemble of best-fit orbits for comet 67P.
We find that the Marsden non-gravitational acceleration parameters reproduce part of the non-gravitational acceleration but need to be augmented by an analysis of the nucleus geometry and surface illumination to draw conclusions about the sublimation process on the surface. The non-gravitational acceleration follows closely the subsolar latitude (seasonal illumination), with a small lag angle with respect to local noon around perihelion. The observed minor changes of the rotation axis do not favor forced precession models for the non-gravitational acceleration.
In contrast to the sublimation induced torques, the non-gravitational acceleration does not put strong constraints on localized active areas on the nucleus. We find a close agreement of the orbit deduced non-gravitational acceleration and the water production independently derived from Rosetta in-situ measurement.
Time- and frequency resolved optical signals provide insights into the properties of light harvesting molecular complexes, including excitation energies, dipole strengths and orientations, as well as in the exciton energy flow through the complex. The hierarchical equations of motion (HEOM) provide a unifying theory, which allows one to study the combined effects of system-environment dissipation and non-Markovian memory without making restrictive assumptions about weak or strong couplings or separability of vibrational and electronic degrees of freedom. With increasing system size the exact solution of the open quantum system dynamics requires memory and compute resources beyond a single compute node. To overcome this barrier, we developed a scalable variant of HEOM. Our distributed memory HEOM, DM-HEOM, is a universal tool for open quantum system dynamics. It is used to accurately compute all experimentally accessible time- and frequency resolved processes in light harvesting molecular complexes with arbitrary system-environment couplings for a wide range of temperatures and complex sizes.
Energy flow in the Photosystem I supercomplex: comparison of approximative theories with DM-HEOM
(2018)
We analyze the exciton dynamics in PhotosystemI from Thermosynechococcus elongatus using the distributed memory implementation of the hierarchical equation of motion (DM-HEOM) for the 96 Chlorophylls in the monomeric unit. The exciton-system parameters are taken from a first principles calculation. A comparison of the exact results with Foerster rates and Markovian approximations allows one to validate the exciton transfer times within the complex and to identify deviations from approximative theories. We show the optical absorption, linear, and circular dichroism spectra obtained with DM-HEOM and compare them to experimental results.
Surface localization of gas sources on comet 67P/Churyumov-Gerasimenko based on DFMS/COPS data
(2019)
We reconstruct the temporal evolution of the source distribution for the four major gas species H2O, CO2, CO, and O2 on the surface of comet 67P/Churyumov-Gerasimenko during its 2015 apparition. The analysis applies an inverse coma model and fits to data between August 6th 2014 and September 5th 2016 measured with the Double Focusing Mass Spectrometer (DFMS) of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) and the COmet Pressure Sensor (COPS). The spatial distribution of gas sources with their temporal variation allows one to construct surface maps for gas emissions and to evaluate integrated productions rates. For all species peak production rates and integrated productions rates per orbit are evaluated separately for the northern and the southern hemisphere. The nine most active emitting areas on the comet’s surface are defined and their correlation to emissions for each of the species is discussed.
Matching the activity of comet 67P/Churyumov-Gerasimenko with long-term ground-based astrometry
(2021)
50 years of astrometric data for comet 67P/C-G (orbital period about 6.45 years) provides a unique opportunity to benchmark non-gravitational acceleration models to the in situ measurements of the volatile release performed from the Rosetta rendezvous mission (2014-2016). Taken together, the Earth-bound and in-situ data yields lower fit errors and serves as a test-case for our ability to deduce thermophysical quantities of cometary nuclei from the Earth-bound observations.
We reconstruct the temporal evolution of surface emissions for the four major gas species H2O, CO2, CO, and O2 emitted during the 2015 apparition of comet 67P/Churyumov-Gerasimenko (67P/C-G). Measurements from the Double Focusing Mass Spectrometer (DFMS) of the Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) and the COmet Pressure Sensor (COPS) are used to determine the gas sources on the surface with an inverse gas model for the entire coma. For all species, peak production rates and integrated production rates per orbit are evaluated separately for the northern and the southern hemisphere. Complemented with the total mass production, this allows us to estimate the dust-to-gas ratio of the
emitted material.
Based on about 1 million of pressure measurements around comet 67P/Churyumov-Gerasimenko we reconstruct the gas emission across the entire nucleus. Dust particles are seeded in the gas model and the resulting dust distribution follows a daily pattern which agrees with observations if a uniform dust release across the entire sunlit surface is assumed.