
LUXZEPLIN (LZ) is a next generation dark matter direct detection experiment
that will operate 4850 feet underground at the Sanford Underground Research
Facility (SURF) in Lead, South Dakota, USA. Using a twophase xenon detector
with an active mass of 7 tonnes, LZ will search primarily for lowenergy
interactions with Weakly Interacting Massive Particles (WIMPs), which are
hypothesized to make up the dark matter in our galactic halo. In this paper,
the projected WIMP sensitivity of LZ is presented based on the latest
background estimates and simulations of the detector. For a 1000 live day run
using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence
level spinindependent WIMPnucleon cross sections above $1.6 \times 10^{48}$
cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$
discovery potential is projected reaching cross sections below the existing and
projected exclusion limits of similar experiments that are currently operating.
For spindependent WIMPneutron(proton) scattering, a sensitivity of $2.7
\times 10^{43}$ cm$^{2}$ ($8.1 \times 10^{42}$ cm$^{2}$) for a 40
$\mathrm{GeV}/c^{2}$ mass WIMP is expected. With construction well underway, LZ
is on track for underground installation at SURF in 2019 and will start
collecting data in 2020.

Momentumresolved photoelectron emission from xenon in colinearly polarized
twocolor laser fields at abovethreshold ionization conditions is studied both
experimentally and theoretically. We utilize phaseofthephase spectroscopy as
recently introduced by Skruszewicz et al., Phys. Rev. Lett. 115, 043001 (2015)
to analyze the dependence of the yields on the relative phase $\varphi$ between
the fundamental and second harmonic laser fields. The resulting phaseofphase
spectra feature a characteristic checkerboard pattern, which can analytically
be described within the strongfield approximation.

Nondipole effects in strongfield photoelectron momentum spectra have been
revealed experimentally [C.T.L. Smeenk et al., Phys. Rev. Lett. 106, 193002
(2011); A. Ludwig et al., Phys. Rev. Lett. 113, 243001 (2014)]. For certain
laser parameters and photoelectron momenta the spectra were found to be shifted
against the laser propagation direction whereas one would naively assume that
the radiation pressure due to the $\vec{v}\times\vec{B}$force pushes electrons
always in propagation direction. Only the interplay between Lorentz and Coulomb
force may give rise to such counterintuitive dynamics. In this work, we
calculate the momentumdependent shift in and against the propagation direction
by extending the quantum trajectorybased Coulombcorrected strongfield
approximation beyond the dipole approximation. A semianalytical treatment
where both magnetic and Coulomb force are treated perturbatively but
simultaneously reproduces the results from the full numerical solution of the
equations of motion.

The $N$particle wavefunction has too many dimensions for a direct time
propagation of a manybody system according to the timedependent Schr\"odinger
equation (TDSE). On the other hand, timedependent density functional theory
(TDDFT) tells us that the singleparticle density is, in principle, sufficient.
However, a practicable equation of motion (EOM) for the accurate time evolution
of the singleparticle density is unknown. It is thus an obvious idea to
propagate a quantity which is not as reduced as the singleparticle density but
less dimensional than the $N$body wavefunction. Recently, we have introduced
timedependent renormalizednaturalorbital theory (TDRNOT). TDRNOT is based on
the propagation of the eigenfunctions of the onebody reduced density matrix
(1RDM), the socalled natural orbitals. In this paper we demonstrate how
TDRNOT is related to the multiconfigurational timedependent HartreeFock
(MCTDHF) approach. We also compare the performance of MCTDHF and TDRNOT vs the
TDSE for singlephoton double ionization (SPDI) of a 1D helium model atom. SPDI
is one of the effects where TDDFT does not work in practice, especially if one
is interested in correlated photoelectron spectra, for which no explicit
density functional is known.

In this Technical Design Report (TDR) we describe the LZ detector to be built
at the Sanford Underground Research Facility (SURF). The LZ dark matter
experiment is designed to achieve sensitivity to a WIMPnucleon
spinindependent cross section of three times ten to the negative fortyeighth
square centimeters.

Photoelectron spectra obtained from the ab initio solution of the
timedependent Schr\"odinger equation can be in striking disagreement with
predictions by the strongfield approximation (SFA) not only at low energy but
also around twice the ponderomotive energy where the transition from the direct
to the rescattered electrons is expected. In fact, the relative enhancement of
the ionization probability compared to the SFA in this regime can be several
orders of magnitude. We show for which laser and target parameters such an
enhancement occurs and for which the SFA prediction is qualitatively good. The
enhancement is analyzed in terms of the Coulombcorrected action along analytic
quantum orbits in the complextime plane, taking soft recollisions under the
Coulomb barrier into account. These recollisions in complex time and space
prevent a separation into subbarrier motion up to the "tunnel exit" and
subsequent classical dynamics. Instead, the entire quantum path up to the
detector determines the ionization probability.

Recently introduced timedependent renormalizednatural orbital theory
(TDRNOT) is extended towards a multicomponent approach in order to describe
H$_2^+$ beyond the BornOppenheimer approximation. Two kinds of natural
orbitals, describing the electronic and the nuclear degrees of freedom are
introduced, and the exact equations of motion for them are derived. The theory
is benchmarked by comparing numerically exact results of the timedependent
Schr\"odinger equation for a H$_2^+$ model system with the corresponding TDRNOT
predictions. Ground state properties, linear response spectra, fragmentation,
and highorder harmonic generation are investigated.

We measure the quantum capacitance and probe thus directly the electronic
density of states of the high mobility, Dirac type of twodimensional electron
system, which forms on the surface of strained HgTe. Here we show that observed
magnetocapacitance oscillations probe, in contrast to magnetotransport,
primarily the top surface. Capacitance measurements constitute thus a powerful
tool to probe only one topological surface and to reconstruct its Landau level
spectrum for different positions of the Fermi energy.

Joint electronion energy spectra for the dissociative ionization of a model
H$_2^+$ in fewcycle, infrared laser pulses are calculated via the numerical ab
initio solution of the timedependent Schr\"odinger equation. A strong,
pulsedependent modulation of the ionization probability for certain values of
the protons' kinetic energy (but almost independent of the electron's energy)
is observed. With the help of models with frozen ions, this featurewhich
mistakenly might be attributed to vibrational excitationsis traced back to
the transient population of electronically excited states, followed by
ionization. This assertion is further corroborated employing a twolevel model
incorporating strongfield ionization from the excited state.

Recently introduced timedependent renormalizednaturalorbital theory
(TDRNOT) is based on the equations of motion for the socalled natural
orbitals, i.e., the eigenfunctions of the onebody reduced density matrix.
Exact TDRNOT can be formulated for any timedependent twoelectron system in
either spin configuration. In this paper, the method is tested against
highorder harmonic generation (HHG) and Fano profiles in absorption spectra
with the help of a numerically exactly solvable onedimensional model He atom,
starting from the spinsinglet ground state. Such benchmarks are challenging
because Fano profiles originate from transitions involving autoionizing states,
and HHG is a strongfield phenomenon well beyond linear response. TDRNOT with
just one natural orbital per spin in the helium spinsinglet case is equivalent
to timedependent HartreeFock or timedependent density functional theory
(TDDFT) in exact exchangeonly approximation. It is not unexpected that TDDFT
fails in reproducing Fano profiles due to the lack of doubly excited,
autoionizing states. HHG spectra, on the other hand, are widely believed to be
wellcaptured by TDDFT. However, HHG spectra of helium may display a second
plateau that originates from simultaneous HHG in He$^+$ and neutral He. It is
found that already TDRNOT with two natural orbitals per spin is sufficient to
capture this effect as well as the Fano profiles on a qualitative level. With
more natural orbitals (68 per spin) quantitative agreement can be reached.
Errors due to the truncation to a finite number of orbitals are identified.

The presence of a weak secondharmonic field in an intenselaser ionization
experiment affects the momentumresolved electron yield, depending on the
relative phase between the $\omega$ and the $2\omega$ component. The proposed
twocolor 'phaseofthephase spectroscopy' quantifies for each final electron
momentum a relativephase contrast (RPC) and a 'phase of the phase' (PP),
describing how much and with which phase lag, respectively, the yield changes
as function of relative phase. Experimental results for RPC and PP spectra for
rare gas atoms and CO$_2$ are presented. The spectra demonstrate a rather
universal structure that is analyzed with the help of a simple model based on
electron trajectories, wavepacket spreading, and (multiple) rescattering.
Details in the PP and RPC spectra are target sensitive and thus may be used to
extract structural (or even dynamical) information with high accuracy.

Recently introduced timedependent renormalized natural orbital theory
(TDRNOT) is tested on nonsequential double ionization (NSDI) of a numerically
exactly solvable onedimensional model He atom subject to fewcycle, 800nm
laser pulses. NSDI of atoms in strong laser fields is a prime example of
nonperturbative, highly correlated electron dynamics. As such, NSDI is an
important "worstcase" benchmark for any timedependent few and manybody
technique beyond linear response. It is found that TDRNOT reproduces the
celebrated NSDI "knee," i.e., a manyorderofmagnitude enhancement of the
double ionization yield (as compared to purely sequential ionization) with only
the ten most significant natural orbitals (NOs) per spin. Correlated
photoelectron spectra  as "more differential" observables  require more NOs.

Natural orbital theory is a computationally useful approach to the few and
manybody quantum problem. While natural orbitals are known and applied since
many years in electronic structure applications, their potential for
timedependent problems is being investigated only since recently. Correlated
twoparticle systems are of particular importance because the structure of the
twobody reduced density matrix expanded in natural orbitals is known exactly
in this case. However, in the timedependent case the natural orbitals carry
timedependent phases that allow for certain timedependent gauge
transformations of the first kind. Different phase conventions will, in
general, lead to different equations of motion for the natural orbitals. A
particular phase choice allows us to derive the exact equations of motion for
the natural orbitals of any (laser) driven twoelectron system explicitly,
i.e., without any dependence on quantities that, in practice, require further
approximations. For illustration, we solve the equations of motion for a model
helium system. Besides calculating the spinsinglet and spintriplet ground
states, we show that the linear response spectra and the results for resonant
Rabi flopping are in excellent agreement with the benchmark results obtained
from the exact solution of the timedependent Schr\"odinger equation.

We survey sufficient degree conditions, for a variety of graph properties,
that are best possible in the same sense that Chvatal's wellknown degree
condition for hamiltonicity is best possible.

The influence of inner electrons on the ionization dynamics in strong laser
fields is investigated in a wavelength regime where the inner electron dynamics
is usually assumed to be negligible. The role of inner electrons is of
particular interest for the application of frozencore approximations and
pseudopotentials in timedependent density functional theory (TDDFT) and the
singleactiveelectron (SAE) approximation in strongfield laser physics.
Results of TDDFT and SAE calculations are compared with exact ones obtained by
the numerical ab initio solution of the threeelectron timedependent
Schr\"odinger equation for a lithium model atom. It is found that dynamical
antiscreening, i.e., a particular form of dynamical core polarization, may
substantially alter the ionization rate in the singlephoton regime.
Requirements for the validity of the approximations in the single and
multiphoton ionization domain are identified.

The strong field approximation (SFA) formulated in terms of socalled
"quantum orbits" led to much insight into intenselaser driven ionization
dynamics. In plain SFA, the emitted electron is treated as a free electron in
the laser field alone. However, with improving experimental techniques and more
advanced numerical simulations it becomes more and more obvious that the plain
SFA misses interesting effects even on a qualitative level. Examples are
holographic side lobes, the lowenergy structure, radial patterns in
photoelectron spectra at low kinetic energies, and strongly rotated angular
distributions. For this reason increasing effort has been recently devoted to
Coulomb corrections of the SFA. In the current paper, we follow a similar line
but consider ionization of metal clusters. It is known that photoelectrons from
clusters can be much more energetic than those emitted from atoms or small
molecules, especially if the Mie resonance of the expanding cluster is evoked.
We develop a SFA that takes the collective field inside the cluster via the
simple rigidsphere model into account. Our approach is based on
fieldcorrected quantum orbits so that the acceleration process (or any other
spectral feature of interest) can be investigated in detail.

Favorably scaling numerical timedependent manyelectron techniques such as
timedependent density functional theory (TDDFT) with adiabatic
exchangecorrelation potentials typically fail in capturing highly correlated
electron dynamics. We propose a method based on natural orbitals, i.e., the
eigenfunctions of the onebody reduced density matrix, that is almost as
inexpensive numerically as adiabatic TDDFT, but which is capable of describing
correlated phenomena such as doubly excited states, autoionization, Fano
profiles in the photoelectron spectra, and strongfield ionization in general.
Equations of motion (EOM) for natural orbitals and their occupation numbers
have been derived earlier. We show that by using renormalized natural orbitals
(RNO) both can be combined into one equation governed by a hermitian effective
Hamiltonian. We specialize on the twoelectron spinsinglet system, known as
being a "worst case" testing ground for TDDFT, and employ the widely used,
numerically exactly solvable, onedimensional helium model atom (in a laser
field) to benchmark our approach. The solution of the full, nonlinear EOM for
the RNO is plagued by instabilities, and resorting to linear response is not an
option for the ultimate goal to study nonperturbative dynamics in intense laser
fields. We therefore make two rather bold approximations: we employ the
initialstate"frozen" effective RNO Hamiltonian for the timepropagation and
truncate the number of RNO to only two per spin. Surprisingly, it turns out
that even with these strong approximations we obtain a highly accurate ground
state, reproduce doublyexcited states, and autoionization.

As part of the Snowmass process, the Cosmic Frontier WIMP Direct Detection
subgroup (CF1) has drawn on input from the Cosmic Frontier and the broader
Particle Physics community to produce this document. The charge to CF1 was (a)
to summarize the current status and projected sensitivity of WIMP direct
detection experiments worldwide, (b) motivate WIMP dark matter searches over a
broad parameter space by examining a spectrum of WIMP models, (c) establish a
community consensus on the type of experimental program required to explore
that parameter space, and (d) identify the common infrastructure required to
practically meet those goals.

We study the ionization dynamics in intense laserdroplet interaction using
threedimensional, relativistic particleincell simulations. Of particular
interest is the laser intensity and frequency regime for which initially
transparent, wavelengthsized targets are not homogeneously ionized. Instead,
the charge distribution changes both in space and in time on a subcycle scale.
One may call this the extreme nonlinear Mieoptics regime. We find that 
despite the fact that the plasma created at the droplet surface is overdense 
oscillating electric fields may penetrate into the droplet under a certain
angle, ionize, and propagate in the just generated plasma. This effect can be
attributed to the local field enhancements at the droplet surface predicted by
standard Mie theory. The penetration of the fields into the droplet leads to
the formation of a highly inhomogeneous charge density distribution in the
droplet interior, concentrated mostly in the polarization plane. We present a
selfsimilar, exponential fit of the fractional ionization degree which depends
only on a dimensionless combination of electric field amplitude, droplet
radius, and plasma frequency with only a weak dependence on the laser frequency
in the overdense regime.

The Floquet theorem allows to reformulate periodic timedependent problems
such as the interaction of a manybody system with a laser field in terms of
timeindependent, fielddressed states, also known as Floquet states. If this
was possible for density functional theory as well, one could reduce in such
cases timedependent density functional theory to a timeindependent Floquet
density functional theory. We analyze under which conditions the Floquet
theorem is applicable in a densityfunctional framework. By employing numerical
{\em ab initio} solutions of the interacting timedependent Schr\"odinger
equation with timeperiodic external potentials we show that the exact
effective potential in the corresponding KohnSham equation is {\em not}
unconditionally periodic. Whenever several Floquet states in the interacting
system are involved in a physical process the corresponding
Hartreeexchangecorrelation potential is not periodic with the external
frequency only. Using an analytically solvable example we demonstrate that, in
general, the periodicity of the timedependent KohnSham Hamiltonian cannot be
restored by choosing a different initial state. Only if the external periodic
potential is sufficiently weak such that the initial state of the interacting
system evolves adiabatically to a single, fielddressed state, the resulting
KohnSham system admits the application of the Floquet theorem.

We use a quantum trajectorybased semiclassical method to account for
Coulomb interaction between the photoelectron and the parent ion in the
classically forbidden, subbarrier region during strongfield tunneling
ionization processes. We show thatbesides the wellknown modification of the
tunneling ionization probabilitythere is also an influence on the
interference pattern in the photoelectron spectra. In the longwavelength
limit, the shift of the intracycle interference fringes caused by subbarrier
Coulomb effects in the laser polarization direction can be derived
analytically. We compare our results with \emph{ab initio} solutions of the
timedependent Schr\"odinger equation and find good agreement in the
longwavelength regime, whereas the standard strong field approximation fails.
We show that the nodal structure along loworder abovethreshold ionization
rings is also affected by subbarrier Coulomb effects.

The screening of a 2p corehole in Na clusters is investigated using density
functional theory applied to an extended jellium model with an allelectron
atom in the center. The study is related to recent experiments at the free
electron laser at DESY in which photoelectron spectra from massselected,
coreshell ionized metal clusters have been recorded. Relaxed and unrelaxed
binding energies as well as KohnSham orbital energies are calculated in
PerdewZunger selfinteractioncorrected exchangeonly local spindensity
approximation for valence and 2p core electrons in Na clusters up to 58 atoms.
The relaxed binding energies follow approximately the metalsphere behavior.
The same behavior is seen in the experiment for sufficiently big clusters,
indicating perfect screening and that the relaxation energy due to screening
goes to the photoelectron. Instead, calculating the kinetic energy of the
photoelectrons using unrelaxed binding energies or KohnSham orbital energies
yields wrong results for coreshell electrons. The screening {\em dynamics} are
investigated using timedependent density functional theory. It is shown that
screening occurs on two time scales, a coreshelldependent inneratomic and an
interatomic valence electron time scale. In the case of Na 2p ionization the
remaining electrons in the 2p shell screen within tens of attoseconds while the
screening due to cluster valence electrons occurs within several hundreds of
attoseconds. The screening timescales may be compared to the photon energy and
cluster sizedependent escape times of the photoelectron in order to estimate
whether the photoelectron is capable of picking up the relaxation energy or
whether the residual system is left in an excited state.

We propose a method to obtain Floquet statesalso known as lightinduced
statesand their quasienergies from realtime wavefunctions without solving
the Floquet equation. This is useful for the analysis of various phenomena in
timedependent quantum dynamics if the Hamiltonian is not strictly periodic, as
in short laser pulses, for instance. There, the population of the Floquet
states depends on the pulse form and is automatically contained in the
realtime wavefunction but not in the standard Floquet approach. Several
examples in the area of intense laseratom interaction are exemplarily
discussed: (i) the observation of even harmonics for an inversionsymmetric
potential with a single bound state; (ii) the dependence of the population of
Floquet states on (gauge) transformations and the emergence of an invariant,
observable photoelectron spectrum; (iii) the driving of resonant transitions
between dressed states, i.e., the dressing of dressed states, and (iv) spectral
enhancements at channel closings due to the ponderomotive shift of
abovethreshold ionization peaks.

We prove a generalization of the van Leeuwen theorem towards quantum
electrodynamics, providing the formal foundations of a timedependent KohnSham
construction for coupled quantized matter and electromagnetic fields. Thereby
we circumvent the symmetrycausality problems associated with the
actionfunctional approach to KohnSham systems. We show that the effective
external fourpotential and fourcurrent of the KohnSham system are uniquely
defined and that the effective fourcurrent takes a very simple form. Further
we rederive the RungeGross theorem for quantum electrodynamics.

The effects of intense electromagnetic fields on the decay of quasistationary
states are investigated theoretically. We focus on the parameter regime of
strong laser fields and nonlinear effects where an essentially nonperturbative
description is required. Our approach is based on the imaginary time method
previously introduced in the theory of strongfield ionization. Spectra and
total decay rates are presented for a test case and the results are compared
with exact numerical calculations. The potential of this method is confirmed by
good quantitative agreement with the numerical results.