Origin of the double-peak structure in longitudinal momentum distributions for single ionization of an He atom in strong laser field

The single ionization of an He atom by intense linearly polarized laser field in the tunneling regime is studied by S- matrix theory. When only the first term of the expansion of the S matrix is considered and time, spatial distribution, and fluctuation of the laser pulse are taken into account, the obtained momentum distribution in the polarization direction of laser field is consistent with the semiclassical calculation, which only considers tunneling and the interaction between the free electron and external field. When the second term, which includes the interaction between the core and the free electron, is considered, the momentum distribution shows a complex multipeak structure with the central minimum and the positions of some peaks are independent of the intensity in some intensity regime, which is consistent with the recent experimental result. Based on our analysis, we found that the structures observed in the momentum distribution of an He atom are attributed to the " soft" collision of the tunneled electron with the core.

Optimal Data Distributions in Machine Learning

In the first part of the thesis we explore three fundamental questions that arise naturally when we conceive a machine learning scenario where the training and test distributions can differ. Contrary to conventional wisdom, we show that in fact mismatched training and test distribution can yield better out-of-sample performance. This optimal performance can be obtained by training with the dual distribution. This optimal training distribution depends on the test distribution set by the problem, but not on the target function that we want to learn. We show how to obtain this distribution in both discrete and continuous input spaces, as well as how to approximate it in a practical scenario. Benefits of using this distribution are exemplified in both synthetic and real data sets.

In order to apply the dual distribution in the supervised learning scenario where the training data set is fixed, it is necessary to use weights to make the sample appear as if it came from the dual distribution. We explore the negative effect that weighting a sample can have. The theoretical decomposition of the use of weights regarding its effect on the out-of-sample error is easy to understand but not actionable in practice, as the quantities involved cannot be computed. Hence, we propose the Targeted Weighting algorithm that determines if, for a given set of weights, the out-of-sample performance will improve or not in a practical setting. This is necessary as the setting assumes there are no labeled points distributed according to the test distribution, only unlabeled samples.

Finally, we propose a new class of matching algorithms that can be used to match the training set to a desired distribution, such as the dual distribution (or the test distribution). These algorithms can be applied to very large datasets, and we show how they lead to improved performance in a large real dataset such as the Netflix dataset. Their computational complexity is the main reason for their advantage over previous algorithms proposed in the covariate shift literature.

In the second part of the thesis we apply Machine Learning to the problem of behavior recognition. We develop a specific behavior classifier to study fly aggression, and we develop a system that allows analyzing behavior in videos of animals, with minimal supervision. The system, which we call CUBA (Caltech Unsupervised Behavior Analysis), allows detecting movemes, actions, and stories from time series describing the position of animals in videos. The method summarizes the data, as well as it provides biologists with a mathematical tool to test new hypotheses. Other benefits of CUBA include finding classifiers for specific behaviors without the need for annotation, as well as providing means to discriminate groups of animals, for example, according to their genetic line.

Stationary absolute distributions for chains of infinite order

Let {Ƶ_n}^∞_{n = -∞} be a stochastic process with state space S₁ = {0, 1, …, D – 1}. Such a process is called a chain of infinite order. The transitions of the chain are described by the functions

Q_i(i⁽⁰⁾) = Ƥ(Ƶ_n = i | Ƶ_{n - 1} = i ⁽⁰⁾₁, Ƶ_{n - 2} = i ⁽⁰⁾₂, …) (i ɛ S₁), where i⁽⁰⁾ = (i⁽⁰⁾₁, i⁽⁰⁾₂, …) ranges over infinite sequences from S₁. If i⁽ⁿ⁾ = (i⁽ⁿ⁾₁, i⁽ⁿ⁾₂, …) for n = 1, 2,…, then i⁽ⁿ⁾ → i⁽⁰⁾ means that for each k, i⁽ⁿ⁾_k = i⁽⁰⁾_k for all n sufficiently large.

Given functions Q_i(i⁽⁰⁾) such that

(i) 0 ≤ Q_i(i⁽⁰) ≤ ξ ˂ 1

(ii)D – 1/Ʃ/i = 0 Q_i(i⁽⁰⁾) Ξ 1

(iii) Q_i(i⁽ⁿ⁾) → Q_i(i⁽⁰⁾) whenever i⁽ⁿ⁾ → i⁽⁰⁾,

we prove the existence of a stationary chain of infinite order {Ƶ_n} whose transitions are given by

Ƥ (Ƶ_n = i | Ƶ_{n - 1}, Ƶ_{n - 2}, …) = Q_i(Ƶ_{n - 1}, Ƶ_{n - 2}, …)

With probability 1. The method also yields stationary chains {Ƶ_n} for which (iii) does not hold but whose transition probabilities are, in a sense, “locally Markovian.” These and similar results extend a paper by T.E. Harris [Pac. J. Math., 5 (1955), 707-724].

Included is a new proof of the existence and uniqueness of a stationary absolute distribution for an Nth order Markov chain in which all transitions are possible. This proof allows us to achieve our main results without the use of limit theorem techniques.

Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

Transverse intensity distributions of a broadband laser modulated by a hard-edged aperture

By means of the Huygens-Fresnel diffraction integral, the field representation of a laser beam modulated by a hard-edged aperture is derived. The near-field and far-field transverse intensity distributions of the beams with different bandwidths are analyzed by using the representation. The numerical calculation results indicate that the amplitudes and numbers of the intensity spikes decrease with increasing bandwidth, and beam smoothing is achieved when the bandwidth takes a certain value in the near field. In the far field, the radius of the transverse intensity distribution decreases as the bandwidth increases, and the physical explanation of this fact is also given. (c) 2005 Optical Society of America.

Influence of the chirp on the intensity distributions of an apertured pulse

Starting from the Huygens-Fresnel diffraction integral and the Fourier transform, the propagation expression of a chirped pulse passing through a hard-edged aperture is derived. Using the obtained expression, the intensity distributions of the pulse with different chirp in the near and far fields are analyzed in detail. Due to the modulation of the aperture, many intensity peaks emerge in the intensity distributions of the chirped pulse in the near field. However, the amplitudes of the intensity peaks decrease on increasing the chirp, which results in the smoothing effect in the intensity distributions. The beam smoothing brought by increasing the chirp is explained physically. Also, it is found that the radius of the intensity distribution of the chirped pulse decreases when the chirp increases in the far field. (c) 2005 Elsevier GmbH. All rights reserved.

板条热容激光器的二维热特性

固体热容激光器（SSHCL）作为高功率固体激光器的一个重要发展方向,引起人们广泛关注。数值模拟激光介质板条在热容方式下工作的温度和应力分布是了解该类激光器工作特性的一种有效手段,采用平面应力近似法导出了半导体激光器抽运热容激光介质板的二维温度和应力分布公式,同时也对二维抽运光吸收密度、介质板温度分布和折射率变化进行了分析与讨论。数值计算的结果表明二维效应的温度分布和应力分布要比一维效应给出的分布更均匀。

Negative refraction via domain wall resonances in a homogeneous mixture of ferro- and nonmagnetic substances

We deliver the general conditions on the synthetic proportions for a homogeneous mixture of ferro- and nonmagnetic substances to become left-handed. As an alternative for left-handed metamaterials, we consider mixing ferromagnetic materials with nonmagnetic microscopic particles. In the mixture, the ferromagnetic material provides the needed permeability via domain wall resonances at high frequencies, whereas the nonmagnetic material gives the required permittivity. Using the effective medium theory, we have found that when the concentration of the nonmagnetic particles falls into a certain range, the refractive index of the mixture is negative, n < 0, which includes the double negative ( epsilon < 0 and mu < 0) and other cases ( e. g. epsilon < 0 and mu > 0). We finally give the requirements on the microscopic material properties for the ferromagnetic materials to reach the domain wall resonances at high frequencies.