Cross-propagating Beam-Deflection measurements of third-order nonlinear optical susceptibility

Similar documents
Periodic Nonlinear Refractive Index of Carbon Disulfide Vapors

Nonlinear Optics (NLO)

Terahertz Kerr effect

Doctor of Philosophy

gives rise to multitude of four-wave-mixing phenomena which are of great

Temporal and polarization dependence of the nonlinear optical response of solvents

BY TEMPORALLY TWO-DIMENSIONAL

requency generation spectroscopy Rahul N

CHAPTER 7 SUMMARY OF THE PRESENT WORK AND SUGGESTIONS FOR FUTURE WORK

White-light continuum Z-scan technique for nonlinear materials characterization

Optics, Light and Lasers

Nanocomposite photonic crystal devices

Optical and Photonic Glasses. Lecture 37. Non-Linear Optical Glasses I - Fundamentals. Professor Rui Almeida

Optical Spectroscopy of Advanced Materials

Advanced Vitreous State The Physical Properties of Glass

LIST OF TOPICS BASIC LASER PHYSICS. Preface xiii Units and Notation xv List of Symbols xvii

Supplemental material for Bound electron nonlinearity beyond the ionization threshold

Supplementary Figures

As a partial differential equation, the Helmholtz equation does not lend itself easily to analytical

Liquid Crystals IAM-CHOON 1(1100 .,4 WILEY 2007 WILEY-INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION. 'i; Second Edition. n z

Time resolved optical spectroscopy methods for organic photovoltaics. Enrico Da Como. Department of Physics, University of Bath

Multidimensional femtosecond coherence spectroscopy for study of the carrier dynamics in photonics materials

Temporal, spectral, and polarization dependence of the nonlinear optical response of carbon disulfide

Linear and Nonlinear Optical Properties of Acridine Dye Doped PMMA Polymer

Lasers and Electro-optics

Elements of Quantum Optics

Nonlinear Optics. Second Editio n. Robert W. Boyd

1 Fundamentals of laser energy absorption

Linear pulse propagation

Time-resolved pump-probe system based on a nonlinear imaging technique with phase object

Multi-cycle THz pulse generation in poled lithium niobate crystals

5.74 Introductory Quantum Mechanics II

Molecular alignment, wavepacket interference and Isotope separation

Mansoor Sheik-Bahae. Class meeting times: Mondays, Wednesdays 17:30-18:45 am; Physics and Astronomy, Room 184

Theoretical simulation of nonlinear spectroscopy in the liquid phase La Cour Jansen, Thomas

Supporting Information

Magnetostatic modulation of nonlinear refractive index and absorption in quantum wires

Chap. 4. Electromagnetic Propagation in Anisotropic Media

Survey on Laser Spectroscopic Techniques for Condensed Matter

Spectral Resolution. Spectral resolution is a measure of the ability to separate nearby features in wavelength space.

Final Report for AOARD grant FA Measurement of the third-order nonlinear susceptibility of graphene and its derivatives

F. Elohim Becerra Chavez

Photoelectron Spectroscopy using High Order Harmonic Generation

SURFACE-DISPLACEMENT IMAGING USING OPTICAL BEAM DEFLECTION. S.E. McBride* and G.C. Wetsel, Jr.* Southern Methodist University Dallas, Texas 75275

Highly Efficient and Anomalous Charge Transfer in van der Waals Trilayer Semiconductors

SUPPLEMENTARY INFORMATION

NONLINEAR OPTICS. Ch. 1 INTRODUCTION TO NONLINEAR OPTICS

Introduction to Nonlinear Optics

Implementation and evaluation of data analysis strategies for time-resolved optical spectroscopy

The generation of terahertz frequency radiation by optical rectification

No. 9 Experimental study on the chirped structure of the construct the early time spectra. [14;15] The prevailing account of the chirped struct

Nonlinear Optics (NLO)

Observation of spectral enhancement in a soliton fiber laser with fiber Bragg grating

SUPPLEMENTARY INFORMATION

3.5 Cavities Cavity modes and ABCD-matrix analysis 206 CHAPTER 3. ULTRASHORT SOURCES I - FUNDAMENTALS

Second Harmonic Generation in Solid-State Materials

Single Emitter Detection with Fluorescence and Extinction Spectroscopy

Lukas Gallmann. ETH Zurich, Physics Department, Switzerland Chapter 4b: χ (2) -nonlinearities with ultrashort pulses.

CONTENTS. 2 CLASSICAL DESCRIPTION 2.1 The resonance phenomenon 2.2 The vector picture for pulse EPR experiments 2.3 Relaxation and the Bloch equations

Model Answer (Paper code: AR-7112) M. Sc. (Physics) IV Semester Paper I: Laser Physics and Spectroscopy

Nanoscale Energy Conversion and Information Processing Devices - NanoNice - Photoacoustic response in mesoscopic systems

Electromagnetic fields and waves

Probing and Driving Molecular Dynamics with Femtosecond Pulses

F.G. Major. The Quantum Beat. The Physical Principles of Atomic Clocks. With 230 Illustrations. Springer

(Introduction) Linear Optics and Nonlinear Optics

Electronic Supplementary Information

Nonlinear Effects in Optical Fiber. Dr. Mohammad Faisal Assistant Professor Dept. of EEE, BUET

Probing the Propagation Dynamics in a Femtosecond Laser Filament

INTRODUCTION TO NONLINEAR OPTICAL EFFECTS IN MOLECULES AND POLYMERS

A tutorial on meta-materials and THz technology

Observation of the waveform of accumulated photon echoes in a dye-doped polymer film by use of an interferometer

Control of dispersion effects for resonant ultrashort pulses M. A. Bouchene, J. C. Delagnes

Tailorable stimulated Brillouin scattering in nanoscale silicon waveguides.

Time Resolved Faraday Rotation Measurements of Spin Polarized Currents in Quantum Wells

Laserphysik. Prof. Yong Lei & Dr. Yang Xu. Fachgebiet Angewandte Nanophysik, Institut für Physik

Richard Miles and Arthur Dogariu. Mechanical and Aerospace Engineering Princeton University, Princeton, NJ 08540, USA

F. Elohim Becerra Chavez

Picosecond z-scan measurements of the two-photon absorption in beta-carotene solution over the nm wavelength range

PRINCIPLES OF PHYSICAL OPTICS

1 Mathematical description of ultrashort laser pulses

Nonlinear effects and pulse propagation in PCFs

χ (3) Microscopic Techniques

ATOMIC AND LASER SPECTROSCOPY

Saturation Absorption Spectroscopy of Rubidium Atom

PROCEEDINGS OF SPIE. Imaging carrier dynamics on the surface of the N-type silicon

Let us consider a typical Michelson interferometer, where a broadband source is used for illumination (Fig. 1a).

Ultrafast Mechanisms of Nonlinear Refraction and Two-photon Photochromism

Supporting Online Material for

Supplementary Figure 1 Schematics of an optical pulse in a nonlinear medium. A Gaussian optical pulse propagates along z-axis in a nonlinear medium

Self-Phase Modulation in Optical Fiber Communications: Good or Bad?

Optics, Optoelectronics and Photonics

Slow Light in Crystals

Grading. Class attendance: (1 point/class) x 9 classes = 9 points maximum Homework: (10 points/hw) x 3 HW = 30 points maximum

Digital Holographic Measurement of Nanometric Optical Excitation on Soft Matter by Optical Pressure and Photothermal Interactions

Two-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO2 gas

Optical Properties of Thin Semiconductor Films

Nonlinear optical properties of nanocomposite materials

Laser Physics OXFORD UNIVERSITY PRESS SIMON HOOKER COLIN WEBB. and. Department of Physics, University of Oxford

Carrier dynamics of rubrene single-crystals revealed by transient broadband terahertz

Slowing Down the Speed of Light Applications of "Slow" and "Fast" Light

Transcription:

Cross-propagating Beam-Deflection measurements of third-order nonlinear optical susceptibility Sepehr Benis, David J. Hagan, Eric W. Van Stryland CREOL, The College of Optics and Photonics, P.O. Box 162700, University of Central Florida, Orlando, FL 32816 USA ABSTRACT We extend our recently developed Beam-Deflection (BD) method which we used to determine the sign, magnitude, and ultrafast dynamics of bound electronic and nuclear nonlinear optical responses, to measure elements of the third-order nonlinear optical susceptibility tensor with longitudinal field components. In these measurements, in contrast to the conventional BD technique where excitation and probe beams are nearly collinear, the interaction geometry involves orthogonal propagation of excitation and probe beams. This cross-propagating BD method enables probing with the electric field polarized parallel to the wavevector of the excitation beam as well as conventional parallel and perpendicular polarizations as in most other experimental methods. This technique may be of use in detecting possible magneto-electric contributions to the nonlinear susceptibility. The ratio between parallel and perpendicular bound-electronic responses in this method is shown to be larger than the factor of 3 predicted by theory and measured in the conventional configuration. Keywords: Ultrafast Nonlinear Optics; Nonlinear Optical Materials; Ultrafast Phenomena; Kerr Effect 1. INTRODUCTION Nonlinear refraction (NLR) and absorption (NLA) cause changes in the optical properties of a material as the irradiance of light increases. These phenomena open a variety of applications such as generation of ultrashort pulses, all-optical signal processing, and ultrafast switching [1]. The ever-increasing role of nonlinear optics has created a pressing demand for the characterization of the nonlinear optical susceptibility of various materials. The nonlinear susceptibility tensor is controlled by symmetry restrictions which are imposed by the spatial symmetry of the nonlinear medium. These limitations may reduce the number of independent coefficients or tensor elements for a given nonlinear susceptibility tensor. The lowestorder nonlinear response with non-vanishing terms is the third-order nonlinearity. Various experimental techniques are devoted to measuring the nonlinear optical (NLO) response. For example, Z-scan [2] measures the magnitude and sign of the NLR and NLA. Also, different time-resolved and polarization dependent techniques have been reported for measuring the third-order nonlinear properties of materials such as the optical Kerr effect (OKE) [3-6], degenerate four-wave mixing (DWFM) [7], as well as nonlinear interferometry [8]. However, the polarization combinations are often limited because of the interaction geometry. To fully exploit the third-order nonlinear optical properties of materials, we modified our Beam- Deflection (BD) technique [9] to study the elements of the third-order susceptibility tensor with field components along the propagation direction. BD is a recently developed technique for measuring the sign, magnitude, and ultrafast dynamics of the third-order NLO response of a material [9]. Nonlinear refraction occurs at any wavelength in different phases of matter, e.g. liquids, gases, and solids, and may exhibit non-instantaneous responses [10-12]. The NLO response of materials includes both bound-electronic and nuclear contributions. The great sensitivity of the BD method (demonstrated sensitivity to nonlinearly-induced phase distortion of /20,000) allows such measurements [9]. Here we present an extension of the BD technique which involves a crossed-beam interaction geometry for excitation and probe beams. Here we use femtosecond laser pulses. This cross-propagating beam geometry was first used for characterizing the change in the refractive index induced by absorption due to thermal nonlinearities in the case of small thermal diffusion lengths compared to the excitation beam radius [13, 14]. In this geometry, the optical interaction length is determined by a combination of the beam sizes and shapes along with the pulsewidths. The cross-propagating geometry allows measurements of the nonlinear susceptibility tensor elements with longitudinal field components in addition to those available with the standard geometries. The control of polarization combinations helps to measure the ratio between different third-order susceptibility tensor elements. This may be of particular interest for the detection of the recently proposed magneto-electric (ME) Nonlinear Frequency Generation and Conversion: Materials and Devices XVI, edited by Konstantin L. Vodopyanov, Kenneth L. Schepler, Proc. of SPIE Vol. 10088, 100880N 2017 SPIE CCC code: 0277-786X/17/$18 doi: 10.1117/12.2251451 Proc. of SPIE Vol. 10088 100880N-1

susceptibility, where a nonlinear magnetic dipole is induced parallel to the magnetic optical field direction [15-17]. As we shall show, the crossed-beam method allows us to separately look at cases where the magnetic vectors in pump and probe are parallel, but the electric fields are perpendicular, and vice versa, perhaps allowing identification of magneto-electric contributions to the nonlinear susceptibility. In this work, we present polarization resolved cross-propagating BD measurements of a few selected organic solvents (CCl 4, SnCl 4, and Anisaldehyde) to study the third-order susceptibility tensor elements and responses of these materials. 2. CROSS-PROPAGATING BEAM-DEFLECTION TECHNIQUE Beam-Deflection is a time-resolved and polarization-dependent technique for measuring the magnitude, sign, and dynamics of NLR and NLA. The principle of the BD method is similar to transient absorption spectroscopy with some modifications for measuring NLR. This method can measure the dynamics of NLR and NLA as opposed to techniques such as Z-scan [2] which only measures the magnitude and sign of the temporally and spatially averaged nonlinearities. Different polarization combinations between excitation and probe beams reveal information about both the instantaneous bound-electronic response and non-instantaneous nuclear motions. Three useful polarization combinations are copolarized (parallel), cross-polarized (perpendicular), and magic angle (54.7 ). At the magic angle, nuclear reorientation does not contribute to the detected signal, thus, in many situations, only the collision-induced and instantaneous bound electron response is detected. For highly symmetric molecules such as carbon tetrachloride, the reorientational response from nuclear reorientation vanishes [18]. When a strong excitation pulse interacts with a Kerr medium, it creates an index gradient that follows the Gaussian spatial profile of the irradiance of the beam, I e (r; t), in both transverse and longitudinal planes. The induced phase gradient appears as a refractive index gradient in its wings and acts as a thin prism. When the probe beam is focused to a spot size ~ 3-5 smaller than that of the excitation beam and displaced to the wings of the induced phase gradient, it experiences a linear index gradient. This deflects the probe beam by a small angle which is measured using a position sensitive segmented detector in the far-field. The nearly collinearly interaction scheme, requires the polarization of the excitation and probe beams to be in the transverse plane and, thus, it limits the polarization combinations to a single plane. By modifying this method to the cross-propagating geometry, we can measure the induced index change for polarization combinations which are not permitted by the collinear geometry, allowing the detection of suspected magneto-electric [15-17] susceptibility contributions to the nonlinear refractive index along the wavevector of the excitation beam. In the present configuration, the experimental geometry entails a complex interaction scheme between excitation and probe beams and restricts the interaction length to the temporal and spatial dimensions of the two pulses, see Fig. 1(a). In other words, the detected signal is not directly proportional to the thickness of the sample and the effective interaction length is determined by the beam pulsewidths (τ e, and τ p ) as well as spot sizes (w e, and w p ) of both the excitation and probe beams, respectively. This can also help to extract the pulsewidth and spot size of the excitation beam from the crosscorrelation signal. The effective interaction length is less than the spatial extent of the pulse; therefore, the spatial overlap of two pulses is relatively small within the interval of temporal overlap, and the deflection signal is significantly smaller than in the collinear case. To eliminate the need for tight focusing along the deflection axis, a cylindrical lens is used to tightly focus the excitation beam in one axis. When using two different wavelengths for the excitation and probe beams, the effects of group velocity mismatch is reduced due to the smaller interaction length comparing to the collinear case. Also, in the case of measurements on liquids, the cuvette walls do not contribute to the detected signal because the excitation and probe beams are interacting inside the liquid sample and do not spatially overlap on the cuvette walls. Since excitation and probe beams have perpendicular propagation vectors, it is possible to vary the polarization of the two beams to measure configurations different from the collinear case. In addition to the two standard parallel and perpendicular cases, another two perpendicular polarization configurations as well as magic angle configurations are measurable, see Fig. 1(b, c). This method enables the possibility of measuring two combinations where the magnetic optical fields of excitation and probe beams are parallel (χ zyyz ) or perpendicular (χ xyyx ) but the electric fields are perpendicular. This may help to measure the effect of nonlinear magnetic dipole moments which are along the magnetic field of the excitation beam [15-17]. The induced irradiance-dependent refractive index in different configurations, is given by: n j (I) n 0 + 3 χ 4n jiij E i,e 2 (i, j = x, y, z), (1) 0 where n 0 is the linear index of refraction, χ jiij is the corresponding tensor element, and E i,e is the amplitude of the excitation field polarized along i. The deflection angle, θ, is also proportional to the change in the refractive index and is given by θ [ n(r, t)]ds, (2) Proc. of SPIE Vol. 10088 100880N-2

where ds is an infinitesimal displacement along the beam path [9]. a) Probe Segmented Detector b) xxxxx xxyyx Ée c) y Excitation Plane e = 54.7' / B Probe Plane e = 54.7' Fig. 1. (a) Schematic diagram of cross-propagating Beam-Deflection experiment and the interaction of excitation and probe beams. The probe is spatially displaced to the wings of the Gaussian profile of the excitation beam. (b) Different polarization configurations for parallel, perpendicular, and (c) magic angle between excitation (E e, B e ) and probe (E p, B p ). Case (b) offers opportunities to study the cases where electric fields of the excitation and probe beams are perpendicular but the magnetic fields are parallel, or vice-versa, or where both are perpendicular. A regeneratively-amplified femtosecond laser (Clark-MXR) with 775 nm output pulse, 1 khz repetition rate, 1 mj pulse energy, and ~150 fs (FWHM) pulsewidth is used. A portion of the beam is focused into a 1 cm quartz cuvette filled with water to generate a white-light continuum, and then spectrally filtered by a narrow band-pass filter ( λ~10 nm) at 650 nm to be used as the probe. A combination of half-wave plate, and polarizer is used to control the energy, and another half-wave plate is used to rotate the polarization of the probe (similarly for the excitation). The probe beam is focused in the sample and propagated to the far field and collected by a quad-segmented Silicon photodiode (OSI QD5-0-SD). To increase the signal-to-noise ratio we modulate the 1 khz repetition rate excitation beam at 285 Hz for lock-in detection. The excitation beam is also filtered at 775 nm by using a narrow band-pass filter (~ 10 nm) and temporally delayed to measure the dynamics. Ee L y Be Proc. of SPIE Vol. 10088 100880N-3

3. RESULTS AND DISCUSSION Cross-propagating Beam-Deflection measurements are performed on a few symmetric molecules such as carbon tetrachloride (CCl 4) and tin tetrachloride (SnCl 4). In the non-resonant limit, the magnitude of the bound-electronic nonlinearity in isotropic molecules is predicted to be a factor of 3 smaller for crossed polarizations than for parallel polarizations. This theoretical prediction for isotropic tensor symmetry was verified by collinear BD geometry measurements in previous works [9-12]. Here, Gaussian fits to the results for carbon tetrachloride and tin tetrachloride are used to find the ratio between parallel and perpendicular polarization contributions, as shown in Fig. 2(a, b). It is intriguing to see a ratio of 4 compared to the collinear geometry, where our experiments and theory show that the ratio is 3 [9]. Measurements on isotropic liquids have not shown any difference between the different off-diagonal cases, (χ zyyz ) or (χ xyyx ) which indicates that the magnetic optical field direction does not seem to affect the result. This technique is also used in measuring dynamics of NLR of meta-anisaldehyde and para- Anisaldehyde organic liquids for parallel, perpendicular, and magic angle configurations as shown in Fig. 2(c, d). In this case, the non-instantaneous contributions to the BD signal are a result of nuclear rotations [9-12]. The polarization dependence of the non-instantaneous response can be observed at longer delay, where for parallel polarization the signal is positive, and for perpendicular polarization it is negative. At the magic angle the nuclear rotations are completely eliminated and the signal is only due to mechanisms with isotropic symmetry which we attribute to the instantaneous bound-electronic response and collision-induced mechanisms [9-12]. 3.0 2.5 2.0 --0- Parallel (xxxx) -0-- Perpendicular (zxxz) Perpendicular (xyyx) -0- Perpendicular (zyyz) Gaussian Fit 3.0 2.5 2.0 (b) SnCl4-0-Parallel (xxxx) -0- P e rpen dic ula r (zxxz) Perpendicular (x)yx) - 0-- Perpendicular (zyyz) - Gaussian Fit w 1.5 u7 1.0 0.5 0.5 0.0-1.5-1.0-0.5 0.0 0.5 1.0 1.5 0.0 0 8 6 _(C) m- Anisaldehyde Parallel (xxxx) Perpendicular (zxxz) -- Perpendicular (xyyx) - Perpendicular (zy)z) Magic Angle 8 6 (d) p- Anisaldehyde Parallel (xxxx) -Perpendicular(zxxz) - -Perpendicular(xyyx) Perpendicular (zyyz) Magic Angle w4 u7 2 0-1 0 1 2 3 4 5-1 0 1 2 3 4 5 Fig. 2. Beam-Deflection measurements of isotropic organic molecules in liquid form in different polarization configurations, for (a) carbon tetrachloride (CCl4), (b) Tin tetrachloride (SnCl4), (c) 3-Methoxybenzaldehyde (m-anisaldehyde), and (d) 4- Methoxybenzaldehyde (p-anisaldehyde). Proc. of SPIE Vol. 10088 100880N-4

4. CONCLUSION We extended our recently developed Beam-Deflection method, which allowed bound-electronic and each nuclear contribution to be explicitly separated, to a geometry where the excitation and probe propagate at 90 0 with respect to each other. Both geometries time-resolve the nonlinear response, allowing various elements of the bound-electronic and nuclear third-order nonlinear optical susceptibility tensor to be measured. In addition to being able to measure the usual tensor components available in the collinear propagation geometry, this cross-propagating BD method is able to measure the deflection signal with probe beam polarized along the wavevector of the excitation beam. This combination of geometries allows all tensor elements to contribute to the NLR and NLA depending on polarization combinations, and will, therefore, help in separately determining their contributions to the overall nonlinearity. The possibility of measuring tensor elements with different orientations of magnetic fields with respect to electric field, helps to investigate the possible contribution of magnetoelectric effects in the third-order nonlinear optical susceptibility which is under further study. The ratio of 4 between parallel and perpendicular polarization rather than the predicted 3 as measured in the collinear geometry, is not yet understood and is the subject of ongoing investigation. 5. ACKNOWLEDGEMENT This work was supported by the Air Force Office of Scientific Research (AFOSR) (FA9550-14-1-0040) MURI Center for Dynamic Magneto-Optics, and Army Research Lab (ARL) (W911NF-15-2-0090). REFERENCES [1] R. W. Boyd, [Nonlinear optics] Academic press, (2003). [2] M. Sheik-Bahae, A. A. Said, T.-H. Wei et al., Sensitive measurement of optical nonlinearities using a single beam, IEEE Journal of Quantum Electronics, 26(4), 760-769 (1990). [3] P. P. Ho, and R. R. Alfano, Optical Kerr effect in liquids, Physical Review A, 20(5), 2170-2187 (1979). [4] D. McMorrow, W. T. Lotshaw, and G. A. Kenney-Wallace, Femtosecond optical Kerr studies on the origin of the nonlinear responses in simple liquids, IEEE Journal of Quantum Electronics, 24(2), 443-454 (1988). [5] D. McMorrow, N. Thantu, V. Kleiman et al., Analysis of Intermolecular Coordinate Contributions to Third- Order Ultrafast Spectroscopy of Liquids in the Harmonic Oscillator Limit, The Journal of Physical Chemistry A, 105(34), 7960-7972 (2001). [6] Q. Zhong, and J. T. Fourkas, Optical Kerr Effect Spectroscopy of Simple Liquids, The Journal of Physical Chemistry B, 112(49), 15529-15539 (2008). [7] T. Steffen, J. T. Fourkas, and K. Duppen, Time resolved four and six wave mixing in liquids. I. Theory, The Journal of Chemical Physics, 105(17), 7364-7382 (1996). [8] Yasuhiro Sato, Ryuji Morita, and Mikio Yamashita, Study on Ultrafast Dynamic Behaviors of Different Nonlinear Refractive Index Components in CS 2 Using a Femtosecond Interferometer, Japanese Journal of Applied Physics, 36(4R), 2109 (1997). [9] M. R. Ferdinandus, H. Hu, M. Reichert et al., Beam deflection measurement of time and polarization resolved ultrafast nonlinear refraction, Optics Letters, 38(18), 3518-3521 (2013). [10] M. Reichert, H. Hu, M. R. Ferdinandus et al., Temporal, spectral, and polarization dependence of the nonlinear optical response of carbon disulfide, Optica, 1(6), 436-445 (2014). [11] M. Reichert, P. Zhao, J. M. Reed et al., Beam deflection measurement of bound-electronic and rotational nonlinear refraction in molecular gases, Optics Express, 23(17), 22224-22237 (2015). [12] P. Zhao, M. Reichert, D. J. Hagan et al., Dispersion of nondegenerate nonlinear refraction in semiconductors, Optics Express, 24(22), 24907-24920 (2016). [13] A. C. Boccara, W. Jackson, N. M. Amer et al., Sensitive photothermal deflection technique for measuring absorption in optically thin media, Optics Letters, 5(9), 377-379 (1980). [14] W. B. Jackson, N. M. Amer, A. C. Boccara et al., Photothermal deflection spectroscopy and detection, Applied Optics, 20(8), 1333-1344 (1981). [15] A. A. Fisher, E. F. Cloos, W. M. Fisher et al., Dynamic symmetry-breaking in a simple quantum model of magneto-electric rectification, optical magnetization, and harmonic generation, Optics Express, 22, 2910-2924 (2014). Proc. of SPIE Vol. 10088 100880N-5

[16] W. Fisher, and S. Rand, Light-induced dynamics in the Lorentz oscillator model with magnetic forces, Physical Review A, 82(1), 013802 (2010). [17] S. C. Rand, W. M. Fisher, and S. L. Oliveira, Optically induced magnetization in homogeneous, undoped dielectric media, Journal of the Optical Society of America B, 25(7), 1106-1117 (2008). [18] P. Zhao, M. Reichert, T. R. Ensley et al., "Nonlinear refraction dynamics of solvents and gases." 9731, 97310F- 97310F-8. Proc. of SPIE Vol. 10088 100880N-6