abstract: The experimental conditions by which electromagnetic signals (EMS) of low frequency can be emitted by diluted aqueous solutions of some bacterial and viral DNAs are described. That the recorded EMS and nanostructures induced in water carry the DNA information (sequence) is shown by retrieval of that same DNA by classical PCR amplification using the TAQ polymerase, including both primers and nucleotides. Moreover, such a transduction process has also been observed in living human cells exposed to EMS irradiation. These experiments suggest that coherent long range molecular interaction must be at work in water so to allow the observed features. The quantum field theory analysis of the phenomenon is presented.
author: Luc Montagnier${}^{a,b}$; Emilio Del Giudice${}^{c,}$[^1]; Jamal Aïssa${}^{b}$; Claude Lavallee${}^{a}$; Steven Motschwiller${}^{d}$; Antonio Capolupo${}^{e,f}$; Albino Polcari${}^{g}$; Paola Romano${}^{g,h}$; Alberto Tedeschi${}^{i}$; Giuseppe Vitiello${}^{e,f}$
title: Transduction of DNA information through water and electromagnetic waves

# INTRODUCTION

This paper is an overview of what we have achieved during the past ten years in this new field of Biology: the role of water memory and electromagnetic waves in biological processes, including pathological conditions. The reported data is not only of theoretical interest, but leads to many medical applications.

This work could not have been done and analyzed without the constant interaction of biologists and physicists. The quantum field theoretical analysis of the phenomenon points to the crucial role played by coherent molecular dynamics.

# ELECTROMAGNETIC SIGNALING OF DNA

##  The detection of electromagnetic signals (EMS)

On 13 July 2005, (the eve of Bastille Day in France), by using a device previously designed by the Jacques Benveniste team to detect electromagnetic signals in water dilutions of biologically active compounds, and with the help of one of his former collaborators, Dr. Jamal Aïssa, two of us (LM, JA) observed for the first time an increase of amplitude and frequency of the recorded electric signals emitted by some high dilutions of filtrates of bacteria (Mycoplasma pirum, then Escherichia coli). This was the beginning of an extensive investigation on the role and the molecular origin of this new phenomenon (Montagnier, Aïssa, Ferris, et al., 2009; Montagnier, Aïssa, Lavallee, et al., 2009; Montagnier, Aïssa, Del Giudice, et al., 2011).

We soon discovered that DNA was the main source of the initiation of electromagnetic signals in water. In contrast to the fresh preparation of biological fluids (blood plasma, culture media) which lose their capacity of inducing EMS in water upon freezing, DNA extraction could be done from frozen material without losing its EMS capacity.

In fact, as we will see below, this property of some bacterial and viral DNA sequences of emitting EMS is like an indelible tag, and is faithfully transmitted to water structures. The bacterial species with pathological potential cultured in standard growth media yield DNA with EMS capacity. However, we noticed that one apathogenic strain of E. coli used for DNA cloning lacks this capacity as does a probiotic bacterium (Lactobacillus).

The size of DNA fragments emitting EMS ranges between 104 base pair (LTR fragment of HIV) to several kilo-bases (adhesin gene of Mycoplasma pirum: 1.5 - 3 kbp). Some PCR (Polymerase Chain Reaction) amplicons[^2], of 400 - 500 bp, have been found good emitters for the transduction experiments (see below) . For the capture of EMS in water dilutions, the conditions are very strict, and are the same for the extracted DNA as for the fresh unfrozen samples of plasma or of culture medium:

: through 450 nm and then 100 nm Millipore filters, for detection of EMS having a bacterial origin; or 450 nm then 20 nm (Whatman anotop) for EMS of viral origin (only tested for some small virus DNA and HIV DNA). The usual starting concentration is 2 ng of DNA/1 ml diluted 100 times (10 mls) for filtration.

: several decimal dilutions are made in conic plastic tubes (Eppendorf) usually 0.1 ml/0.9 ml of water under a laminar flow hood. (Fig. 1) Strong vortex shaking (for 15 seconds) is made at each dilution at room temperature. Water is purchased from commercial firms (usually 5 prime water, DNase - RNase free from Sigma). Usually, the EMS-emitting dilutions are between $10^{-7}$ to $10^{-13}$ for bacteria, $10^{-6}$ to $10^{-10}$ for small viruses (as HIV1). The first lower dilutions are apparently "silent", not emitting detectable signals.

: As the EMS are produced by resonance,(Montagnier, Aïssa, Del Giudice, et al., 2011) we use either an artificial excitation (7 Hz which was found to be the minimally active frequency) in a closed environment shielded by mu-metal, or in open air exposed to the background ambient noise. Usually, this background noise predominates in the 50 - 300 Hz range, so that the positive EMS, which are in the range of 500 - 3000 Hz, are easily detected.

In the measurement room, cell phones should be turned off (battery removed) as some phones are regulated by low frequency signals.

## Evidence that EMS emission depends on specific modification of the DNA molecule

The diversity of the DNA sequences emitting EMS does not indicate any clue as to attribute this EMS emission property to specific sequences.

However, we have studied an interesting situation in the case of HIV-infected patients: here, besides EMS produced by HIV DNA (nanostuctures filtering at 20 nm pores), we detected EMS filtering at 100nm pores and retained at 20 nm pores produced by DNA of an intracellular bacterium present in red blood cells. Surprisingly, these "bacterial" EMS were found to be produced in part by human DNA sequences integrated in or strongly associated with the bacterial DNA. The same DNA sequences belonging to the chromosomal genome of the same patient never produced EMS.

Moreover, the same sequence was found present in the red blood cells of some healthy individuals, HIV negative; but in these HIV negative individuals this sequence was found not to emit signals. This would indicate that the modification of this DNA resulting in EMS emission occurred only under pathogenic conditions. This modification was maintained in all molecules derived by PCR amplification (amplicon).

As EMS are so far only detected in patients suffering of various chronic diseases, it is tempting to speculate that there is a common biochemical modification of the DNA of infectious bacteria and/or viruses present in such diseases. This modification, which remains to be determined, should be different from a base change (mutation), since there is no difference in the base sequence of the previously mentioned amplicons in diseased patients, compared to healthy individuals.

##  Water nanostructures and EMS do carry DNA information

Our formerly reported experiments (Fig. 2 and ref. (Montagnier, Aïssa, Lavallee, et al., 2009)) indicate that the ability of EMS production can be transmitted from tube 1 containing an emitter DNA dilution to tube 2 of "naive" water, provided the system is excited overnight by electromagnetic waves of a minimal frequency of 7 Hz. Presumably tube 1 transmits waves to the water in tube 2, which did not originally contain any trace of the DNA at the origin of the signals.

The emission of EMS by the exposed tubes is thus a resonance phenomenon, dependent on external wave input. More importantly, these EMS carry specific information of the initial DNA, as shown by retrieving the DNA by PCR in the recipient tube.

This experiment has been repeated many times in our laboratory, with extraordinary precautions taken to avoid contamination in the PCR step, and many controls were always done. Omission of any of the main parameters of the procedure (7 Hz excitation, mu metal shield, time of exposure to the 7 Hz excitation, any ingredient of the PCR) as well as any minor detail of the protocol will result in failure of the experiment.

To further make the demonstration unassailable, the EMS carrying the DNA information were recorded as a digital file and sent via Internet to a recipient laboratory where work on this DNA or on the bacterium or virus which was the source of that DNA had never been done (Fig. 3). Several labs in Italy and Germany accepted the challenge.

Here, as an example, we show the results obtained at Göttingen University using a recorded file (digitized in a lap top computer in our laboratory) of ribosomal 16S DNA from Borrelia burgdorferi (Fig. 4).

In the German laboratory, the electric current resulting from the digital file communicated by our lab was converted to analog and was amplified. The current was then connected to a solenoid. A tube of water was inserted in the solenoid and in this way was submitted to the induced modulated magnetic field for one hour. Then the PCR ingredients were introduced in an aliquot of water from the tube, and after 40 PCR cycles of amplification the original DNA was detected, as shown by a specific band in gel electrophoresis of the expected molecular weight.

These intriguing results raised several questions :

1\) How a DNA polymerase (the TAQ polymerase of a thermophilic bacterium) can "read" a genetic code on water structures?

2\) What about other DNA polymerases of procaryotic and eucaryotic cells? Do they have the same capacity?

Although still at its early stages, the theoretical study of how water structures can store molecular information and transport it by electromagnetic waves gives a crucial role to the coherent molecular dynamics in the formation of water nanostructures (see below and (Montagnier, Aïssa, Del Giudice, et al., 2011)). We need, however, further theoretical analysis for a complete understanding of the phenomenon, especially because we have recent evidence that some other DNA polymerases have the same capacity as the TAQ polymerase to read water messages and can act in living cells.

## Transduction of DNA in living cells

The modified DNA transduction system is shown in figure 5.

Instead of magnetizing water in a tube placed inside a solenoid reading the modulated current from the recorded EMS signature, we placed inside the solenoid a flask containing cultured cells; the flask was placed in a vertical position for cells growing in suspension, and in a horizontal position for cells adhering to the surface of the flask. The voltage (between 2 - 4 volts) applied to the solenoid was adjusted in order to not generate heat damaging the cultured cells. This weak intensity was compensated by the duration time of exposure, between 5 to 10 days. A control flask was always placed outside the solenoid in the same 37 ${}^{\circ} C$ incubator as the exposed flask.

We used several recorded EMS files, including the 16S Borrelia and the 194 bp HIV1 LTR amplicon all having been previously shown to be good at transducing their DNA through water.

We tested several immortalized human cell lines derived from leukemias, or cancers: the HL60, originated from a myeloblastic leukemia, the U937, derived from a lung lymphoma, the MCF7, derived from a breast adenocarcinoma. In addition, we tested normal cells: the MRC5 diploid fibroblast cell line, derived from the lung of a human embryo, T lymphocytes from a healthy blood donor activated with PHA, and interleukin 2.

Results were striking: All cells of tumor origin synthesized Borrelia 16S DNA after they were exposed for several days to magnetic field modulated by the EMS of Borellia 16S DNA. At the same time, cell growth was inhibited, ending in cell death. DNA was extracted from the dying cells and the Borrelia amplicon was detected by PCR.

Remarkably the resulting amplicon was found to be EMS emitter, showing that this initial property was not lost during the complex transmission of DNA information. The normal differenciated MRC5 cells and the T lymphocytes were not affected in their growth under the same culture conditions and the Borrelia amplicon could not be detected in these cells. The 194 bp HIV LTR amplicon had no effect on the tumor cells.

These preliminary results indicate that the tumoral cell lines so far investigated do possess the enzymatic ability of reading the water nanostructures carrying the DNA information. It remains to be determined whether or not normal embryonic totipotent stem cells have the same ability to read the DNA sequence signals.

# THEORETICAL ANALYSIS

In the previous Section we have reported the experimental observation that EMS can be emitted by diluted aqueous solutions of bacterial and viral DNA under proper conditions. Moreover, it has been observed that duplication of the emitting DNA segment can be obtained by using pure water exposed to the corresponding DNA EMS and, upon addition of enzymes, primers, etc., submitted to PCR cycles. Such a transduction process has been observed to occur also in EMS exposed living cells of tumoral origin. These experimental observations suggest that long range molecular interaction must be at work in water so to allow the observed properties. Indeed, since in the transduction process the high level of sequential ordering among several hundreds of nucleotides entering the transduced DNA chain is obtained, we are clearly in the presence of collective molecular dynamical behavior of water. In quantum field theory (QFT) it is known that the ordering of the elementary components of a system is achieved as a result of the spontaneous breakdown of symmetry and constitutes the observable manifestation of coherence (Blasone, Jizba and Vitiello, 2011; Fröhlich, 1977; Umezawa, 1993; Vitiello, 1998;). Ordering is thus not the result of short range forces, but of long range collective coherent correlation. The classical behavior of the ordered pattern derives from the fact that in coherent states the ratio between the quantum fluctuation $\Delta n$ in the correlation modes and their condensate number $n$ is $\Delta n/n = 1/{|\alpha|}$ and quantum fluctuations are thus negligible for high $|\alpha|$, which denotes the coherent strength. In the present case, the symmetry which gets broken is the rotational symmetry of the electrical dipoles of the water molecules and correlation modes are the ones associated to the dipole waves (similar to spin waves in ferromagnets) (Del Giudice, Doglia, Milani and Vitiello, 1985, 1986).

We thus conclude that the observed properties of the DNA-water system provide an indication of (may be accounted by) the coherent molecular dynamics. The theoretical analysis based on quantum electrodynamics (QED) shows (Montagnier, Aïssa, Del Giudice, et al., 2011) that liquid water appears to behave as an active medium able to perform through very low frequency electromagnetic fields (e.m.f.). Short range H-bond and electric dipole-dipole static interactions among liquid water molecules set in as the consequence of the molecule interaction with time-dependent radiative e.m.f. over an extended region called coherence domain (CD) (Del Giudice, Doglia, Milani and Vitiello, 1985, 1986; Del Giudice, Preparata and Vitiello, 1988; Del Giudice and Vitiello, 2006; Bono, Del Giudice, Gamberale and Henry, 2012). Short range H-bond and electric dipole-dipole static interactions are themselves the dynamical effects caused by the most fundamental long range molecular and radiative e.m.f. interaction. This last one is thus responsible for the dynamic origin of short range interactions. This can be better understood by recalling a few points of the discussion presented in (Montagnier, Aïssa, Del Giudice, et al., 2011).

Above a density threshold and below a critical temperature, an ensemble of molecules interacting with the e.m.f. undergoes a transition to a dynamical regime characterized by a minimum energy state where the phase oscillations of the molecules are no longer uncorrelated. Such a minimum energy state implies a configuration of the system where all molecules enclosed within the CD oscillate in unison in tune with the e.m.f. trapped within the CD (phase locking). The linear size of the CD is determined by the wavelength $\lambda$ of the trapped e.m.f. (typically of the order of 100 nm). The dynamical mechanism ruling the CD formation is the one of the spontaneous breakdown of symmetry and it is described in (Del Giudice, Doglia, Milani and Vitiello, 1985, 1986; Del Giudice, Preparata and Vitiello, 1988; Del Giudice, Spinetti and Tedeschi, 2010; Del Giudice and Tedeschi, 2009; Del Giudice and Vitiello, 2006). Its mathematical formulation (Matsumoto, Papastematiou, Umezawa and Vitiello, 1975) is similar to the one of the Anderson-Higgs-Kibble mechanism (Anderson, 1958; Higgs, 1966; Kibble, 1967) which has led to the recent discovery of the Higgs particle. One important aspect of such a general QFT mechanism is that the transition to the coherent dynamical regime can be triggered by a vanishingly weak external input. Due to the weakness of the input, the system does not get "slaved" by it, but reacts to it according to its own internal dynamics and, provided that the mentioned conditions of temperature and density are satisfied, the system sets in a coherent state, whose phase is determined by the phase of the triggering input (Blasone, Jizba and Vitiello, 2011; Umezawa, 1993). Its coherence strength, however, does not depend on the input strength[^3]. In the Appendix we will comment on the question whether coherence is not destroyed by, or theoretically incompatible with the decoherence phenomenon in quantum mechanics (QM). This is an important issue, which shows that the DNA transduction process is indeed a *quantum field theory* process; it could not be understood and described in QM, where the decoherence phenomenon does in fact occur (Alfinito, Viglione and Vitiello, 2001). Our framework, however, is the one of QFT.

These features already help us in the understanding of some of the experimental observations. In particular, it is immediately recognized the observed relevance of extremely low signal (ELS) in the phenomena under study. The stimulation caused by the electromagnetic background of very low frequency is indeed observed to be essential in order for the DNA-water system to emit the EMS. In the experiments, the background ELS is either produced from natural sources (the Schumann resonances which start at 7.83 Hz (Montagnier, Aïssa, Del Giudice, et al., 2011; Nickolaenko and Hayakawa, 2002)) or from artificial sources.

The fragment of DNA embedded in the water also acts as a trigger in the surrounding water, causing the spontaneous formation of CDs, which appear as a self-produced cavity with molecular coherence strength independent from the specific input strength. However, there is "phase locking" between the specific DNA molecular structure and the water molecules. Such a specific feature of the DNA-water coherent coupling accounts for the experimental observations. In the part of the experiment concerning the DNA transduction, the dynamical phase locking is shared with pure water in a tube when it is irradiated by the EMS emitted by the aqueous DNA solution system. For brevity, we omit to report further analysis of the interplay between the size of the CD, the e.m.f. wavelength and the e.m.f. self-trapping and frequency inside the CD. For the reader convenience, we report in the Appendix a brief summary of some other features of CD discussed in (Montagnier, Aïssa, Del Giudice, et al., 2011).

The question now naturally arises: does the EMS have any specific property related to the coherent dynamical structure discussed above? The question is particularly relevant because the emitted EMS, acting on water molecular dynamics, produces coherent structures such that in PCR processes the DNA transduction occurs with the same nucleotide sequence as the one of the parent DNA. The answer to the question is provided by observing that the EMS appears to carry not only the specific information of its frequency spectrum, amplitude and phase modulation (the syntactic level), but it also describes the dynamics out of which it is generated. In other words, beside the syntactic level of pure information (à la Shannon), there is a *semantic* content, which manifests itself in the underlying coherent dynamics of the DNA-water system responsible of the polymerization (highly ordered sequence) of hundreds of nucleotides. We refer to such a semantic content as to the "meaning" of the EMS. In some recent papers it has been shown that an isomorphism exists between (squeezed) coherent states and self-similar fractal properties (Vitiello, 2009a, 2009b, 2012, 2014). Here, for brevity we do not proceed further with our analysis on this point. Results will be presented elsewhere.

Above we have mentioned the mechanism of phase locking and phase content in water molecular dynamics. Let us close this Section by stressing indeed the crucial role played by the phase in the considered processes. Due to the relation between phase and electromagnetic potential, non-trivial topological properties, with associated Bohm-Aharonov-like effects, have a non-secondary part in the molecular dynamical properties of water (see e.g. (Del Giudice, Fuchs and Vitiello, 2010; Del Giudice and Vitiello, 2006) and discussions there reported). Such a remark may turn out to be important when considering charge transfer along the double helix, produced by oxidative agents observed in (Genereux and Barton, 2010). The question if such a current could contribute to the production of EMS may have an affirmative answer since the charge transfer along the DNA produces magnetic field in the surroundings. On the other hand, the oscillations of electric dipoles of the DNA macromolecule may propagate on the DNA in wave form so to contribute to the electromagnetic signal emission. These EMS are considered in (Del Giudice, Doglia, Milani and Vitiello, 1985, 1986) where it is shown that they produce symmetry breakdown in the water in which the dipole chain (DNA or protein chain) is embedded and the mathematical details of the proof are reported.

Finally, concerning the question 1 in Section II.C (how a DNA polymerase, the TAQ polymerase of a thermophilic bacterium, can "read" a genetic code on water structures) the whole dynamical scenario above presented provides the answer to it. It is a complex scenario founded on QFT of coherent systems and is then not surprising that those who are unaware of it could not conceive the positive answer, with all of its complex but clear details, to that question.

# CONCLUSION AND PERSPECTIVES

This 10 year long collaborative work has yielded some scientific facts and concepts in a new domain of Science at the frontier of Biology and quantum field Physics.

A new property of some DNA molecules has been discovered, that of emitting low frequency electromagnetic waves in water dilutions. These DNAs originate from pathogenic agents or agents endowed with pathogenic potential. It may be not pure coincidence that such EMS are associated with diseases, particularly chronic diseases.

Under natural conditions EMS and water nanostructures may play a role of stealthy elements carrying DNA information of infectious agents while being undetected by the immune system or being insensitive to conventional therapies (Fig. 6). However, one cannot discard the possibility that DNA waves can play a role in the physiology of living entities.

Moreover, in the laboratory, we have shown for the first time that EMS can be re-transcribed into DNA in living cells. These cells are so far of tumoral origin, opening the way to non-invasive treatments of cancers, assuming that normal stem cells are not affected, or less affected.

Thus, this new biology that we can call after Jacques Benveniste, Digital Biology, has a very promising future, both at the level of quantum Physics, and in numerous medical applications.

From the point of view of the theoretical understanding of the observed phenomena, the discussion above presented suggests that the dynamical law of coherence acts as a *law of form* inducing guided polymerization processes (controlling morphogenesis): the specific polymerization so obtained is the expression of the semantic content of the EMS mentioned in the previous Section, by us referred to as the signal meaning. The dynamics of coherence appears to play the role of dynamic paradigm ruling the ordering (polymerization) processes through dissipative non-equilibrium dynamics controlled by entropy variations and the consequent appearance of the arrow of time (breakdown of time-reversal symmetry) (Celeghini, Rasetti and Vitiello, 1992; Vitiello, 2012, 2014).

The experiments discussed in this paper suggest that also in the *usual* PCR processes the DNA duplication is obtained due to the EMS emitted by the parent DNA in the environment of reciprocal inter-actions with water molecules, enzymes, primers and nucleotides in the solution.

The EMS appears thus to be the carrier of the coherence expressed in the DNA code. One might conjecture (Vitiello, 2014) that modifications are induced in the properties of the EMS resulting in the "deformation"[^4] of coherence (e.g. such as those, but not only those, induced by the observed bacterial actions; cf. Section II.B). This may play a role in epigenetic modifications, thus revealing the appearance of "new meanings" (in the above mentioned sense) associated to deformed properties of EMS. DNA appears to be the *vehicle* through which coherence and its dynamical deformations propagate in living matter (Vitiello, 2014).

# Acknowledgements

Partial financial support from MIUR and INFN is acknowledged. We thank Pr. E. Schutz, and Dr. H. Urnovitz, Chronix Biomedical, for allowing transduction experiments in their laboratory. Mrs. Laila Aïssa is acknowledged for her skilled technician assistance. Mrs. S. McDonnell is acknowledged for her constant participation in the management of this work.

# On some properties of water molecular dynamics

One of the peculiarities of water consists in the fact that water molecules in the CD coherently oscillate between the ground state and an excited state lying at 12.06 eV, just below the 12.60 eV ionization threshold (Bono, Del Giudice, Gamberale and Henry, 2012; Montagnier, Aïssa, Del Giudice, et al., 2011). The almost free electrons in the CD can be excited by external inputs so to form coherent excitations (vortices), whose entropy is lower than the entropy of the energy supplied by the inputs. Vortices, due to coherence and their non-trivial topology, are not easily destroyed by small, additional external supplies of energy. On the contrary, additional energetic inputs may add up to form a unique vortex, thus storing in the CDs an amount of energy which may be large enough to activate chemical reactions among molecules, otherwise below the activation energy threshold. Thus, small energy contributions coming from many high entropy inputs add up to form low entropy ordered patterns of upgraded high energy (Del Giudice, Fuchs and Vitiello, 2010; Voeikov and Del Giudice, 2009).

An important remark is that DNA and proteins are polyelectrolytes, and are surrounded by positive counterions. Ions having a cyclotron frequency, $\nu_{c} = q \, B/(2\, \pi\, m)$, where $q$ and $m$ are the electric charge and the mass of the ion, respectively, and $B$ is the magnetic field (Liboff, 1997; Liboff, Smith and McLeod, 1987), may play an important role in obtaining a collective performance of water CDs, a coherence of coherent domains. The observed dependence of the signal emission on the aqueous dilution may be understood (Montagnier, Aïssa, Del Giudice, et al., 2011) as follows: suppose that a low magnetic field (for example a natural or artificially produced background magnetic field) matches the ion cyclotron frequency; suppose it may be then able to extract $n$ ions per CD. Then, due to angular momentum conservation, the plasma of $N$ quasi-free electrons in the CDs starts to counter-rotate with a frequency much higher than the ion cyclotron frequency since electron mass is much smaller than the ion mass. This frequency depends on the number of involved ions, namely on their concentration, which therefore is the only relevant variable. This occurs on all the CDs of the system, whose number is irrelevant for the frequency purpose, in agreement with observations. The magnetic component of EMS is produced by the so induced rotation of the plasma of the quasi-free electrons in the CDs. As a further effect, a co-resonating field appears in the surroundings of the rotating CDs depending on the ion concentration, i.e. on the DNA solution dilution. It could be at the origin of an extended coherence among CDs. The existence of the observed window of dilutions for the occurrence of the EMS emission could be understood by presuming that the signal produced by the lower dilutions could have a frequency higher than the interval of the values detectable by the used instruments. Higher dilutions, on the contrary, could produce no signal because the ion concentration is decreased below the threshold able to excite the CDs (Montagnier, Aïssa, Del Giudice, et al., 2011).

We observe that thermal collisions could be in competition with electrodynamic attraction of molecules inside the CD and produce permanent fluxes of molecules between a coherent regime and a non-coherent one, and vice-versa, although the total number of coherent and non-coherent molecules are constant for a given temperature T. Water is thus not a homogeneous liquid, rather it appears as a two fluid system, with coexisting coherent and non-coherent phases, like in the Landau theory of liquid Helium (Landau, and Lifshitz,1959). We have thus a mixed structure system, consistent also with experimental findings (Taschin, Bartolini, Eramo et. al. 2013), which may appear in observations only when observation time is very short with respect to the time of flickering between the two phases. Near surfaces, the coherent phase may be more stable due to the interaction between water molecules and the surface (Pollack, 2001, 2013). For example, in living matter, water, which is bound to membranes or to biomolecules, could more easily manifest the properties of coherence.

Let us finally consider the question whether the phenomenon of decoherence in quantum mechanics (QM) might contradict our analysis based on the formation and non-vanishing life-time of coherent structures in water. We remark that the *belief* that coherence is possible only at very low temperatures is disproved by the fact that coherence is observed in a wide range of temperature, from very low to very high ones: the diamond crystal loses its coherence (it melts) at a temperature of about $+3545 ~{}^{0}C$ in the absence of oxigen; the common kitchen salt $NaCl$ melts at $+804 ~{}^{0}C$; in the iron the coherence of the elementary magnets is lost at $+770 ~{}^{0}C$. In superconductors the critical temperature $T_c$ are much lower, in some compound of niobium they are not higher than $- 252 ~{}^{0}C$ and for some high $T_c$ superconductors it is a little above $-153 ~{}^{0}C$. Also the so called BEC systems (which are mostly condensates of atoms) need very low temperatures and are not so stable. All these systems are macroscopic systems and they are described in quantum field theory in terms of coherent Bose-Einstein condensates. The bosons that condense in a crystal are called the phonons, i.e. the quanta of the elastic waves responsible of the ordering in crystals; in the magnets, they are called the magnons, namely the quanta of the spin waves in magnets; in the water, they are called "dipole wave quanta" (DWQ), the quanta of the fluctuating molecular dipole waves; and so on. This teaches us that the ordered patterns we observe at a macroscopic scale in these systems are sustained and generated by long range correlations maintained by these waves (Alfinito, Viglione and Vitiello, 2001). One would never be able to construct any of these systems by using short range interaction among the nearest neighbours. Short range interaction, if it is there, is made possible by the long range one which brings "near" the components (e.g., making possible the formation of $H$-bonds in water). Decoherence in QM would forbid the existence of crystals, magnets, superconductors, etc.. However, these systems do exist and are observed since they are QFT systems.

# Bibliography

Alfinito, E., Viglione, R. and Vitiello, G. (2001). The decoherence criterion. Mod. Phys. Lett. B15: 127-135.

Anderson, P.W. (1958). Coherent Excited States in the Theory of Superconductivity: Gauge Invariance and the Meissner Effect. Phys. Rev. 110:827-835.

Blasone, M., Jizba, P. and Vitiello, G. (2011). Quantum field theory and its macroscopic manifestations. London: Imperial College Press.

Bono, I., Del Giudice, E., Gamberale, L. and Henry, M. (2012). Emergence of the Coherent Structure of Liquid Water. Water 4(3): 510-532. doi:10.3390/w4030510

Celeghini, E., De Martino, S., De Siena, S., Rasetti, M. and Vitiello, G. (1995). Quantum groups, coherent states, squeezing and lattice quantum mechanics. Annals of Phys.(N.Y.), 241:50-67.

Celeghini, E., Rasetti, M. and Vitiello, G. (1992). Quantum dissipation. Nucl. Phys. B215:156-170.

Del Giudice, E. Doglia, S. Milani, M. and Vitiello, G. (1985). A quantum field theoretical approach to the collective behavior of biological systems. Nucl. Phys. B251 (FS 13):375-400.

Del Giudice, E. Doglia, S. Milani, M. and Vitiello, G. (1986). Electromagnetic field and spontaneous symmetry breakdown in biological matter. Nucl. Phys. B275 (FS 17):185-199.

Del Giudice, E., Fuchs, E. C. and Vitiello, G. (2010). Collective Molecular Dynamics of a Floating Water Bridge. Water J. 2:69-82. doi.org/10.14294/WATER.2010.5

Del Giudice, E., Preparata, G. and Vitiello, G. (1988). Water as a free electron laser. Phys. Rev. Lett. 61:1085-1088.

Del Giudice, E. and Tedeschi, A. (2009). Water and the Autocatalysis in Living Matter. Electr. Biol. Med. 28(1):46-52.

Del Giudice, E., Spinetti, P. R. and Tedeschi, A. (2010). Water dynamics at the root of metamorphosis in living organisms. Water. Water J. 2:566-586.

Del Giudice, E. and Vitiello, G. (2006). Role of the electromagnetic field in the formation of domains in the process of symmetry-breaking phase transitions. Phys. Rev. A 74:022105.

Fröhlich, H. (1977). Long-range coherence in biological systems. Rivista del Nuovo Cim. 7:399-418.

Genereux, J. C. and Barton, J. K. (2010). Mechanisms for DNA Charge Transport. Chem Rev. 110(3):1642-1662.

Higgs, P. W. (1966). Spontaneous symmetry breakdown without massless bosons. Phys. Rev. 145:1156.

Kibble, T.W. B. (1967). Symmetry breaking in non-Abelian gauge theories. Phys. Rev. 155:1554

Landau, L.D. and Lifshitz, E. M. (1959). Fluids Mechanics. Reading, Mass.: Addison-Wesley, Ch. XVI.

Liboff, A. R. (1997). Electric-field ion cyclotron resonance. Bioelectromagnetics 18(1):85-7.

Liboff, A. R. Smith, S. D., McLeod, B. R. (1987). Experimental Evidence for Ion Cyclotron Resonance Mediation of Membrane Transport. In Mechanistic Approaches to Interactions of Electric and Electromagnetic Fields with Living Systems, Editors: Blank, M., Findl, E.. New York:Springer-Verlag, pp 109-132.

Matsumoto, H., Papastamatiou, L., Umezawa, H. and Vitiello, G. (1975). Dynamical rearrangement of Symmetry in the Anderson-Higgs-Kibble Mechanism. Nucl. Phys. B97:61-89.

Montagnier, L., Aïssa, J., Ferris, S., Montagnier, J-L., Lavallee, C. (2009). Electromagnetic Signals Are Produced by Aqueous Nanostructures Derived from Bacterial DNA Sequences. Interdiscip. Sci. Comput. Life Sci. 1:81-90.

Montagnier, L., Aïssa, J., Lavallee, C., Mireille Mbamy, M., Varon, J., Chenal, H. (2009). Electromagnetic detection of HIV DNA in the blood of AIDS patients treated by antiretroviral therapy. Interdiscip Sci. Comput. Life Sci. 1:245-253.

Montagnier, L., Aïssa, J., Del Giudice, E. Lavallee, C., Tedeschi, A., Vitiello, G. (2011). DNA waves and water. Journal of Physics: Conference Series 306:012007.

Nickolaenko, A. P. and Hayakawa M. (2002). Resonances in the Earth-ionosphere Cavity. Dordrecht-Boston-London: Kluwer Academic Publishers.

Pollack, G.H. (2001). Cells, Gels and the Engines of Life. Seattle, WA: Ebner and Sons.

Pollack, G.H. (2013). The fourth phase of water: Beyond Solid, Liquid, and Vapor. Seattle, WA: Ebner and Sons.

Taschin, A., Bartolini, P., Eramo, R. et. al. (2013). Evidence of two distinct local structures of water from ambient to supercooled conditions. Nature Comm. 4:2041-2045.
doi:10.1038/ncomms3041

Umezawa, H. (1993). Advanced Field theory: micro, macro and thermal concepts. New York: American Institute of Physics.

Vitiello, G. (1998). Structure and function. An open letter to Patricia Churchland. In Hameroff, S.R. Kaszniak, A. W. and Scott, A.C. Eds.. Toward a science of consciousness II. The second Tucson Discussions and debates. Cambridge: MIT Press. p. 231-236.

Vitiello, G. (2009a). Coherent states, Fractals and brain waves. New Mathematics and Natural Computation 5: 245-264.

Vitiello, G. (2009b). Fractals and the Fock-Bargmann representation of coherent states. In P. Bruza, D. Sofge, et al. Eds., Quantum Interaction. Lecture Notes in Artificial Intelligence, Edited by R.Goebel, J. Siekmann, W.Wahlster. Berlin: Springer-Verlag, p. 6-16.

Vitiello, G. (2012). Fractals, coherent states and self-similarity induced noncommutative geometry. Phys. Lett. A 376:2527-2532.

Vitiello, G. (2014). On the Isomorphism between Dissipative Systems, Fractal Self-Similarity and Electrodynamics. Toward an Integrated Vision of Nature. Systems 2:203-216. doi:10.3390/systems2020203

Voeikov, V. L. and Del Giudice, E. (2009). Water respiration: the base of the living state. Water J. 1:52-75.

Yuen, H. P. (1976). Two-photon coherent states of the radiation field. Phys. Rev. A 13:2226.

[^1]: deceased 31 January 2014

[^2]: An amplicon is a piece of DNA or RNA that is the source and/or product of natural or artificial amplification or replication events, usually PCR.

[^3]: We stress that a strong input may drive the system shielding its own internal dynamics. In such a case, the symmetry is said to be explicitly broken and one has a substantial modification of the original system by inclusion of the strong perturbing agent. However, this is not what we are interested in in the present case, and in general in Biology, where small perturbing inputs may trigger relevant reaction of the system driven by its own internal dynamics

[^4]: In the jargon of coherent states, the word deformation, or also $q$-deformation, refers to the technically well defined process of "squeezing" of the coherent state. For technical details see (Yuen, 1976; Celeghini, De Martino, De Siena, et al. 1995).


abstract: The Adaptive Gain Integrating Pixel Detector (AGIPD) is an x-ray imager, custom designed for the European x-ray Free-Electron Laser (XFEL). It is a fast, low noise integrating detector, with an adaptive gain amplifier per pixel. This has an equivalent noise of less than 1 keV when detecting single photons and, when switched into another gain state, a dynamic range of more than 10$^4$ photons of 12 keV. In burst mode the system is able to store 352 images while running at up to 6.5 MHz, which is compatible with the 4.5 MHz frame rate at the European XFEL. The AGIPD system was installed and commissioned in August 2017, and successfully used for the first experiments at the Single Particles, Clusters and Biomolecules (SPB) experimental station at the European XFEL since September 2017. This paper describes the principal components and performance parameters of the system.
author: Aschkan; Julian; Annette; Roberto; Peter; Dominic; Beat; Helmut; Manuela; Robert; Alexander; Hans; Sabine; Torsten; Alessandro; Davide; Aldo; Magdalena; Jennifer; Joern; Igor; Xintian; Sergej; Lothar; Jolanta; Ulrich; Qingqing; Mourad; Jiaguo; Manfred; Bernd
title: The Adaptive Gain Integrating Pixel Detector at the European XFEL

Allahgholi

Becker

Delfs

Dinapoli

Goettlicher

Greiffenberg

Henrich

Hirsemann

Kuhn

Klanner

Klyuev

Krueger

Lange

Laurus

Marras

Mezza

Mozzanica

Niemann

Poehlsen

Schwandt

Sheviakov

Shi

Smoljanin

Steffen

Sztuk-Dambietz

Trunk

Xia

Zeribi

Zhang

Zimmer

Schmidt

# Introduction

With the start of the European XFEL, a new milestone is set in the field of x-ray research and many related fields due to the high coherence, pulse intensity and repetition rate of the x-ray pulses available at this facility. The superconducting accelerator provides up to 600 s long trains with up to 2700 pulses followed by an inter-train gap of 99.4 ms. Inside each train, consecutive pulses of typically less than 100 fs duration are spaced approximately 220 ns apart. This corresponds to an effective repetition rate of 4.5 MHz during a train. Each pulse contains up to 10$^{12}$ photons , which in many cases is sufficient to produce a complete scattering pattern from the sample with a single pulse. This means that the area detectors at the European XFEL not only have to be compatible with the high repetition rate of the source, but also need to have a dynamic range from single photons to 10$^4$ photons/pixel/pulse. More detailed and complete requirements can be found in . Dedicated detector development programs for the European XFEL were started more than 10 years before its inauguration, as it was clear that no existing detector would be able to meet the requirements imposed by the new facility. The Adaptive Gain Integrating Pixel Detector (AGIPD) was designed to fulfill as many of the general requirements as possible with a focus on the requirements most important for scattering experiments in the energy range from 7 to 15 keV, for which is has been successfully used in first user experiments and . There has been interest to use the system at other bright sources to study fast dynamics down to the microsecond scale.

The AGIPD is not the only camera developed for use at the European XFEL. The LPD system is currently in use at the FXE beamline and the DSSC system will be available soon. Other FELs are using other custom developed camera systems like the CSPAD and ePIX detectors at LCLS or the JUNGFRAU detector at the SwissFEL.

# System layout

The AGIPD camera consists of four individually moveable quadrants, each having four detector tiles with 512 x 128 pixels per tile, giving a total of 1024 x 1024, or roughly 1 million pixels. Fig. a shows a CAD design of the AGIPD 1 million pixel detector with cuts to expose the arrangement of the electronics inside and outside of the vacuum vessel.

Each detector tile consists of a front end module (FEM), an in-vacuum board that provides power to the FEM and routes signals, two ADC boards and a control and data IO board. Fig. b and c show a CAD drawing and a photograph of these boards, respectively.

Being glued to a ceramic board, each FEM consists of a monolithic pixelated silicon sensor responsible for absorbing the x-ray photons and creating an electrical signal per pixel which proportional to the sum of the energies of all simultaneously absorbed photons in that pixel. The silicon sensor is bump-bonded to 2 x 8 pixelated Application Specific Integrated Circuits (ASICs), which are responsible for signal integration and intermediate image storage.

The front-end modules protrude into the attached sample interaction chamber and can be cooled to stabilize their temperature and improve performance. Although there are many components in the vacuum chamber, vacuum levels better than $10^{-5}$ mbar have been reached in our labs. Note that this was without attaching the system to the sample interaction chambers at the European XFEL, so the vacuum levels obtained during XFEL experiments might be different.

The detector logically and electrically divides into two wings. Each wing consists of the ADC boards and control and data IO board of each tile (one set each per tile, 8 tiles per wing) as well as a vacuum backplane board, which acts as a vacuum barrier and routes signals in and out of the vacuum vessel, a micro controller board for slow control and a master FPGA board. These boards are located outside the vacuum chamber[^1] in a thermally sealed, water cooled housing.

The two master FPGA boards, one for each side, provide the interface to the clock and control system, and control the detector tiles, including the FEMs. In the following sections the individual components and operational concepts are described in more detail.

# Mechanics

Fig. a shows the front view of the AGIPD 1 million pixel system installed at the SPB station of the European XFEL. The four independently movable quadrants allow the formation of a horizontal slit or a rectangular central hole with user selectable size for the direct beam to pass through. The movement is realized by mounting each cooling block on a motion stage formed by two wedges (Fig. b).

The cooling blocks were made by a combination of milling and electroforming techniques. In the first step, the basic shape of the cooling block was milled out of a solid block of HCP copper. This basic shape omitted details that would be included in the final milling step, but already included the cooling channels. To enhance the turbulence of the flow of the silicone oil coolant, holes were drilled into the channels and pins were inserted into the holes. Afterwards, the channels were covered by copper using an electroforming process. Initially omitted details, like connector feedthroughs, were defined in a final milling run, which also ensured the overall dimensions and tolerances of the cooling block. Fig. c and d show the channels with inserted pins before the cover deposition.

# Front end modules

An AGIPD 1 million pixel detector incorporates 16 front end modules (FEMs). Each front end module uses a bump bonded hybrid of a monolithic silicon sensor with 128 x 512 pixels and 2 x 8 AGIPD ASICs. The power and signal contacts of the ASICs are wire bonded to gold plated pads on an LTCC (Low Temperature Co-fired Ceramics) board. The bias voltage of the sensor is provided via wire bond connections between sensor and LTCC on all four corners of the assembly. A PT1000 temperature sensor is located on the backside of the LTCC, next to a Samtec 500 pin connector, which connects the FEM to the interface electronics. The hybrid assembly is glued to a silicon heat spreader, which reduces thermal gradients across the tile. For the bonding between heat spreader and LTCC an adhesive with high thermal conductance is used that is able to handle the different coefficients of thermal expansion of silicon and ceramics.

For mounting and handling, the wire bonded assembly is screwed onto an interposer made from a copper alloy which provides sufficient rigidity and heat conductivity. The interposer features an insertion pin for defined mounting to the cooling block. To reduce the thermal resistance between the layers, the interface between LTCC and interposer is filled with vacuum compatible liquid gap filler and the interface between interposer and cooling block with graphite. Fig. a shows a macro photograph of the edge of a FEM, detailing sensor, ASICs, heat spreader, LTCC and interposer. Fig. b shows two fully assembled FEMs and a bare interposer.

In our lab the system has been tested at temperatures as low as -20 C on the LTCC, but the system can run at temperatures above that, which was also the case for early user experiments at the SPB beamline. The system is also compatible with operation at ambient pressure and temperature.

## The ASIC

The AGIPD 1.1 ASIC forms the heart of the system; each ASIC incorporates 64 x 64 pixels and the necessary readout and control circuitry. It is manufactured in 130 nm CMOS technology, using radiation hardened layout techniques in most parts of the circuitry.

Fig. shows a schematic diagram of the pixel circuitry: the input of each pixel is formed by a resettable charge sensitive preamplifier, built around an inverter core. Its output feeds a discriminator and a correlated double-sampling (CDS) stage with two globally adjustable gain settings. Once the discriminator output exceeds the globally defined threshold, additional feedback capacitors of 3 or 10 pF are added to the 60 fF of the un-switched preamplifier feedback loop. This way the sensitivity of the preamplifier is adaptively decreased and the dynamic range is extended in two steps , where each pixel automatically and independently adapts its gain to the incoming signal. The CDS stage is used to remove noise from the reset switch and to suppress low frequency noise components .

The pixel response is recorded to an in-pixel memory at high speed, compatible with the 4.5 MHz requirement for the European XFEL . Each memory address contains two separate pieces of information: a) the output voltage of the CDS stage, which is proportional to the detected signal, is stored on a 200 fF capacitor and b) the gain state of the pixel, encoded as a voltage, is stored on a 30 fF capacitor.

The memory matrix occupies about 80% of the pixel area and can store up to 352 images consisting of signal and gain information. The pixel size of (200 m)$^2$ is a compromise between resolution, analog performance and number of memory cells.

The memory can be randomly accessed, providing the option of overwriting images or frame selective readout. At the European XFEL it is used to implement a veto system .

During readout of the chip another charge sensitive amplifier in each pixel is used to read the memory, which happens in parallel for each row of pixels. The further readout uses two interleaved column buses and four multiplexers, each serializing data from a block of 16 x 64 pixels. This parallelizing onto 4 outputs reduces the power consumption of the readout circuitry by reducing its speed. Instrumentation amplifiers convert the signal to differential levels, which are driven off-chip for subsequent digitization.

A command based control circuit provides all the signals for memory access, read and write operations to the pixels. It uses a 3 line serial current mode logic interface and also provides slow control tasks, like the programming of internal timings and on-chip biases generated by digital to analog converters .

## The Sensor

Much like the ASIC, the sensor design was driven by the requirements set forth by the expected experiments and performance of the European XFEL . The main requirements were a thickness of 500 m in order to reach sufficiently high quantum efficiency and a tolerance of a total dose of 1 GGy during the expected lifetime of the detector. All of this should be accomplished while making large area sensors that allow minimizing the dead area of the final detector system.

In addition to the challenge of radiation tolerance, the impact of plasma effects due to the high number of instantaneously absorbed photons was investigated . When many photons are locally absorbed, sufficient electron-hole pairs are created to form a plasma. This plasma shields its core from the electric field required to drift the charges to the readout ASICs.

Diffusion will eventually expand the plasma and lower its density enough for the drift field to take over, but this process takes time. As a result the charge cloud will spread laterally, potentially degrading the spatial resolution, and the total charge collection time will increase, potentially piling up with the next photon pulse.

These studies concluded that the sensor should be constructed using p$^+$ electrodes in a highly resistive n-type bulk, thus collecting holes. Also a voltage of at least 500 V should be applied to the AGIPD sensor to suppress the consequences of the plasma effects as much as possible. At this voltage the time to collect at least 95% of the deposited charge is less than 60 ns for tightly focused spots (3 m rms) of up to 10$^5$ 12 keV photons.

Surface damages are the dominant type of radiation damage for the sensor, as the damage threshold for silicon bulk damage is far above 12 keV, the original design photon energy of the European XFEL.

Therefore, surface damages, namely the creation of positive charges in the SiO$_2$, and the introduction of traps at the Si-SiO$_2$ interface, were studied in detail with numerous irradiation campaigns to establish the relevant parameters, which in turn were used for sensor optimization studies .

These studies showed that several design parameters, i.e. oxide thickness, pixel implant depth, and metal overhang, have significant influence on the ability to operate the sensor at a high voltage. Some of the findings show conflicting results for the situation before and after irradiation, i.e., a thick oxide is preferred before irradiation and a thin oxide after irradiation .

Finally, a compromise layout was chosen that fulfilled the design specifications and showed a breakdown voltage above 800 V in simulation . This layout was produced and its radiation tolerance was studied , finding sufficient performance for operation at the European XFEL.

The dead area of the final system is minimized for each front end module by using a monolithic sensor bump bonded to 2 x 8 individual ASICs. The large sensor guarantees that there are no blind spots within the sensitive surface (0% dead area), however since the ASICs cannot physically touch each other the pixels horizontally in between two ASICs are twice as wide (400 m x 200 m), covering the area of two ordinary pixels. Pixels vertically in between two ASICs have normal size.

The entire sensitive area is surrounded by a guard ring that is 1.2 mm (or 6 pixels) wide. A detailed description of the sensitive and non-sensitive areas and the implications for coherent diffraction imaging at the SPB instrument of the European XFEL can be found in . Disregarding intentional gap between quadrants caused by the moving apparatus each quadrant has less than 15% insensitive area, which increases to less than 18% for the whole system, still disregarding intentional gaps.

# Signal handling, data preparation and control

As indicated earlier, the AGIPD 1 million pixel detector is made of two electronically independent halves. Each half consists of a master FPGA board, a micro controller based slow control unit, an 8 port vacuum barrier board and eight tiles consisting of front end modules and their readout boards.

The 8-port vacuum barrier board, which interconnects between the different boards, is realized using a multi-layer printed circuit board which acts as a vacuum barrier and as a backplane for signal distribution.

The master FPGA acts as an interface device between the detector, the clock and control system of the European XFEL and the control PC. It receives configuration, 'start' and 'stop' signals from the control computer; receives bunch synchronized clocks and classification flags; synchronizes the operation of the ASICs and fast ADCs; and triggers the read out FPGAs of the tiles.

The boards constituting the readout electronics of each tile are located either inside the vacuum vessel or in the enclosed air environment of the external housing.

From each front end module 64 differential analog signal lines are guided via the in-vacuum board and vacuum barrier board to the analog boards. The flexible part of the in-vacuum board (indicated in Fig. a and b) allows the movement of the quadrants.

The ADCs (Analog Devices AD9257) sample at a resolution of 14 bit and operate at a frequency of up to 33 MHz, which results in a minimum possible read out time of approximately 22 ms for 352 frames (signal amplitude and gain state information are read separately, no overheads), well within the 99.4 ms inter train spacing at the European XFEL. At 14 bit resolution the ADC noise and the quantization errors do not contribute significantly to the overall noise of the system.

Each readout FPGA orders the information of the 64 ADC channels of its tile into a single frame and sends the data via a 10 GbE optical link to the data acquisition system of the European XFEL.

The slow control board monitors the status of the detector by collecting data on supply currents and voltages, as well as temperature, humidity and cooling fan information.

During the detector start-up phase the experimental control computer sends commands to the micro controller to power the electronics sequentially. Collected monitor and status information can be queried by the slow control computer of the experiment which adds time stamps and makes the data available to the XFEL system.

In addition, the slow control board serves as a second level interlock, transmitting hard wired flags to the experimental programmable logic controller (PLC) to initiate a shutdown, if operational conditions are outside pre-defined safety margins. The architecture of the electronics is discussed in more detail in .

# Calibration

The performance of the AGIPD ASIC has been extensively tested and documented over the years . Since every pixel can be in any of its three gain states in any image (i.e., memory cell) the current calibration procedure calibrates each memory cell individually for all three gain settings requiring more than 2.8 billion calibration parameters for a 1 million pixel detector system.

In a big system of one million pixels or more, calibrating all relevant parameters using only external sources is not feasible in a time efficient manner. Therefore the ASIC includes two internal calibration sources - a current source and a pulsed capacitor - to enable dynamic range scans without the need for external stimulus. Both sources combine elements that are global to the chip and local in each pixel. The current source is implemented as a distributed current mirror and each pixel features a calibration capacitor while the voltage step is generated in the periphery.

Each data point from the detector contains two values per pixel per memory cell: the analog signal value and the encoded gain state information. Each value is expressed as Analog to Digital converter Units (ADU) in the raw data stream of the detector.

The internal sources allow the determination of valuable information. For the analog signal value, that is the relative gain of all 3 states within each pixel and the relative gain of the pixels with respect to each other. For the encoded gain state used, that is the discrimination threshold that indicates the different states of the gain.

The offsets for each gain state are measured without illumination (dark frames) and the gain of the high gain state is measured for each pixel using a flood illumination with characteristic x-rays. The characteristic x-rays were generated by illuminating a metal foil with an x-ray tube. The distance between foil and sensor was approximately 10-20 cm, depending on the foil, to ensure sufficient count rate in all pixels without creating a strong gradient in count rate over the module.

In a last step the calibration data are matched and merged, resulting in independent calibration constants for each memory cell in each pixel. Each memory cell is characterized by 8 parameters: 3 offsets (in ADU, one per state), 3 gain values (in ADU/keV, one per state), and 2 thresholds (in ADU) for state discrimination. This totals more than 2.8 billion calibration parameters in a 1 million pixel detector system.

Figure shows examples of intermediate results during the calibration procedure. Fig. a shows the histogram of 10,000 analog signal values measured for a single memory cell of a single pixel. The data were acquired using illumination with characteristic x-rays from a molybdenum foil. The absolute gain of this memory cell of this pixel is given by the distance between the peaks in the histogram (noise to single photon or photon to photon). The average gain value of a typical FEM is 7.7 ADU/keV +- 3.4 %.

The integration time during data taking was increased to 50 s, a value significantly beyond the 130 ns typically used for experiments at the European XFEL. This was done to increase the number of detected photons per frame to at least 1000 and thereby reduce the statistical uncertainty of the fit. Due to the long integration time the noise of the system (width of the peaks) is much higher than for typical experimental conditions. The noise measurement of the high gain state is extracted from dark frames taken under experimental conditions (not shown) and typical values for the Equivalent Noise Charge (ENC) are 320 and 240 electrons in CDS gain low or high, respectively. For the other gain states the noise is higher, but still below the Poisson limit .

The internal current sources allow extrapolation of the absolute calibration of the high gain state to other states. The procedure using the internal current source utilizes injecting a constant current to the input of the pre-amplifier while sweeping the integration time.

An example of a current source scan is shown in Fig. b. From linear fits to the data the ratios of low, medium and high gain are determined[^2]. The high- and medium gain regime can also be scanned with a pulsed capacitor as shown in Fig. c.

The pulsed capacitor scans the dynamic range by gradually increasing the height of the applied voltage step while keeping the integration time constant. The encoded gain state values and the corresponding discrimination thresholds can be extracted from both current source and pulsed capacitor scans, and are shown by the green dots in Fig. b and c.

Recording the required data to extract the calibration data typically takes 12-14 hours and occupies more than 20 TB of disk space. The largest contribution to this is the current source scan, which typically takes 10-11 hours and generates roughly half of the data volume. The smallest contribution to this has the x-ray illumination, which is the only time the detector must be illuminated during the calibration data taking process, with about 30 min duration and 0.4 TB of data volume.

Analysis of the data on the DESY MAXWELL cluster using routines that have been optimized for parallel data processing can be performed in less than a day.

The electrical calibration described so far is essential for every experiment. In addition to this, the mechanical calibration, which determines the position of each pixel in all three dimensions and especially in relation to the beam, is of great importance for many experiments as well.

The detector mechanics is designed and built with high precision, but the movability of the four quadrants requires a quick, robust and accurate method of calibration for the position and orientation of the FEMs at any time. On top of the uncertainties introduced by the movements, small displacements and tilts of the FEMs[^3] with respect to ideal positions are unavoidable as their manufacture involves gluing steps which are naturally limited in their positioning accuracy.

A commonly used approach to calibrate the absolute position of each FEM in space is to take diffraction data of a well-known sample and fit the detector geometry .

The detector has just entered routine operation at the European XFEL and first crystal structures have been successfully determined . For these experiments the detector was calibrated using the procedures described here. However the procedures and recommended intervals are likely to evolve as the system gets used more often.

Currently, we recommend to take dark data for offset determination before and after each practical block of scientific data taking to account for any small drifts, e.g., of the temperature, that might occur during the scientific data taking process. Rechecking of the absolute gain and the gain ratios between the gain states using x-rays and the internal sources should be done periodically, as these might change over longer periods of time, e.g., due to accumulation of radiation damage.

While the system has been designed to be radiation hard, as detailed in earlier sections, forecasting the point at which radiation damage effects do show up is currently not possible for us.

Further, we recommend an electrical recalibration of the detector every time there has been a beam damage incident and if modules were exchanged. Mechanical position calibration is recommended every time the detector position is changed unless knowledge of the exact module positions is not required for the experiment.

Finding a quicker way to reliably determine high quality calibration constants is an ongoing development effort together with the detector group from the European XFEL.

# Performance data and imaging example

In order to test the noise performance over the entire dynamic range a pulsed IR laser was used, where the deposited energy in a pixel was varied between 1 and $10^4$ 12.4 keV photon equivalents by reducing the intensity of the IR pulses with calibrated filters.

The noise of the system is higher in medium and low gain mode (approximately equivalent to 3.5 keV and 18 keV, respectively), but it was shown that it always remained significantly below the Poisson noise of the incoming signal . The same IR laser was used to scan the dynamic range of the system, which was measured to be 34.4 x 10$^6$ electrons, corresponding to approximately 10$^4$ photons of 12.4 keV. The non-linearity of the low gain state proved to be better than 0.5% up to 5 x 10$^3$ photons .

Fig. a shows the x-ray image of a printed circuit board taken with the setup shown in Fig. b. The shown x-ray image is the average of 10,000 individual images.

Each individual image is corrected for pixel offset and gain on a per pixel basis. These corrections compensate fixed pattern 'zero-level' variations and pixel-to-pixel sensitivity variations that are commonly caused by process variations during the production of the ASIC, but can have many causes.

Some artifacts remain after the correction. Especially in the medium intensity region vertical stripes can be observed. These originate in the double sized pixels (400 m x 200 m) between ASICs. These pixels collect, on average, twice the amount photons, hence they appear brighter.

Fig. c shows the x-ray image of a pen drive. For this image 30,000 individual images were averaged.

Each image was offset and gain corrected as explained above, in addition a flat field correction was applied. The flat field correction accounts for the effective size of the pixels, removing the effect of the double size pixels, and removes artifacts from non-uniform illumination. Compensating the double sized pixels with a flat field correction is neither the only nor necessarily the best approach to correct for these pixels[^4]. The image not only shows the structure of the plastic cover of the pen drive, but also the sticky tape which was used to hold the pen drive to the acrylic glass in front of the FEM, demonstrating the high sensitivity of the system.

Lastly, the results of first user experiments show, that the detector is capable of determining protein crystal structures for both known and previously unknown proteins . Of course this success is possible in combination with all the other infrastructure of the beamline at the European XFEL.

# Summary

The adaptive gain integrating pixel detector, AGIPD, is an x-ray camera developed for use at the European XFEL. It was officially inaugurated together with the European XFEL in August 2017.

This paper reviewed the complete system of the AGIPD 1 million pixel camera currently installed at the SPB beamline of the European XFEL.

The system has a complex mechanical mounting that includes an in-vacuum movement system and many electronic boards, some of which are inside the vacuum vessel, some outside in an external housing. Its four independently movable quadrants can be arranged to form a horizontal slit or a rectangular hole with user selectable size.

The system is built from monolithic blocks of 2 x 8 ASICs forming a matrix of 128 x 512 pixels of (200 m)$^2$ size[^5]. Each pixel automatically adjusts to the incoming signal such that it can detect any number of photons from single photons to 10$^4$ photons of 12.4 keV above its noise floor. Images are stored in one of 352 memory cells during the XFEL pulse train and read out in between trains.

The detector noise is approximately 1 keV (750 eV for high CDS gain), which is sufficient to detect single photons in many experiments and most of the early experiments at the SPB beamline have used the AGIPD with great success .

The MID beamline of the European XFEL is scheduled to have a similar system installed in 2018.

The authors are deeply indebted to the seemingly countless number of people that have contributed to this project. Beyond the former team members that are not included on the author list these are the people 'behind the scenes' at DESY, PSI, Hamburg University, Bonn University and the European XFEL, who, with their tireless efforts, created a joint environment that made the development of this camera possible.

We would like to explicitly thank the European XFEL detector group and the involved XFEL groups CAS and ITDM. Integrating a system as complex as the AGIPD into the environment at the European XFEL was an enormous undertaking that required a huge collaborative effort and was successfully done in time to operate the detector at the inauguration of the European XFEL facility.

We are thankful for the help from many of the DESY infrastructure groups (esp. ZM1, ZM2, ZM3, ZE, FEA and FEB) during the development of the AGIPD.

We tested many prototypes, sometimes with help from external facilities. We would like to thank all of the people involved in these tests, which were crucial to the development process.

Last, but not least, we are very thankful for the valuable discussions with and the proof-reading of this manuscript by David Pennicard.

[^1]: The vacuum backplane forms a vacuum interface.

[^2]: The deviations from the ideal behavior at the beginning and end of each gain state are excluded from this fit.

[^3]: Since each FEM uses a monolithic silicon sensor which is defined by photolithography to precisions much better than 1 m the displacements of the pixels in each module w.r.t. their ideal positions in the pixel matrix is negligible compared to the other displacements described here.

[^4]: For most experiments at the European XFEL these pixels are currently excluded from the data analysis.

[^5]: If the double sized pixels are logically split the matrix is 128 x 526 pixels.


abstract: Recently there has been tremendous increase in the number of identified extra-solar planetary systems. Our understanding of their formation is tied to exoplanet internal structure models, which rely upon equations of state of light elements and compounds like water. Here we present shock compression data for water with unprecedented accuracy that shows water equations of state commonly used in planetary modeling significantly overestimate the compressibility at conditions relevant to planetary interiors. Furthermore, we show its behavior at these conditions, including reflectivity and isentropic response, is well described by a recent first-principles based equation of state. These findings advocate this water model be used as the standard for modeling Neptune, Uranus, and "hot Neptune" exoplanets, and should improve our understanding of these types of planets.
author: M.D. Knudson; M.P. Desjarlais; R.W. Lemke; T.R. Mattsson; M. French; N. Nettelmann; R. Redmer
date: 2024-10-01
title: Probing the interiors of the ice giants: Shock compression of water to 700 GPa and 3.8 g/cm$^3$

The past several years have seen a virtual explosion in the number of extra-solar planets discovered. Two rapidly growing populations of exoplanets are ice giants referred to as "hot Neptunes" and "mini-Neptunes";  planets roughly the same size as or, respectively, smaller than Neptune and Uranus that transit their host stars at significantly smaller radii, resulting in higher temperatures than the ice giants in our solar system. Understanding of the composition and formation of these planets, and thus development of these planetary systems, relies on our knowledge of the equation of state (EOS) of light elements and compounds like water, over a wide pressure and temperature range. To date much of the modeling of ice giants has employed the ANEOS  and Sesame  models for water that were developed decades ago . Discrepancies between these EOS models lead to significant differences in predicted radius evolution of Neptune-mass planets. Depending upon the total amount of heavy elements, and their distribution within the planetary interior, the resulting variation in predicted radius at a given age due to the water EOS can range between 5 and 30%  . This is a major factor in preventing accurate determination of exoplanet internal composition from their observed radius.

Recent quantum molecular dynamics (QMD) calculations of water  suggest an EOS that differs significantly from ANEOS and Sesame. Notably, when incorporated into planetary models, this first-principles (FP) based EOS predicts a $\sim$`<!-- -->`{=html}20% cooler core temperature for Neptune and Uranus . The conductivity properties of this FP model are also noteworthy , suggesting that water is super-ionic  at high densities, $\rho$, and low temperatures, $T$, relevant to planets such as Uranus and Neptune. This predicted property plays a key role in dynamo models to explain the enigmatic magnetic field structure of these planets . Another important result is derived from the predicted phase diagram of water: the icy giants Uranus and Neptune perhaps contain no "ice" but dissociated water at a high ionic conductivity, even less so would close-in exoplanets. Hot and mini-Neptunes may even comprise water plasma with substantial electronic conduction. However, the FP EOS for water has not been widely accepted due to its inability to reproduce results from laser driven shock wave experiments in the Mbar regime .

We present results of magnetically accelerated flyer-plate experiments on water performed at the Sandia Z machine , a pulsed power accelerator capable of producing extremely large current ($\sim$`<!-- -->`{=html}20 MA) and magnetic field densities ($\sim$`<!-- -->`{=html}10 MG) within a short circuit load. These data, in the range of 100-450 GPa along the Hugoniot – the locus of end states achievable through compression by large amplitude shock waves – have considerably higher precision than data obtained with previously used methods, and support the FP EOS for water. The high precision stems from the ability to perform well-defined flyer-plate experiments on Z; the magnetic pressure (\>500 GPa) can propel the outer anode to velocities approaching 30 km/s, enabling high-precision, plate-impact EOS measurements in the TPa regime . Furthermore, and more significantly, the present work obtained re-shock data of water in the range of 200-700 GPa. These data, at high $\rho$ and low $T$, provide a stringent test of the isentropic response of water in the several Mbar regime, which is directly relevant to the conditions of interest for planetary modeling of Neptune, Uranus , and presumably water-rich exoplanets such as the hot Neptune GJ436b . Finally, reflectivity on the Hugoniot was measured and compared to FP calculations for water .

An aluminum flyer-plate  was magnetically accelerated to peak velocities of 12-27 km/s across a 3-4 mm vacuum gap . The flyer-plate velocity was monitored throughout the entire trajectory using a Velocity Interferometer System for Any Reflector (VISAR ), at locations above and below an aluminum water cell . A rear quartz window in the cell provided optical access to the sample. In some cases an additional quartz plate was placed between the aluminum drive plate and the water sample, enabling data to be obtained using two different materials as the high-pressure standard, thereby increasing confidence in the measurements. Impact with the cell generated a strong, multi-Mbar shock wave in the aluminum drive plate. This shock was then transmitted either directly into the water sample, or into a quartz plate and then into the water sample. Upon reaching the rear quartz window, the shock was transmitted into the window and reflected back into the water,which re-shocked the water to a higher $P$ and $\rho$. In all cases the shock waves in the water and quartz were of sufficient amplitude that the resulting shocked material was reflecting , enabling the shock velocities to be directly measured using the VISAR. A total of 18 diagnostic channels were utilized for each experiment, enabling multiple, redundant measurements to be made, resulting in an overall uncertainty in the measured flyer-plate and shock velocities of a few tenths of a percent .

The shocked state of the water was determined using the impedance matching technique and the Rankine-Hugoniot (RH) jump relations , a set of conditions derived by considering conservation of mass, momentum, and energy across a steady propagating shock wave. The shocked state of the aluminum (quartz) drive plate was determined from the known Hugoniot of aluminum  (quartz ) and the measured flyer-plate (quartz shock) velocity; this defined a point in the pressure - particle velocity ($P-u_p$) plane, as shown in Fig. . When the shock transits into the water, a release wave propagates back toward the flyer-plate, and thus the state of the drive plate is constrained to lie on a release adiabat from this point in the $P-u_p$ plane, shown in Fig.  as the green line. The shocked state of the water is constrained to lie along a chord in the $P-u_p$ plane with slope given by the product of the measured shock velocity of water, $U_{sw}$, and the known initial density. The intersection of these two curves provides $P$ and $u_p$, shown in Fig.  as $(P_1,u_{p1})$; The RH jump relations then provide $\rho$ in the shocked state. Uncertainties in all kinematic values were determined through a Monte Carlo technique, which uses a statistical process for propagation of all random measurement errors and systematic errors in the standards . Using this technique, the one-sigma uncertainties in $P$ and $\rho$ were found to be  0.5% and  1%, respectively.

A total of 8 Hugoniot experiments were performed over the range of 100 to 450 GPa. Results of these experiments are shown as the red symbols in Fig. (a). Also shown are Hugoniot data of Mitchell and Nellis , Volkov *et al.* , Celliers *et al.* , and Podurets *et al.* , and the predicted Hugoniot response from ANEOS , Sesame 7150 , and the recent FP EOS model of French *et al.* . Note that a reanalysis of the nuclear driven datum of Podurets *et al.*, using an improved aluminum standard for impedance matching , resulted in a slight decrease in $\rho$. The low-$P$ end of our data is in good agreement with the gas gun data of Mitchell and Nellis and the explosively driven shock data of Volkov *et al.* In contrast, our data are significantly less compressible than the laser driven data of Celliers *et al.*, which tend to support the much more compressible ANEOS and Sesame Hugoniots, albeit with significantly large uncertainty and scatter. The vastly reduced uncertainty in $\rho$ for our data, roughly an order of magnitude, strongly suggest that water is much less compressible than the ANEOS and Sesame models predict, and that water is instead very accurately described by the FP EOS of French *et al.* Furthermore, the reanalyzed Podurets *et al.* datum is also in very good agreement with the FP EOS. Thus, with the exception of the Celliers *et al.* data, the FP based model for water matches all experimental Hugoniot data up to  1.4 TPa.

In all 8 of the Hugoniot experiments described above, the reflected shock from the rear quartz window drove the water from the Hugoniot state to a re-shocked state at higher $P$ and $\rho$. The measured shock velocity in the water immediately prior to reflection from the rear quartz window defined the initial shocked state of the water. The measured shock velocity in the rear quartz window and the known Hugoniot of quartz provided the double-shocked $P$ and $u_p$ for water, shown in Fig.  as $(P_2,u_{p2})$. The velocity of the second shock in the water, $U_{sw2}$, was then determined by the RH jump relations using the change in $P$, $(P_2-P_1 )$, and $u_p$, $(u_{p2}-u_{p1} )$. The re-shock $\rho$ was then determined from $U_{sw2}$, the first shock $\rho$, and $(u_{p2}-u_{p1})$. Using the Monte Carlo technique, the one-sigma uncertainties in $P$ and $\rho$ for the re-shock states were found to be  0.5-1% and  1-2%, respectively. Although the uncertainty for the re-shock data is larger than that for the principal Hugoniot data (entirely due to the larger uncertainty in the initial state), the accuracy of the present data is a significant improvement over previous re-shock data of Mitchell and Nellis  (uncertainty in $\rho$ of 4-14%) and the pre-compressed Hugoniot data of Lee *et al.* , (uncertainty in $\rho$ of 5-10%).

The re-shock data for water are shown in Fig. (b), where first and second shock states are correlated by like symbols. Also shown are several FP re-shock Hugoniots (thin red lines) and isentropes (thin black lines) for comparison . These re-shock Hugoniots along with the known Hugoniot of quartz  were used to determine the double-shock envelopes –- the locus of end states achievable through shock and re-shock using a quartz anvil: FP (orange line), ANEOS  (pink line), and Sesame  (gray line). These re-shock data further confirm the less compressible response of water above  100 GPa.

Note that the FP re-shock Hugoniots (red) and isentropes (black) are nearly coincident over the $\rho$ range accessed through the re-shock experiments. This is due to a second order contact for the Hugoniot and isentrope at the initial state , which is most easily seen by expanding the entropy as a function of volume in a Taylor series. This implies that the Hugoniot and isentrope are very close in $P$ and $\rho$ until, at large compression, the rise in $T$ associated with the irreversible shock becomes large enough that thermal pressures become significant. In the range investigated in this study, the difference in $T$ between the re-shock Hugoniot states and the isentrope at the re-shock $\rho$, as determined by the FP EOS , ranged from 200K (out of 6800K) to 330K (out of 40000K) at the lowest and highest $P$, respectively. This makes such a re-shock measurement the best possible test of the isentropic response of the EOS model in this range of $P$ and $\rho$. Thus the present data validates the isentropic response of the FP EOS in the $P$ and $\rho$ regime that is intersected by the water-rich models of Neptune and Uranus , shown in green, and the exoplanet GJ436b , shown in blue.

The VISAR was also used to infer reflectivity, $R$, of water (at 532 nm) along the Hugoniot. A quadrature VISAR was used for all experiments, which provides four measures of the interference signal at $90\,^{\circ}$ intervals. The signals at $180\,^{\circ}$ intervals can be subtracted, ensuring the remaining signal only includes coherent reflected laser light (incoherent light, such as self-emission from the hot plasma, would equally contribute to all four quadrature signals). Comparison of the magnitude of these subtracted signals before and after shock breakout from the water to the quartz rear window provides a relative measure of the shocked water $R$ with respect to shocked quartz . The uncertainty in $R$ was taken to be the linear sum of the standard deviation of the inferred $R$ from the nine independent VISAR signals obtained from each water cell and the reported uncertainty in $R$ of shocked quartz .

$R$ data along the Hugoniot are shown in Fig. . Also shown are data from Celliers *et al.*  and the predicted $R$ from FP calculations of French and Redmer  using both the Perdew, Burke, and Ernzerhof (PBE) and Heyd, Scuseria, and Ernzerhof (HSE) functionals for exchange and correlation. It was anticipated that the HSE functional, which includes the nonlocal Fock exchange, would prove to be more accurate in the calculation of $R$, as this functional has been shown to better reproduce the band gap in semiconductor materials (PBE is known to significantly underestimate the band gap). In comparison with $R$ data of Celliers *et al.*  it would appear that the HSE calculations are less accurate. However, our data suggest a much lower peak $R$, which is in significantly better agreement with the HSE calculations. We note that two recent data points ($\sim$`<!-- -->`{=html}140 and 260 GPa) from a group  at the Gekko laser in Japan also suggest lower $R$, in very good agreement with our results. These new results lend confidence to the FP calculations, which also predict a super-ionic phase of water at low $T$ and high $\rho$ conditions relevant to planetary interiors. Furthermore, these results strongly suggest that at these conditions water is in a plasma phase, which would imply that a $T$=0 K EOS for water is not sufficient for modeling of hot and mini-Neptunes, and that water would be expected to mix in the H/He envelope rather than form an ice shell separate from an outer H/He envelope.

We presented data with unprecedented accuracy for shock compression of water to 0.7 TPa and 3.8 g/cc in a regime relevant to water-rich models of Uranus, Neptune and the exoplanet GJ436b. The experimental $P$, $\rho$, and $R$ are in excellent agreement with density functional theory predictions, thereby validating first-principles thermodynamic calculations as a sound basis for planetary modeling, and strongly advocating the FP EOS be the standard in modeling water in Neptune, Uranus, and "hot Neptune" exoplanets. In particular this work supports the prediction of a $\sim$`<!-- -->`{=html}20% cooler core temperature for Neptune and Uranus . As the calculated amount of H and He in the planets decreases with the stiffness of the water EOS, confidence in the presence of a few percent H and He in the deep interior of Neptune and Uranus, as derived from the (rather stiff) FP EOS based models  , is strengthened by this work. As H would be metallic, this might influence the generation of the magnetic field. Furthermore, the validation of the FP EOS in the regime relevant to planetary interiors all but eliminates one significant source of uncertainty in the predicted radius evolution of Neptune-mass planets within assumed composition models. This will improve our understanding of the interior structure of these planets, and perhaps our understanding of these planetary systems.

We acknowledge the crew of the Sandia Z facility for their contributions to these experiments, AB and MB for assistance in numerical calculations, and support from the DFG via the SFB 652 and the grant Re 882/11-1. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.