Nonlinear parameter retrieval from three- and four-wave mixing in metamaterials

We present a generalized nonlinear susceptibility retrieval method for metamaterials based on transfer matrices and valid in the nondepleted pump approximation. We construct a general formalism to describe the transfer matrix method for nonlinear media and apply it to the processes of three- and four-wave mixing. The accuracy of this approach is verified via finite element simulations. The method is then reversed to give a set of equations for retrieving the nonlinear susceptibility. Finally, we apply the proposed retrieval operation to a three-wave mixing transmission experiment performed on a varactor loaded split ring resonator metamaterial sample and find quantitative agreement with an analytical effective medium theory model. © 2010 The American Physical Society.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Millimeter-wave compressive holography.

We describe an active millimeter-wave holographic imaging system that uses compressive measurements for three-dimensional (3D) tomographic object estimation. Our system records a two-dimensional (2D) digitized Gabor hologram by translating a single pixel incoherent receiver. Two approaches for compressive measurement are undertaken: nonlinear inversion of a 2D Gabor hologram for 3D object estimation and nonlinear inversion of a randomly subsampled Gabor hologram for 3D object estimation. The object estimation algorithm minimizes a convex quadratic problem using total variation (TV) regularization for 3D object estimation. We compare object reconstructions using linear backpropagation and TV minimization, and we present simulated and experimental reconstructions from both compressive measurement strategies. In contrast with backpropagation, which estimates the 3D electromagnetic field, TV minimization estimates the 3D object that produces the field. Despite undersampling, range resolution is consistent with the extent of the 3D object band volume.

Effect of lithotripter focal width on stone comminution in shock wave lithotripsy.

Using a reflector insert, the original HM-3 lithotripter field at 20 kV was altered significantly with the peak positive pressure (p(+)) in the focal plane increased from 49 to 87 MPa while the -6 dB focal width decreased concomitantly from 11 to 4 mm. Using the original reflector, p(+) of 33 MPa with a -6 dB focal width of 18 mm were measured in a pre-focal plane 15-mm proximal to the lithotripter focus. However, the acoustic pulse energy delivered to a 28-mm diameter area around the lithotripter axis was comparable ( approximately 120 mJ). For all three exposure conditions, similar stone comminution ( approximately 70%) was produced in a mesh holder of 15 mm after 250 shocks. In contrast, stone comminution produced by the modified reflector either in a 15-mm finger cot (45%) or in a 30-mm membrane holder (14%) was significantly reduced from the corresponding values (56% and 26%) produced by the original reflector (no statistically significant differences were observed between the focal and pre-focal planes). These observations suggest that a low-pressure/broad focal width lithotripter field will produce better stone comminution than its counterpart with high-pressure/narrow focal width under clinically relevant in vitro comminution conditions.

Far-field analysis of axially symmetric three-dimensional directional cloaks.

Axisymmetric radiating and scattering structures whose rotational invariance is broken by non-axisymmetric excitations present an important class of problems in electromagnetics. For such problems, a cylindrical wave decomposition formalism can be used to efficiently obtain numerical solutions to the full-wave frequency-domain problem. Often, the far-field, or Fraunhofer region is of particular interest in scattering cross-section and radiation pattern calculations; yet, it is usually impractical to compute full-wave solutions for this region. Here, we propose a generalization of the Stratton-Chu far-field integral adapted for 2.5D formalism. The integration over a closed, axially symmetric surface is analytically reduced to a line integral on a meridional plane. We benchmark this computational technique by comparing it with analytical Mie solutions for a plasmonic nanoparticle, and apply it to the design of a three-dimensional polarization-insensitive cloak.

Improving the lens design and performance of a contemporary electromagnetic shock wave lithotripter.

The efficiency of shock wave lithotripsy (SWL), a noninvasive first-line therapy for millions of nephrolithiasis patients, has not improved substantially in the past two decades, especially in regard to stone clearance. Here, we report a new acoustic lens design for a contemporary electromagnetic (EM) shock wave lithotripter, based on recently acquired knowledge of the key lithotripter field characteristics that correlate with efficient and safe SWL. The new lens design addresses concomitantly three fundamental drawbacks in EM lithotripters, namely, narrow focal width, nonidealized pulse profile, and significant misalignment in acoustic focus and cavitation activities with the target stone at high output settings. Key design features and performance of the new lens were evaluated using model calculations and experimental measurements against the original lens under comparable acoustic pulse energy (E+) of 40 mJ. The -6-dB focal width of the new lens was enhanced from 7.4 to 11 mm at this energy level, and peak pressure (41 MPa) and maximum cavitation activity were both realigned to be within 5 mm of the lithotripter focus. Stone comminution produced by the new lens was either statistically improved or similar to that of the original lens under various in vitro test conditions and was significantly improved in vivo in a swine model (89% vs. 54%, P = 0.01), and tissue injury was minimal using a clinical treatment protocol. The general principle and associated techniques described in this work can be applied to design improvement of all EM lithotripters.

Waveguide QED: Many-body bound-state effects in coherent and Fock-state scattering from a two-level system

Strong coupling between a two-level system (TLS) and bosonic modes produces dramatic quantum optics effects. We consider a one-dimensional continuum of bosons coupled to a single localized TLS, a system which may be realized in a variety of plasmonic, photonic, or electronic contexts. We present the exact many-body scattering eigenstate obtained by imposing open boundary conditions. Multiphoton bound states appear in the scattering of two or more photons due to the coupling between the photons and the TLS. Such bound states are shown to have a large effect on scattering of both Fock- and coherent-state wave packets, especially in the intermediate coupling-strength regime. We compare the statistics of the transmitted light with a coherent state having the same mean photon number: as the interaction strength increases, the one-photon probability is suppressed rapidly, and the two- and three-photon probabilities are greatly enhanced due to the many-body bound states. This results in non-Poissonian light. © 2010 The American Physical Society.

Intracardiac acoustic radiation force impulse (ARFI) and shear wave imaging in pigs with focal infarctions.

Four pigs, three with focal infarctions in the apical intraventricular septum (IVS) and/or left ventricular free wall (LVFW), were imaged with an intracardiac echocardiography (ICE) transducer. Custom beam sequences were used to excite the myocardium with focused acoustic radiation force (ARF) impulses and image the subsequent tissue response. Tissue displacement in response to the ARF excitation was calculated with a phase-based estimator, and transverse wave magnitude and velocity were each estimated at every depth. The excitation sequence was repeated rapidly, either in the same location to generate 40 Hz M-modes at a single steering angle, or with a modulated steering angle to synthesize 2-D displacement magnitude and shear wave velocity images at 17 points in the cardiac cycle. Both types of images were acquired from various views in the right and left ventricles, in and out of infarcted regions. In all animals, acoustic radiation force impulse (ARFI) and shear wave elasticity imaging (SWEI) estimates indicated diastolic relaxation and systolic contraction in noninfarcted tissues. The M-mode sequences showed high beat-to-beat spatio-temporal repeatability of the measurements for each imaging plane. In views of noninfarcted tissue in the diseased animals, no significant elastic remodeling was indicated when compared with the control. Where available, views of infarcted tissue were compared with similar views from the control animal. In views of the LVFW, the infarcted tissue presented as stiff and non-contractile compared with the control. In a view of the IVS, no significant difference was seen between infarcted and healthy tissue, whereas in another view, a heterogeneous infarction was seen to be presenting itself as non-contractile in systole.

A complex intronic enhancer regulates expression of the CFTR gene by direct interaction with the promoter.

Genes can maintain spatiotemporal expression patterns by long-range interactions between cis-acting elements. The cystic fibrosis transmembrane conductance regulator gene (CFTR) is expressed primarily in epithelial cells. An element located within a DNase I-hypersensitive site (DHS) 10 kb into the first intron was previously shown to augment CFTR promoter activity in a tissue-specific manner. Here, we reveal the mechanism by which this element influences CFTR transcription. We employed a high-resolution method of mapping DHS using tiled microarrays to accurately locate the intron 1 DHS. Transfection of promoter-reporter constructs demonstrated that the element displays classical tissue-specific enhancer properties and can independently recruit factors necessary for transcription initiation. In vitro DNase I footprinting analysis identified a protected region that corresponds to a conserved, predicted binding site for hepatocyte nuclear factor 1 (HNF1). We demonstrate by electromobility shift assays (EMSA) and chromatin immunoprecipitation (ChIP) that HNF1 binds to this element both in vitro and in vivo. Moreover, using chromosome conformation capture (3C) analysis, we show that this element interacts with the CFTR promoter in CFTR-expressing cells. These data provide the first insight into the three- dimensional (3D) structure of the CFTR locus and confirm the contribution of intronic cis-acting elements to the regulation of CFTR gene expression.