530 Physik
Refine
Year of publication
Document Type
- Article (868)
- Doctoral Thesis (361)
- Postprint (95)
- Other (47)
- Preprint (36)
- Habilitation Thesis (22)
- Review (11)
- Master's Thesis (10)
- Monograph/Edited Volume (4)
- Course Material (2)
Keywords
- diffusion (43)
- anomalous diffusion (31)
- gamma rays: general (20)
- synchronization (18)
- Synchronisation (16)
- organic solar cells (16)
- stochastic processes (14)
- ISM: supernova remnants (13)
- cosmic rays (13)
- data analysis (12)
Institute
- Institut für Physik und Astronomie (1458) (remove)
We report on the adsorption kinetics of azoben-zene-containing surfactants on solid surfaces of different hydrophobicity. The understanding of this processes is of great importance for many interfacial phenomena that can be actuated and triggered by light, since the surfactant molecules contain a photoresponsive azobenzene group in their hydrophobic tail. Three surfactant types are studied, differing in the spacer connecting the headgroup and the azobenzene unit by between 6 and 10 CH2 groups. Under irradiation with light of a suitable wavelength, the azobenzene undergoes reversible photoisomerization between two states, a nonpolar trans-state and a highly polar cis-state. Consequently, the surfactant molecule changes its hydrophobicity and thus affinity to a surface depending on the photoisomerization state of the azobenzene. The adsorption behavior on hydrophilic (glass) and hydrophobic (TeflonAF) surfaces is analyzed using quartz crystal microbalance with dissipation (QCM-D) and zeta-potential measurements. At equilibrium, the adsorbed surfactant amount is almost twice as large on glass compared to TeflonAF for both isomers. The adsorption rate for the trans-isomers on both surfaces is similar, but the desorption rate of the trans-isomers is faster at the glass-water interface than at the Teflon-water interface. This result demonstrates that the trans-isomers have higher affinity for the glass surface, so the trans-to-cis ratios on glass and TeflonAF are 80/1 and 2/1, respectively, with similar trends for all three surfactant types.
Chemotactic motion in a chemical gradient is an essential cellular function that controls many processes in the living world. For a better understanding and more detailed modelling of the underlying mechanisms of chemotaxis, quantitative investigations in controlled environments are needed. We developed a setup that allows us to separately address the dependencies of the chemotactic motion on the average background concentration and on the gradient steepness of the chemoattractant. In particular, both the background concentration and the gradient steepness can be kept constant at the position of the cell while it moves along in the gradient direction. This is achieved by generating a well-defined chemoattractant gradient using flow photolysis. In this approach, the chemoattractant is released by a light-induced reaction from a caged precursor in a microfluidic flow chamber upstream of the cell. The flow photolysis approach is combined with an automated real-time cell tracker that determines changes in the cell position and triggers movement of the microscope stage such that the cell motion is compensated and the cell remains at the same position in the gradient profile. The gradient profile can be either determined experimentally using a caged fluorescent dye or may be alternatively determined by numerical solutions of the corresponding physical model. To demonstrate the function of this adaptive microfluidic gradient generator, we compare the chemotactic motion of Dictyostelium discoideum cells in a static gradient and in a gradient that adapts to the position of the moving cell. Published by AIP Publishing.
Earth’s surface temperature will continue to rise for another 20 to 30 years even with the strongest carbon emission reduction currently considered. The associated changes in rainfall patterns can result in an increased flood risk worldwide. We compute the required increase in flood protection to keep high-end fluvial flood risk at present levels. The analysis is carried out worldwide for subnational administrative units. Most of the United States, Central Europe, and Northeast and West Africa, as well as large parts of India and Indonesia, require the strongest adaptation effort. More than half of the United States needs to at least double their protection within the next two decades. Thus, the need for adaptation to increased river flood is a global problem affecting industrialized regions as much as developing countries.
Gravitational-wave (GW) astrophysics is a field in full blossom. Since the landmark detection of GWs from a binary black hole on September 14th 2015, fifty-two compact-object binaries have been reported by the LIGO-Virgo collaboration. Such events carry astrophysical and cosmological information ranging from an understanding of how black holes and neutron stars are formed, what neutron stars are composed of, how the Universe expands, and allow testing general relativity in the highly-dynamical strong-field regime. It is the goal of GW astrophysics to extract such information as accurately as possible. Yet, this is only possible if the tools and technology used to detect and analyze GWs are advanced enough. A key aspect of GW searches are waveform models, which encapsulate our best predictions for the gravitational radiation under a certain set of parameters, and that need to be cross-correlated with data to extract GW signals. Waveforms must be very accurate to avoid missing important physics in the data, which might be the key to answer the fundamental questions of GW astrophysics. The continuous improvements of the current LIGO-Virgo detectors, the development of next-generation ground-based detectors such as the Einstein Telescope or the Cosmic Explorer, as well as the development of the Laser Interferometer Space Antenna (LISA), demand accurate waveform models. While available models are enough to capture the low spins, comparable-mass binaries routinely detected in LIGO-Virgo searches, those for sources from both current and next-generation ground-based and spaceborne detectors must be accurate enough to detect binaries with large spins and asymmetry in the masses. Moreover, the thousands of sources that we expect to detect with future detectors demand accurate waveforms to mitigate biases in the estimation of signals’ parameters due to the presence of a foreground of many sources that overlap in the frequency band. This is recognized as one of the biggest challenges for the analysis of future-detectors’ data, since biases might hinder the extraction of important astrophysical and cosmological information from future detectors’ data. In the first part of this thesis, we discuss how to improve waveform models for binaries with high spins and asymmetry in the masses. In the second, we present the first generic metrics that have been proposed to predict biases in the presence of a foreground of many overlapping signals in GW data.
For the first task, we will focus on several classes of analytical techniques. Current models for LIGO and Virgo studies are based on the post-Newtonian (PN, weak-field, small velocities) approximation that is most natural for the bound orbits that are routinely detected in GW searches. However, two other approximations have risen in prominence, the post-Minkowskian (PM, weak- field only) approximation natural for unbound (scattering) orbits and the small-mass-ratio (SMR) approximation typical of binaries in which the mass of one body is much bigger than the other. These are most appropriate to binaries with high asymmetry in the masses that challenge current waveform models. Moreover, they allow one to “cover” regions of the parameter space of coalescing binaries, thereby improving the interpolation (and faithfulness) of waveform models. The analytical approximations to the relativistic two-body problem can synergically be included within the effective-one-body (EOB) formalism, in which the two-body information from each approximation can be recast into an effective problem of a mass orbiting a deformed Schwarzschild (or Kerr) black hole. The hope is that the resultant models can cover both the low-spin comparable-mass binaries that are routinely detected, and the ones that challenge current models. The first part of this thesis is dedicated to a study about how to best incorporate information from the PN, PM, SMR and EOB approaches in a synergistic way. We also discuss how accurate the resulting waveforms are, as compared against numerical-relativity (NR) simulations. We begin by comparing PM models, whether alone or recast in the EOB framework, against PN models and NR simulations. We will show that PM information has the potential to improve currently-employed models for LIGO and Virgo, especially if recast within the EOB formalism. This is very important, as the PM approximation comes with a host of new computational techniques from particle physics to exploit. Then, we show how a combination of PM and SMR approximations can be employed to access previously-unknown PN orders, deriving the third subleading PN dynamics for spin-orbit and (aligned) spin1-spin2 couplings. Such new results can then be included in the EOB models currently used in GW searches and parameter estimation studies, thereby improving them when the binaries have high spins. Finally, we build an EOB model for quasi-circular nonspinning binaries based on the SMR approximation (rather than the PN one as usually done). We show how this is done in detail without incurring in the divergences that had affected previous attempts, and compare the resultant model against NR simulations. We find that the SMR approximation is an excellent approximation for all (quasi-circular nonspinning) binaries, including both the equal-mass binaries that are routinely detected in GW searches and the ones with highly asymmetric masses. In particular, the SMR-based models compare much better than the PN models, suggesting that SMR-informed EOB models might be the key to model binaries in the future. In the second task of this thesis, we work within the linear-signal ap- proximation and describe generic metrics to predict inference biases on the parameters of a GW source of interest in the presence of confusion noise from unfitted foregrounds and from residuals of other signals that have been incorrectly fitted out. We illustrate the formalism with simple (yet realistic) LISA sources, and demonstrate its validity against Monte-Carlo simulations. The metrics we describe pave the way for more realistic studies to quantify the biases with future ground-based and spaceborne detectors.
This Thesis puts its focus on the physics of neutron stars and its description with methods of numerical relativity. In the first step, a new numerical framework the Whisky2D code will be developed, which solves the relativistic equations of hydrodynamics in axisymmetry. Therefore we consider an improved formulation of the conserved form of these equations. The second part will use the new code to investigate the critical behaviour of two colliding neutron stars. Considering the analogy to phase transitions in statistical physics, we will investigate the evolution of the entropy of the neutron stars during the whole process. A better understanding of the evolution of thermodynamical quantities, like the entropy in critical process, should provide deeper understanding of thermodynamics in relativity. More specifically, we have written the Whisky2D code, which solves the general-relativistic hydrodynamics equations in a flux-conservative form and in cylindrical coordinates. This of course brings in 1/r singular terms, where r is the radial cylindrical coordinate, which must be dealt with appropriately. In the above-referenced works, the flux operator is expanded and the 1/r terms, not containing derivatives, are moved to the right-hand-side of the equation (the source term), so that the left hand side assumes a form identical to the one of the three-dimensional (3D) Cartesian formulation. We call this the standard formulation. Another possibility is not to split the flux operator and to redefine the conserved variables, via a multiplication by r. We call this the new formulation. The new equations are solved with the same methods as in the Cartesian case. From a mathematical point of view, one would not expect differences between the two ways of writing the differential operator, but, of course, a difference is present at the numerical level. Our tests show that the new formulation yields results with a global truncation error which is one or more orders of magnitude smaller than those of alternative and commonly used formulations. The second part of the Thesis uses the new code for investigations of critical phenomena in general relativity. In particular, we consider the head-on-collision of two neutron stars in a region of the parameter space where two final states a new stable neutron star or a black hole, lay close to each other. In 1993, Choptuik considered one-parameter families of solutions, S[P], of the Einstein-Klein-Gordon equations for a massless scalar field in spherical symmetry, such that for every P > P⋆, S[P] contains a black hole and for every P < P⋆, S[P] is a solution not containing singularities. He studied numerically the behavior of S[P] as P → P⋆ and found that the critical solution, S[P⋆], is universal, in the sense that it is approached by all nearly-critical solutions regardless of the particular family of initial data considered. All these phenomena have the common property that, as P approaches P⋆, S[P] approaches a universal solution S[P⋆] and that all the physical quantities of S[P] depend only on |P − P⋆|. The first study of critical phenomena concerning the head-on collision of NSs was carried out by Jin and Suen in 2007. In particular, they considered a series of families of equal-mass NSs, modeled with an ideal-gas EOS, boosted towards each other and varied the mass of the stars, their separation, velocity and the polytropic index in the EOS. In this way they could observe a critical phenomenon of type I near the threshold of black-hole formation, with the putative solution being a nonlinearly oscillating star. In a successive work, they performed similar simulations but considering the head-on collision of Gaussian distributions of matter. Also in this case they found the appearance of type-I critical behaviour, but also performed a perturbative analysis of the initial distributions of matter and of the merged object. Because of the considerable difference found in the eigenfrequencies in the two cases, they concluded that the critical solution does not represent a system near equilibrium and in particular not a perturbed Tolmann-Oppenheimer-Volkoff (TOV) solution. In this Thesis we study the dynamics of the head-on collision of two equal-mass NSs using a setup which is as similar as possible to the one considered above. While we confirm that the merged object exhibits a type-I critical behaviour, we also argue against the conclusion that the critical solution cannot be described in terms of equilibrium solution. Indeed, we show that, in analogy with what is found in, the critical solution is effectively a perturbed unstable solution of the TOV equations. Our analysis also considers fine-structure of the scaling relation of type-I critical phenomena and we show that it exhibits oscillations in a similar way to the one studied in the context of scalar-field critical collapse.
The unceasing impact of intense sunlight on earth constitutes a continuous source of energy fueling countless natural processes. On a molecular level, the energy contained in the electromagnetic radiation is transferred through photochemical processes into chemical or thermal energy. In the course of such processes, photo-excitations promote molecules into thermally inaccessible excited states. This induces adaptations of their molecular geometry according to the properties of the excited state. Decay processes towards energetically lower lying states in transient molecular geometries result in the formation of excited state relaxation pathways. The photo-chemical relaxation mechanisms depend on the studied system itself, the interactions with its chemical environment and the character of the involved states. This thesis focuses on systems in which photo-induced deprotonation processes occur at specific atomic sites.
To detect these excited-state proton dynamics at the affected atoms, a local probe of molecular electronic structure is required. Therefore, site-selective and orbital-specific K-edge soft X-ray spectroscopy techniques are used here to detect photo-induced proton dynamics in gaseous and liquid sample environments. The protonation of nitrogen (N) sites in organic molecules and the oxygen (O) atom in the water molecule are probed locally through transitions between 1s orbitals and the p-derived molecular valence electronic structure. The used techniques are X-ray absorption spectroscopy (XAS) and resonant inelastic X-ray scattering (RIXS). Both yield access to the unoccupied local valence electronic structure, whereas the latter additionally probes occupied states.
We apply these probes in optical pump X-ray probe experiments to investigate valence excited-state proton transfer capabilities of aqueous 2-thiopyridone. A characteristic shift of N K-edge X-ray absorption resonances as well as a distinct X-ray emission line are established by us as spectral fingerprints of N deprotonation in the system. We utilize them to identify photo-induced N deprotonation of 2-thiopyridone on femtosecond timescales, in optical pump N K-edge RIXS probe measurements. We further establish excited state proton transfer mechanisms on picosecond and nanosecond timescales along the dominant relaxation pathways of 2-thiopyridone using transient N K-edge XAS.
Despite being an excellent probe mechanism for valence excited-state proton dynamics, the K-edge core-excitation itself also disturbs the electronic structure at specific sites of a molecule. The rapid reaction of protons to 1s photo-excitations can yield directional structural distortions within the femtosecond core-excited state lifetime. These directional proton dynamics can change the energetic separation of eigenstates of the system and alter probabilities for radiative decay between them. Both effects yield spectral signatures of the dynamics in RIXS spectra.
Using these signatures of RIXS transitions into electronically excited states, we investigate proton dynamics induced by N K-edge excitation in the amino-acid histidine. The minor core-excited state dynamics of histidine in basic and neutral chemical environments allow us to establish XAS and RIXS spectral signatures of different N protonation states at its imidazole N sites. Based on these signatures, we identify an excitation-site-independent N-H dissociation for N K-edge excitation under acidic conditions.
Such directional structural deformations, induced by core-excitations, also make proton dynamics in electronic ground states accessible through RIXS transitions into vibrationally excited states. In that context, we interpret high resolution RIXS spectra of the water molecule for three O K-edge resonances based on quantum-chemical wave packet propagation simulations. We show that highly oriented ground state vibrational modes of coupled nuclear motion can be populated through RIXS processes by preparation of core-excited state nuclear wave packets with the same directionality. Based on that, we analytically derive the possibility to extract one-dimensional directional cuts through potential energy surfaces of molecular systems from the corresponding RIXS spectra. We further verify this concept through the extraction of the gas-phase water ground state potential along three coordinates from experimental data in comparison to quantum-chemical simulations of the potential energy surface.
This thesis also contains contributions to instrumentation development for investigations of photo-induced molecular dynamics at high brilliance X-ray light sources. We characterize the setup used for the transient valence-excited state XAS measurements of 2-thiopyridone. Therein, a sub-micrometer thin liquid sample environment is established employing in-vacuum flat-jet technology, which enables a transmission experimental geometry. In combination with a MHz-laser system, we achieve a high detection sensitivity for photo-induced X-ray absorption changes. Additionally, we present conceptual improvements for temporal X-ray optical cross-correlation techniques based on transient changes of multilayer optical properties, which are crucial for the realization of femtosecond time-resolved studies at synchrotrons and free-electron lasers.
An einigen CT-Modellkomplexen in verschiedenen Lösungsmitteln und bei Temperaturen von 113-300 K sollte der Einfluß der Umgebung auf die Form und Lage der Absorption von CT-Komplexen unterschiedlicher Bindungsfestigkeit untersucht werden.
Dazu wurden bekannte Bandenprofilfunktionen auf ihre Anwendbar-keit geprüft. Da eine optimale Anpassung nicht möglich war, wurde eine neue Profilfunktion entwickelt, die eine bessere Beschreibung ergab.
Nach der Bestimmung der Gleichgewichtskonstante und des Extink-tionskoeffizienten konnte mit der Profilfläche das Übergangsmoment berechnet werden.
Die Lösungsmittelabhängigkeit wurde bei verschiedenen Brechzahlen und Dielektrizitätskonstanten untersucht.
Für feste Komplexe wurde eine spezielle Präparationstechnik gewählt. Die beobachteten Feinstrukturen und der auftretende Streuuntergrund werden diskutiert.
Absorption Tails of Donor
(2017)
In disordered organic semiconductors, the transfer of a rather localized charge carrier from one site to another triggers a deformation of the molecular structure quantified by the intramolecular relaxation energy. A similar structural relaxation occurs upon population of intermolecular charge-transfer (CT) states formed at organic electron donor (D)-acceptor (A) interfaces. Weak CT absorption bands for D A complexes occur at photon energies below the optical gaps of both the donors and the C-60 acceptor as a result of optical transitions from the neutral ground state to the ionic CT state. In this work, we show that temperature-activated intramolecular vibrations of the ground state play a major role in determining the line shape of such CT absorption bands. This allows us to extract values for the relaxation energy related to the geometry change from neutral to ionic CT complexes. Experimental values for the relaxation energies of 20 D:C-60 CT complexes correlate with values calculated within density functional theory. These results provide an experimental method for determining the polaron relaxation energy in solid-state organic D-A blends and show the importance of a reduced relaxation energy, which we introduce to characterize thermally activated CT processes.
We study the diffusive motion of a particle in a subharmonic potential of the form U(x) = |x|( c ) (0 < c < 2) driven by long-range correlated, stationary fractional Gaussian noise xi ( alpha )(t) with 0 < alpha <= 2. In the absence of the potential the particle exhibits free fractional Brownian motion with anomalous diffusion exponent alpha. While for an harmonic external potential the dynamics converges to a Gaussian stationary state, from extensive numerical analysis we here demonstrate that stationary states for shallower than harmonic potentials exist only as long as the relation c > 2(1 - 1/alpha) holds. We analyse the motion in terms of the mean squared displacement and (when it exists) the stationary probability density function. Moreover we discuss analogies of non-stationarity of Levy flights in shallow external potentials.
We have numerically studied the bifurcation properties of a sheet pinch with impenetrable stress-free boundaries. An incompressible, electrically conducting fluid with spatially and temporally uniform kinematic viscosity and magnetic diffusivity is confined between planes at x1=0 and 1. Periodic boundary conditions are assumed in the x2 and x3 directions and the magnetofluid is driven by an electric field in the x3 direction, prescribed on the boundary planes. There is a stationary basic state with the fluid at rest and a uniform current J=(0,0,J3). Surprisingly, this basic state proves to be stable and apparently to be the only time-asymptotic state, no matter how strong the applied electric field and irrespective of the other control parameters of the system, namely, the magnetic Prandtl number, the spatial periods L2 and L3 in the x2 and x3 directions, and the mean values B¯2 and B¯3 of the magnetic-field components in these directions.
For hybrid organic-inorganic as well as all-inorganic lead halide perovskites a Rashba effect has been invoked to explain the high efficiency in energy conversion by prohibiting direct recombination. Both a bulk and surface Rashba effect have been predicted. In the valence band of methylammonium (MA) lead bromide a Rashba effect has been reported by angle-resolved photoemission and circular dichroism with giant values of 7-11 eV angstrom. We present band dispersion measurements of MAPbBr(3) and spin-resolved photoemission of CsPbBr3 to show that a large Rashba effect detectable by photoemission or circular dichroism does not exist and cannot be the origin of the high effciency.
In Allefeld & Kurths [2004], we introduced an approach to multivariate phase synchronization analysis in the form of a Synchronization Cluster Analysis (SCA). A statistical model of a synchronization cluster was described, and an abbreviated instruction on how to apply this model to empirical data was given, while an implementation of the corresponding algorithm was (and is) available from the authors. In this letter, the complete details on how the data analysis algorithm is to be derived from the model are filled in.
The nitrogen and sulfur K-edge X-ray absorption spectra of aqueous 2-thiopyridone, a model system for excited-state proton transfer in several recent time-resolved measurements, have been simulated from ab initio molecular dynamics. Spectral signatures of the local intra- and inter-molecular structure are identified and rationalized, which facilitates experimental interpretation and optimization. In particular, comparison of aqueous and gas phase spectrum simulations assesses the previously unquantified solvatization effects, where hydrogen bonding is found to yield solvatochromatic shifts up to nearly 1 eV of the main peak positions. Thereby, while each K-edge can still decisively determine the local protonation of its core-excited site, only their combined, complementary fingerprints allow separating all of the three relevant molecular forms, giving a complete picture of the proton transfer.
Gamma-ray bursts (GRBs) are brief flashes of gamma-rays and are considered to be the most energetic explosive phenomena in the Universe(1). The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow-produced by the interaction between the ejected matter and the circumburst medium-slows down, and a gradual decrease in brightness is observed(2). GRBs typically emit most of their energy via.-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments(3). However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elusive(4). Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow-ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies.
We present results from deep observations toward the Cygnus region using 300 hr of very high energy (VHE)gamma-ray data taken with the VERITAS Cerenkov telescope array and over 7 yr of high-energy.-ray data taken with the Fermi satellite at an energy above 1 GeV. As the brightest region of diffuse gamma-ray emission in the northern sky, the Cygnus region provides a promising area to probe the origins of cosmic rays. We report the identification of a potential Fermi-LAT counterpart to VER J2031+415 (TeV J2032+4130) and resolve the extended VHE source VER J2019+368 into two source candidates (VER J2018+367* and VER J2020+368*) and characterize their energy spectra. The Fermi-LAT morphology of 3FGL J2021.0+4031e (the Gamma Cygni supernova remnant) was examined, and a region of enhanced emission coincident with VER J2019+407 was identified and jointly fit with the VERITAS data. By modeling 3FGL J2015.6+3709 as two sources, one located at the location of the pulsar wind nebula CTB 87 and one at the quasar QSO J2015+371, a continuous spectrum from 1 GeV to 10 TeV was extracted for VER J2016+371 (CTB 87). An additional 71 locations coincident with Fermi-LAT sources and other potential objects of interest were tested for VHE gamma-ray emission, with no emission detected and upper limits on the differential flux placed at an average of 2.3% of the Crab Nebula flux. We interpret these observations in a multiwavelength context and present the most detailed gamma-ray view of the region to date.
We demonstrate a tilted pulse-front transient grating (TG) technique that allows to optimally utilize time resolution as well as TG line density while probing under grazing incidence as typically done in extreme ultraviolet (EUV) or soft x-ray (SXR) experiments. Our optical setup adapts the pulse front tilt of the two pulses that create the TG to the grazing incident pulse. We demonstrate the technique using all 800 nm femtosecond laser pulses for TG generation on a vanadium dioxide film. We probe that grating via diffraction of a third 800 nm pulse. The time resolution of 90 fs is an improvement by a factor of 30 compared to our previous experiments on the same system. The scheme paves the way for EUV and SXR probing of optically induced TGs on any material.
The high resolution near edge X-ray absorption fine structure spectrum of nitrogen displays the vibrational structure of the core-excited states. This makes nitrogen well suited for assessing the accuracy of different electronic structure methods for core excitations. We report high resolution experimental measurements performed at the SOLEIL synchrotron facility. These are compared with theoretical spectra calculated using coupled cluster theory and algebraic diagrammatic construction theory. The coupled cluster singles and doubles with perturbative triples model known as CC3 is shown to accurately reproduce the experimental excitation energies as well as the spacing of the vibrational transitions. The computational results are also shown to be systematically improved within the coupled cluster hierarchy, with the coupled cluster singles, doubles, triples, and quadruples method faithfully reproducing the experimental vibrational structure. Published by AIP Publishing.
A task-based parallel elliptic solver for numerical relativity with discontinuous Galerkin methods
(2022)
Elliptic partial differential equations are ubiquitous in physics. In numerical relativity---the study of computational solutions to the Einstein field equations of general relativity---elliptic equations govern the initial data that seed every simulation of merging black holes and neutron stars. In the quest to produce detailed numerical simulations of these most cataclysmic astrophysical events in our Universe, numerical relativists resort to the vast computing power offered by current and future supercomputers. To leverage these computational resources, numerical codes for the time evolution of general-relativistic initial value problems are being developed with a renewed focus on parallelization and computational efficiency. Their capability to solve elliptic problems for accurate initial data must keep pace with the increasing detail of the simulations, but elliptic problems are traditionally hard to parallelize effectively.
In this thesis, I develop new numerical methods to solve elliptic partial differential equations on computing clusters, with a focus on initial data for orbiting black holes and neutron stars. I develop a discontinuous Galerkin scheme for a wide range of elliptic equations, and a stack of task-based parallel algorithms for their iterative solution. The resulting multigrid-Schwarz preconditioned Newton-Krylov elliptic solver proves capable of parallelizing over 200 million degrees of freedom to at least a few thousand cores, and already solves initial data for a black hole binary about ten times faster than the numerical relativity code SpEC. I also demonstrate the applicability of the new elliptic solver across physical disciplines, simulating the thermal noise in thin mirror coatings of interferometric gravitational-wave detectors to unprecedented accuracy. The elliptic solver is implemented in the new open-source SpECTRE numerical relativity code, and set up to support simulations of astrophysical scenarios for the emerging era of gravitational-wave and multimessenger astronomy.
On 2015 March 23, the Very Energetic Radiation Imaging Telescope Array System (VERITAS) responded to a Swift-Burst Alert Telescope (BAT) detection of a gamma-ray burst, with observations beginning 270 s after the onset of BAT emission, and only 135 s after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40-minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant because the very-high-energy (VHE) observation started only similar to 2 minutes after the prompt emission peaked, and Fermi-Large Area Telescope observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB 150323A (z = 0.593) limits the attenuation by the extragalactic background light to similar to 50% at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below similar to 100 GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be A greater than or similar to 3 x 10(11) g . cm(-1), consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the interstellar medium, which therefore cannot be ruled out as the environment of GRB 150323A.
Context. Galactic globular clusters (GCs) are now known to host multiple populations displaying particular abundance variations. The different populations within a GC can be well distinguished following their position in the pseudo two-colors diagrams, also referred to as "chromosome maps". These maps are constructed using optical and near-UV photometry available from the Hubble Space Telescope (HST) UV survey of GCs. However, the chemical tagging of the various populations in the chromosome maps is hampered by the fact that HST photometry and elemental abundances are both only available for a limited number of stars. Aims. The spectra collected as part of the MUSE survey of globular clusters provide a spectroscopic counterpart to the HST photometric catalogs covering the central regions of GCs. In this paper, we use the MUSE spectra of 1115 red giant branch (RGB) stars in NGC 2808 to characterize the abundance variations seen in the multiple populations of this cluster. Methods. We used the chromosome map of NGC 2808 to divide the RGB stars into their respective populations. We then combined the spectra of all stars belonging to a given population, resulting in one high signal-to-noise ratio spectrum representative of each population. Results. Variations in the spectral lines of O, Na, Mg, and Al are clearly detected among four of the populations. In order to quantify these variations, we measured equivalent width differences and created synthetic populations spectra that were used to determine abundance variations with respect to the primordial population of the cluster. Our results are in good agreement with the values expected from previous studies based on high-resolution spectroscopy. We do not see any significant variations in the spectral lines of Ca, K, and Ba. We also do not detect abundance variations among the stars belonging to the primordial population of NGC 2808. Conclusions. We demonstrate that in spite of their low resolution, the MUSE spectra can be used to investigate abundance variations in the context of multiple populations.