
Electromagnetic radiation at terahertz (THz) frequencies is a useful tool in
probing and controlling matter and light in new and interesting ways,
especially at high peakfields and pulse energies. Generating THz radiation
often employs nonlinear optical processes, for which the overlapping of
stretched, broadband nearinfrared (NIR) pulse copies within nonlinear crystals
is common. Here we show that for narrowband THz generation, the higherorder
phase present on the NIR pulses offers control of the properties of the THz,
for example creating temporally complex THz with multiplexed NIR pulses. We
manipulate the phase of two NIR pump pulses independently to remove higher
order effects and generate record mJlevel THz in two crystals simultaneously,
with an average total energy of 604 microjoule at 361 GHz with 1% bandwidth.
This high pulse energy combined with such a narrow bandwidth has broad
implications for accelerator applications, resonant driven material studies,
and nonlinear THz spectroscopy.

We perform a detailed experimental and theoretical study of the transition
from multiphoton to opticalfield photoemission from ndoped, singlecrystal
silicon nanotips. Around this transition, we measure an enhanced emission rate
as well as intensitydependent structure in the photoelectron yield from the
illuminated nanostructures. Numerically solving the timedependent
Schr\"odinger equation (TDSE), we demonstrate that the excess emission derives
from the buildup of standing electronic wavepackets near the surface of the
silicon, and the intensity dependent structure in this transition results from
the increased ponderomotive potential and channel closing effects. By way of
timedependent perturbation theory (TDPT), we then show that the visibility of
intensity dependent structure, the transition rate from multiphoton to
opticalfield emission, and scaling rate at high intensities are all consistent
with a narrow band of groundstate energies near the conduction band dominating
the emission process in silicon. These results highlight the importance of
considering both coherent electron wavepacket dynamics at the emitter surface
as well as the groundstate energy distribution when interpreting strongfield
photoemission from solids.

A velocitymapimaging spectrometer is demonstrated to characterize the
normalized transverse emittance of photoemitted electron bunches. The
twodimensional (2D) projected velocity distribution images of photoemitted
electrons are recorded by the detection system and analyzed to obtain the
normalized transverse emittance. With the presented distribution function of
the electron photoemission angles a mathematical method is implemented to
reconstruct the threedimensional (3D) velocity distribution curve. As a first
example, multiphoton emission from a planar Au surface is studied via
irradiation at a glancing angle by intense 45 fs laser pulses at a central
wavelength of 800 nm. The reconstructed energy distribution agrees very well
with the BerglundSpicer theory of photoemission. The normalized transverse
emittance of the intrinsic electron bunch is characterized to be 0.52 and 0.05
$\pi \cdot mm \cdot mrad$ in $X$ and $Y$directions, respectively.

The strong ellipticity dependence of highharmonic generation in gases
enables numerous experimental techniques that are nowadays routinely used, for
instance, to create isolated attosecond pulses. Extending such techniques to
highharmonic generation in solids requires a fundamental understanding of the
microscopic mechanism of the highharmonic generation. Here, using extensive
firstprinciples simulations within a timedependent densityfunctional
framework, we show how intraband and interband mechanisms are strongly and
differently affected by the ellipticity of the driving laser field. The complex
interplay between intraband and interband effects can be used to tune and
improve harmonic emission in solids. In particular, we show that the energy
cutoff of the highharmonic plateau can be increased by as much as 30\% using a
finite ellipticity of the driving field, opening a new avenue for better
understanding and control of HHG in solids based on ellipticity. Also, we
demonstrate the possibility to generate, from a single circularly polarized
driving field, circularly polarized harmonics with alternating helicity. Our
work shows that ellipticity provides an additional knob to experimentally
control highorder harmonic generation in solids.

Highenergy, carrierenvelope phase (CEP)stable, subcycle, midinfrared
(midIR) pulses can provide unique opportunities of exploring phasesensitive
strongfield lightmatter interactions in atoms, molecules, and solids. In the
midIR wavelength, the ponderomotive energy of laser pulses is dramatically
increased (versus the visible/nearinfrared) and, therefore, the Keldysh
parameter is much smaller than unity even at relatively modest laser
intensities. This enables to study the subcycle electron dynamics in solids
via highharmonic generation (HHG) without damage. One can also control the
electron emissions from nanodevices in the subcycle time scale. These efforts
are opening a great opportunity towards petahertz electronics. Here, we present
a highenergy, subcycle pulse synthesizer based on a midIR optical parametric
amplifier (OPA), pumped by CEPstable, 2.1 um femtosecond pulses, and its
application to HHG in solids. The signal and idler combined spectrum spans from
2.5 to 9.0 um, which covers the whole midwaveinfrared (MWIR) region. We
coherently synthesize the passively CEPstable fewcycle signal and idler
pulses to generate 33 uJ, 0.88cycle (12.4 fs), multiGW pulses centered at
~4.2 um, which is further energy scalable. The inline synthesis of the
CEPstable subcycle pulse is realized through the typeI collinear OPA with
minimal temporal walkoff. The MWIR subcycle pulse is used for driving HHG in
thin silicon samples, producing harmonics up to ~19th order with a continuous
spectral coverage due to the isolated emission by the subcycle driver. Our
demonstration offers an energy scalable and technically simple platform of
laser sources generating CEPstable subcycle pulses in the whole MWIR region
for investigating isolated phasesensitive strongfield interactions in solids
and gases.

Free Electron Lasers (FELs) are a solution for providing intense, coherent
and bright radiation in the hard Xray regime. Due to the low wallplug
efficiency of FEL facilities, it is crucial and additionally very useful to
develop complete and accurate simulation tools for better optimizing a FEL
interaction. The highly sophisticated dynamics involved in a FEL process was
the main obstacle hindering the development of general simulation tools for
this problem. We present a numerical algorithm based on finite difference time
domain/Particle in cell (FDTD/PIC) in a Lorentz boosted coordinate system which
is able to fulfil a fullwave simulation of a FEL process. The developed
software offers a suitable tool for the analysis of FEL interactions without
considering any of the usual approximations. A coordinate transformation to
bunch rest frame makes the very different length scales of bunch size, optical
wavelengths and the undulator period transform to values with the same order.
Consequently, FDTD/PIC simulations in conjunction with efficient
parallelization techniques make the fullwave simulation feasible using the
available computational resources. Several examples of free electron lasers are
analyzed using the developed software, the results are benchmarked based on
standard FEL codes and discussed in detail.

Linear acceleration in free space is a topic that has been studied for over
20 years, and its ability to eventually produce highquality, high energy
multiparticle bunches has remained a subject of great interest. Arguments can
certainly be made that such an ability is very doubtful. Nevertheless, we chose
to develop an accurate and truly predictive theoretical formalism to explore
this remote possibility in a computational experiment. The formalism includes
exact treatment of Maxwell's equations, exact relativistic treatment of the
interaction among the multiple individual particles, and exact treatment of the
interaction at near and far field. Several surprising results emerged. For
example, we find that 30 keV electrons (2.5% energy spread) can be accelerated
to 7.7 MeV (2.5% spread) and to 205 MeV (0.25% spread) using 25 mJ and 2.5 J
lasers respectively. These findings should hopefully guide and help develop
compact, highquality, ultrarelativistic electron sources, avoiding
conventional limits imposed by material breakdown or structural constraints.

An accurate analytic model describing highharmonic generation (HHG) in
solids is derived. Extensive firstprinciples simulations within a
timedependent densityfunctional framework corroborate the conclusions of the
model. Our results reveal that: (i) the emitted HHG spectra are highly
anisotropic and laserpolarization dependent even for cubic crystals, (ii) the
harmonic emission is enhanced by the inhomogeneity of the electronnuclei
potential, the yield is increased for heavier atoms, and (iii) the cutoff
photon energy is driverwavelength independent. Moreover, we show that it is
possible to predict the laser polarization for optimal HHG in bulk crystals
solely from the knowledge of their electronic band structure. Our results pave
the way to better control and optimize HHG in solids by engineering their band
structure.

Ultrashort electron beams with narrow energy spread, high charge, and low
jitter are essential for resolving phase transitions in metals, semiconductors,
and molecular crystals. These semirelativistic beams, produced by
phototriggered electron guns, are also injected into accelerators for xray
light sources. The achievable resolution of these timeresolved electron
diffraction or xray experiments has been hindered by surface field and timing
jitter limitations in conventional RF guns, which thus far are <200 MV/m and
>96 fs, respectively. A gun driven by opticallygenerated singlecycle THz
pulses provides a practical solution to enable not only GV/m surface fields but
also absolute timing stability, since the pulses are generated by the same
laser as the phototrigger. Here, we demonstrate an alloptical THz gun yielding
peak electron energies approaching 1 keV, accelerated by 300 MV/m THz fields in
a novel micronscale waveguide structure. We also achieve quasimonoenergetic,
subkeV bunches with 32 fC of charge, which can already be used for
timeresolved lowenergy electron diffraction. Such ultracompact, easy to
implement guns driven by intrinsically synchronized THz pulses that are pumped
by an amplified arm of the already present photoinjector laser provide a new
tool with potential to transform accelerator based science.

We refine the method towards extraction of subcycle transients in the 0.11
THz frequency (mmwavelength) range from optical rectification in lithium
niobate using tilted pulse fronts. Our scheme exploits previously unexplored
spatiotemporal shaping of the pump pulses, resulting in highly efficient and
near diffractionlimited subTHz transients reaching 0.2 GV/m electric field
strengths in freespace using only mJlevel optical pspulses. We address
experimentally and theoretically the means to producing aboveGV/m relativistic
subTHz transients with proper beam confinement and under moderate pumping
conditions, thereby bringing widespread access to strongfield and nonlinear
terahertz applications and devices.

Femtosecond electron bunches with keV energies and eV energy spread are
needed by condensed matter physicists to resolve state transitions in carbon
nanotubes, molecular structures, organic salts, and charge density wave
materials. These semirelativistic electron sources are not only of interest for
ultrafast electron diffraction, but also for electron energyloss spectroscopy
and as a seed for xray FELs. Thus far, the output energy spread (hence pulse
duration) of ultrafast electron guns has been limited by the achievable
electric field at the surface of the emitter, which is 10 MV/m for DC guns and
200 MV/m for RF guns. A singlecycle THz electron gun provides a unique
opportunity to not only achieve GV/m surface electric fields but also with
relatively low THz pulse energies, since a singlecycle transformlimited
waveform is the most efficient way to achieve intense electric fields. Here,
electron bunches of 50 fC from a flat copper photocathode are accelerated from
rest to tens of eV by a microjoule THz pulse with peak electric field of 72
MV/m at 1 kHz repetition rate. We show that scaling to the readilyavailable
GV/m THz field regime would translate to monoenergetic electron beams of ~100
keV.

The cost, size and availability of electron accelerators is dominated by the
achievable accelerating gradient. Conventional highbrightness radiofrequency
(RF) accelerating structures operate with 3050 MeV/m gradients. Electron
accelerators driven with optical or infrared sources have demonstrated
accelerating gradients orders of magnitude above that achievable with
conventional RF structures. However, laserdriven wakefield accelerators
require intense femtosecond sources and direct laserdriven accelerators and
suffer from low bunch charge, submicron tolerances and subfemtosecond timing
requirements due to the short wavelength of operation. Here, we demonstrate the
first linear acceleration of electrons with keV energy gain using
opticallygenerated terahertz (THz) pulses. THzdriven accelerating structures
enable highgradient electron or proton accelerators with simple accelerating
structures, high repetition rates and significant charge per bunch. Increasing
the operational frequency of accelerators into the THz band allows for greatly
increased accelerating gradients due to reduced complications with respect to
breakdown and pulsed heating. Electric fields in the GV/m range have been
achieved in the THz frequency band using all optical methods. With recent
advances in the generation of THz pulses via optical rectification of slightly
subpicosecond pulses, in particular improvements in conversion efficiency and
multicycle pulses, increasing accelerating gradients by two orders of
magnitude over conventional linear accelerators (LINACs) has become a
possibility. These ultracompact THz accelerators with extremely short electron
bunches hold great potential to have a transformative impact for free electron
lasers, future linear particle colliders, ultrafast electron diffraction,
xray science, and medical therapy with xrays and electron beams.

Compact laserdriven accelerators are versatile and powerful tools of
unarguable relevance on societal grounds for the diverse purposes of science,
health, security, and technology because they bring enormous practicality to
stateoftheart achievements of conventional radiofrequency accelerators.
Current benchmarking laserbased technologies rely on a medium to assist the
lightmatter interaction, which impose material limitations or strongly
inhomogeneous fields. The advent of few cycle ultraintense radially polarized
lasers has materialized an extensively studied novel accelerator that adopts
the simplest form of laser acceleration and is unique in requiring no medium to
achieve strong longitudinal energy transfer directly from laser to particle.
Here we present the first observation of direct longitudinal laser acceleration
of nonrelativistic electrons that undergo highlydirectional multiGeV/m
accelerating gradients. This demonstration opens a new frontier for direct
laserdriven particle acceleration capable of creating well collimated and
relativistic attosecond electron bunches and xray pulses.

We experimentally investigate the limits to 800 nmtoterahertz (THz) energy
conversion in lithium niobate at room temperature driven by amplified
Ti:Sapphire laser pulses with tiltedpulsefront. The influence of the pump
central wavelength, pulse duration, and fluence on THz generation is studied.
We achieved a high peak efficiency of 0.12% using transform limited 150 fs
pulses and observed saturation of the optical to THz conversion efficiency at a
fluence of 15 mJ/cm2. We experimentally identify two main limitations for the
scaling of opticaltoTHz conversion efficiencies: (i) the large spectral
broadening of the optical pump spectrum in combination with large angular
dispersion of the tiltedpulsefront and (ii) freecarrier absorption of THz
radiation due to multiphoton absorption of the 800 nm radiation.

We derive an asymptotically accurate formula for the beam waist of
ultrashort, tightlyfocused fundamental linearlypolarized, radiallypolarized,
and azimuthallypolarized modes in free space. We compute the exact beam waist
via numerical cubature to ascertain the accuracy with which our formula
approximates the exact beam waist over a broad range of parameters of practical
interest. Based on this, we describe a method of choosing parameters in the
model given the beam waist and pulse duration of a laser pulse.

We report on a multimJ 2.1$\mu$m OPCPA system operating at a 1kHz
repetition rate, pumped by a picosecond cryogenic Yb:YAG pump laser, and the
phasematched highflux highharmonic soft Xray generation.

We numerically investigate the acceleration and bunch compression
capabilities of 20 mJ, 0.6 THzcentered coherent terahertz pulses in optimized
metallic dielectricloaded cylindrical waveguides. In particular, we
theoretically demonstrate the acceleration of 1.6 pC and 16 pC electron bunches
from 1 MeV to 10 MeV over an interaction distance of 20mm, the compression of a
1.6 pC 1 MeV bunch from 100 fs to 2 fs (50 times compression) over an
interaction distance of about 18mm, and the compression of a 1.6 pC 10 MeV
bunch from 100 fs to 1.61 fs (62 times) over an interaction distance of 42 cm.
The obtained results show the promise of coherent THz pulses in realizing
compact electron acceleration and bunch compression schemes.

We propose a direct electron acceleration scheme that uses a twocolor pulsed
radiallypolarized laser beam. The twocolor scheme achieves electron
acceleration exceeding 90% of the theoretical energy gain limit, over twice of
what is possible with a onecolor pulsed beam of equal total energy and pulse
duration. The scheme succeeds by exploiting the Gouy phase shift to cause an
accelerationfavoring interference of fields only as the electron enters its
effectively final accelerating cycle. Optimization conditions and power scaling
characteristics are discussed.

We study the direct acceleration of a free electron in infinite vacuum along
the axis of a pulsed radiallypolarized laser beam. We find that net energy
transfer from laser pulse to electron is maximized with the tightest focusing.
We show that the net energy gain of an electron initially moving at a
relativistic velocity may exceed more than half the theoretical limit of energy
transfer, which is not possible with an initially stationary electron in the
parameter space studied. We determine and analyze the power scaling of maximum
energy gain, extending our study to include a relatively unexplored regime of
low powers and revealing that substantial acceleration is already possible
without the use of petawatt peakpower laser technology.

We hypothesize that a charged particle in unbounded vacuum can be
substantially accelerated by a force linear in the electric field of a
propagating electromagnetic wave only if the accelerating field is capable of
bringing the particle to a relativistic energy in its initial rest frame during
the interaction. We consequently derive a general formula for the acceleration
threshold of such schemes and support our conclusion with the results of
numerical simulations over a broad range of parameters for different kinds of
pulsed laser beams.

We experimentally and numerically study the atomic response and pulse
propagation effects of highorder harmonics generated in Xe, Kr, and Ar driven
by a 2.1\mu m infrared femtosecond light source. The light source is an
optical parametric chirpedpulse amplifier, and a modified strongfield
approximation and 3dimensional pulse propagation code are used for the
numerical simulations. The extended cutoff in the longwavelength driven
highharmonic generation has revealed the spectral shaping of highorder
harmonics due to the atomic structure (or photorecombination crosssection)
and the macroscopic effects, which are the main factors of determining the
conversion efficiency besides the driving wavelength. Using precise numerical
simulations to determine the macroscopic electron wavepacket, we are able to
extract the photorecombination crosssections from experimental highorder
harmonic spectra in the presence of macroscopic effects. We have experimentally
observed that the macroscopic effects shift the observed Cooper minimum of Kr
from 80 eV to 6070 eV and wash out the Cooper minimum of Ar. Measured
highharmonic conversion efficiencies per harmonic near the cutoff are ~10^{9}
for all three gases.

We present a study of direct laser driven electron acceleration and scaling
of attosecond bunch compression in unbound vacuum. Simple analytical
expressions and detailed threedimensional numerical calculations including
space charge reveal the conditions for compression to attosecond electron
sheets. Intermediate emittance minima suitable for brilliant xray generation
via coherent inverse Compton scattering (ICS) are predicted. We verify the
coherent emission properties of the resulting xray fields and demonstrate
feasability for realistic laser parameters.

We theoretically study the carrierenvelope phase dependent inversion
generated in a twolevel system by excitation with a fewcycle pulse. Based on
the invariance of the inversion under time reversal of the exciting field,
parameters are introduced to characterize the phase sensitivity of the induced
inversion. Linear and nonlinear phase effects are numerically studied for
rectangular and sincshaped pulses. Furthermore, analytical results are
obtained in the limits of weak fields as well as strong dephasing, and by
nearly degenerate perturbation theory for sinusoidal excitation. The results
show that the phase sensitive inversion in the ideal twolevel system is a
promising route for constructing carrierenvelope phase detectors.

Using the Kantorovitch method in combination with a Gaussian ansatz, we
derive the equations of motion for spatial, temporal and spatiotemporal optical
propagation in a dispersive Kerr medium with a general transverse and spectral
gain profile. By rewriting the variational equations as differential equations
for the temporal and spatial Gaussian q parameters, optical ABCD matrices for
the Kerr effect, a general transverse gain profile and nonparabolic spectral
gain filtering are obtained. Further effects can easily be taken into account
by adding the corresponding ABCD matrices. Applications include the temporal
pulse dynamics in gain fibers and the beam propagation or spatiotemporal pulse
evolution in bulk gain media. As an example, the steadystate spatiotemporal
Gaussian pulse dynamics in a Kerrlens modelocked laser resonator is studied.

We propose and analyze a new approach to generate a broadband astrocomb by
spectral broadening of a narrowband astrocomb inside a highly nonlinear
optical fiber. Numerical modeling shows that cascaded fourwavemixing
dramatically degrades the input comb's sidemode suppression and causes
sidemode amplitude asymmetry. These two detrimental effects can systematically
shift the centerofgravity of astrocomb spectral lines as measured by an
astrophysical spectrograph with resolution \approx100,000; and thus lead to
wavelength calibration inaccuracy and instability. Our simulations indicate
that this performance penalty, as a result of nonlinear spectral broadening,
can be compensated by using a filtering cavity configured for doublepass. As
an explicit example, we present a design based on an Ybfiber source comb (with
1 GHz repetition rate) that is filtered by doublepassing through a low finesse
cavity (finesse = 208), and subsequent spectrally broadened in a 2cm,
SF6glass photonic crystal fiber. Spanning more than 300 nm with 16 GHz line
spacing, the resulting astrocomb is predicted to provide 1 cm/s (~10 kHz)
radial velocity calibration accuracy for an astrophysical spectrograph. Such
extreme performance will be necessary for the search for and characterization
of Earthlike extrasolar planets, and in direct measurements of the change of
the rate of cosmological expansion.