The needle in the hundred square degree haystack: The hunt for binary neutron star mergers with LIGO and Palomar Transient Factory

The Advanced LIGO and Virgo experiments are poised to detect gravitational waves (GWs) directly for the first time this decade. The ultimate prize will be joint observation of a compact binary merger in both gravitational and electromagnetic channels. However, GW sky locations that are uncertain by hundreds of square degrees will pose a challenge. I describe a real-time detection pipeline and a rapid Bayesian parameter estimation code that will make it possible to search promptly for optical counterparts in Advanced LIGO. Having analyzed a comprehensive population of simulated GW sources, we describe the sky localization accuracy that the GW detector network will achieve as each detector comes online and progresses toward design sensitivity. Next, in preparation for the optical search with the intermediate Palomar Transient Factory (iPTF), we have developed a unique capability to detect optical afterglows of gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM). Its comparable error regions offer a close parallel to the Advanced LIGO problem, but Fermi's unique access to MeV-GeV photons and its near all-sky coverage may allow us to look at optical afterglows in a relatively unexplored part of the GRB parameter space. We present the discovery and broadband follow-up observations (X-ray, UV, optical, millimeter, and radio) of eight GBM-IPTF afterglows. Two of the bursts (GRB 130702A / iPTF13bxl and GRB 140606B / iPTF14bfu) are at low redshift (z=0.145 and z = 0.384, respectively), are sub-luminous with respect to "standard" cosmological bursts, and have spectroscopically confirmed broad-line type Ic supernovae. These two bursts are possibly consistent with mildly relativistic shocks breaking out from the progenitor envelopes rather than the standard mechanism of internal shocks within an ultra-relativistic jet. On a technical level, the GBM--IPTF effort is a prototype for locating and observing optical counterparts of GW events in Advanced LIGO with the Zwicky Transient Facility.

Simulation of shock-induced bubble collapse with application to vascular injury in shockwave lithotripsy

Shockwave lithotripsy is a noninvasive medical procedure wherein shockwaves are repeatedly focused at the location of kidney stones in order to pulverize them. Stone comminution is thought to be the product of two mechanisms: the propagation of stress waves within the stone and cavitation erosion. However, the latter mechanism has also been implicated in vascular injury. In the present work, shock-induced bubble collapse is studied in order to understand the role that it might play in inducing vascular injury. A high-order accurate, shock- and interface-capturing numerical scheme is developed to simulate the three-dimensional collapse of the bubble in both the free-field and inside a vessel phantom. The primary contributions of the numerical study are the characterization of the shock-bubble and shock-bubble-vessel interactions across a large parameter space that includes clinical shockwave lithotripsy pressure amplitudes, problem geometry and tissue viscoelasticity, and the subsequent correlation of these interactions to vascular injury. Specifically, measurements of the vessel wall pressures and displacements, as well as the finite strains in the fluid surrounding the bubble, are utilized with available experiments in tissue to evaluate damage potential. Estimates are made of the smallest injurious bubbles in the microvasculature during both the collapse and jetting phases of the bubble's life cycle. The present results suggest that bubbles larger than 1 μm in diameter could rupture blood vessels under clinical SWL conditions.

Gauge Invariant Spectral Cauchy Characteristic Extraction of Gravitational Waves in Computational General Relativity

We present a complete system for Spectral Cauchy characteristic extraction (Spectral CCE). Implemented in C++ within the Spectral Einstein Code (SpEC), the method employs numerous innovative algorithms to efficiently calculate the Bondi strain, news, and flux.

Spectral CCE was envisioned to ensure physically accurate gravitational wave-forms computed for the Laser Interferometer Gravitational wave Observatory (LIGO) and similar experiments, while working toward a template bank with more than a thousand waveforms to span the binary black hole (BBH) problem’s seven-dimensional parameter space.

The Bondi strain, news, and flux are physical quantities central to efforts to understand and detect astrophysical gravitational wave sources within the Simulations of eXtreme Spacetime (SXS) collaboration, with the ultimate aim of providing the first strong field probe of the Einstein field equation.

In a series of included papers, we demonstrate stability, convergence, and gauge invariance. We also demonstrate agreement between Spectral CCE and the legacy Pitt null code, while achieving a factor of 200 improvement in computational efficiency.

Spectral CCE represents a significant computational advance. It is the foundation upon which further capability will be built, specifically enabling the complete calculation of junk-free, gauge-free, and physically valid waveform data on the fly within SpEC.

Baroclinic waves in an air-filled thermally driven rotating annulus.

In this study an experimental investigation of baroclinic waves in air in a differentially heated rotating annulus is presented. Air has a Prandtl number of 0.707, which falls within a previously unexplored region of parameter space for baroclinic instability. The flow regimes encountered include steady waves, periodic amplitude vacillations, modulated amplitude vacillations, and either monochromatic or mixed wave number weak waves, the latter being characterized by having amplitudes less than 5% of the applied temperature contrast. The distribution of these flow regimes in parameter space are presented in a regime diagram. It was found that the progression of transitions between different regimes is, as predicted by recent numerical modeling results, in the opposite sense to that usually found in experiments with high Prandtl number liquids. No hysteresis in the flow type, with respect to variations in the rotation rate, was found in this investigation.

Fluid-structure interaction with mean flow: Over-scattering and unstable resonance growth

It is known theoretically [1-3] that infinitely long fluid loaded plates in mean flow exhibit a range of unusual phenomena in the 'long time' limit. These include convective instability, absolute instability and negative energy waves which are destabilized by dissipation. However, structures are necessarily of finite length and may have discontinuities. Moreover, linear instability waves can only grow over a limited number of cycles before non-linear effects become dominant. We have undertaken an analytical and computational study to investigate the response of finite, discontinuous plates to ascertain if these unusual effects might be realized in practice. Analytically, we take a "wave scattering" [2,4] - as opposed to a "modal superposition" [5] - view of the fluttering plate problem. First, we solve for the scattering coefficients of localized plate discontinuities and identify a range of parameter space, well outside the convective instability regime, where over-scattering or amplified reflection/transmission occurs. These are scattering processes that draw energy from the mean flow into the plate. Next, we use the Wiener-Hopf technique to solve for the scattering coefficients from the leading and trailing edges of a baffled plate. Finally, we construct the response of a finite, baffled plate by a superposition of infinite plate propagating waves continuously scattering off the plate ends and solve for the unstable resonance frequencies and temporal growth rates for long plates. We present a comparison between our computational results and the infinite plate theory. In particular, the resonance response of a moderately sized plate is shown to be in excellent agreement with our long plate analytical predictions. Copyright © 2010 by ASME.

Experimental observation of von Kármán vortices during drop impact

The observation of von Kármán type vortices during the impact of water droplets onto a pool of water is reported. Shadowgraph imaging and laser-sheet visualization are used to document these events. The appearance of these vortices occurs within theoretically predicted regions in a Reynolds-splash number parameter space. In addition, and also in agreement with theoretical predictions, smooth splashing, with vortices absent, is found for smaller Reynolds number. © 2012 American Physical Society.

Experimental observation of von Kármán vortices during drop impact.

The observation of von Kármán type vortices during the impact of water droplets onto a pool of water is reported. Shadowgraph imaging and laser-sheet visualization are used to document these events. The appearance of these vortices occurs within theoretically predicted regions in a Reynolds-splash number parameter space. In addition, and also in agreement with theoretical predictions, smooth splashing, with vortices absent, is found for smaller Reynolds number.

Density stratified environments: The double-tank method

We present an alternative method of producing density stratifications in the laboratory based on the 'double-tank' method proposed by Oster (Sci Am 213:70-76, 1965). We refer to Oster's method as the 'forced-drain' approach, as the volume flow rates between connecting tanks are controlled by mechanical pumps. We first determine the range of density profiles that may be established with the forced-drain approach other than the linear stratification predicted by Oster. The dimensionless density stratification is expressed analytically as a function of three ratios: the volume flow rate ratio n, the ratio of the initial liquid volumes λ and the ratio of the initial densities ψ. We then propose a method which does not require pumps to control the volume flow rates but instead allows the connecting tanks to drain freely under gravity. This is referred to as the 'free-drain' approach. We derive an expression for the density stratification produced and compare our predictions with saline stratifications established in the laboratory using the 'free-drain' extension of Oster's method. To assist in the practical application of our results we plot the region of parameter space that yield concave/convex or linear density profiles for both forced-drain and free-drain approaches. The free-drain approach allows the experimentalist to produce a broad range of density profiles by varying the initial liquid depths, cross-sectional and drain opening areas of the tanks. One advantage over the original Oster approach is that density profiles with an inflexion point can now be established. © 2008 Springer-Verlag.

Control of material transport through pulse shape manipulation-a development toward designer pulses

The variety of laser systems available to industrial laser users is growing and the choice of the correct laser for a material target application is often based on an empirical assessment. Industrial master oscillator power amplifier systems with tuneable temporal pulse shapes have now entered the market, providing enormous pulse parameter flexibility in an already crowded parameter space. In this paper, an approach is developed to design interaction parameters based on observations of material responses. Energy and material transport mechanisms are studied using pulsed digital holography, post process analysis techniques and finite-difference modelling to understand the key response mechanisms for a variety of temporal pulse envelopes incident on a silicon (1/1/1) substrate. The temporal envelope is shown to be the primary control parameter of the source term that determines the subsequent material response and the resulting surface morphology. A double peak energy-bridged temporal pulse shape designed through direct application of holographic imaging data is shown to substantially improve surface quality. © 2014 IEEE.

CALCULATION OF SPIN-REORIENTATION TEMPERATURE IN ND2FE14B

The spin-reorientation phenomenon in Nd2Fe14B has been investigated using an angular dependent free energy approach. A magnetic Hamiltonian which includes the crystal electric field term and the exchange term has been established using realistic band structure results. The temperature dependence of the molecular field is accounted for by introducing the Brillouin function and the magnetic Hamiltonian is diagonalized within the ground state multiplet of the Nd ion. The eigenstates are then used to form the partition function for the free energy. At each temperature, the direction of the molecular field is obtained by searching for the minimum in the angular parameter space of the free energy. Our calculations show that for Nd2Fe14B, the net magnetic anisotropy direction is canted away from the c axis at a temperature close to the experimentally reported spin-reorientation temperature of 150 K. The temperature dependence of the magnetic structure is found to be very sensitive to the size of the second order crystal field parameter B20.

Molecular weight dependence of phase Behavior of PEO/P(EO-b-DMS) blends: Application of Sanchez-Lacombe lattice fluid theory

Molecular weight dependence of phase separation behavior of the Poly (ethylene oxide) (PEO)/Poly(ethylene oxide-block-dimethylsiloxane) (P(EO-b-DMS)) blends was investigated by both experimental and theoretical methods. The cloud point curves of PEO/P(EO-b-DMS) blends were obtained by turbidity method. Based on Sanchez-Lacombe lattice fluid theory (SLLFT), the adjustable parameter, epsilon*(12)/k (quantifying the interaction energy between different components), was evaluated by fitting the experimental data in phase diagrams. To calculate the spinodals, binodals, and the volume changes of mixing for these blends, three modified combining rules of the scaling parameters for the block copolymer were introduced.

Solvent-induced Morphology of the Binary Mixture of Diblock Copolymer in Thin Film: The Block Length and Composition Dependence of Morphology

A series of binary SB blend samples with various overall volume fraction of PS (Phi(PS)) and different discrete distribution of the block length (denoted as d(PS) or d(PB)) were prepared by mixing various asymmetric poly(styrene)-block-poly(butadiene) (SB) block copolymers with a symmetric SB block copolymer. The influences of the external solvent field, composition, and the block length distribution on the morphologies of the blends in the thin films were investigated by atomic force microscopy (AFM) and transmission electron microscopy (TEM). The experimental results revealed that after solvent annealing, the interface of the blend thin films depended mainly on the cooperative effects of the annealing solvent and the inherently interfacial curvature of the blends. Upon exposure to the saturated vapor of cyclohexane, which has preferential affinity for the PB block, a "threshold" of Phi(PS) (approximate 0.635-0.707) was found. Below such threshold, the influence of the annealing solvent played an important role on the interfacial curvature of the blend thin film.

Formalizing Triggers: A Learning Model for Finite Spaces

In a recent seminal paper, Gibson and Wexler (1993) take important steps to formalizing the notion of language learning in a (finite) space whose grammars are characterized by a finite number of parameters. They introduce the Triggering Learning Algorithm (TLA) and show that even in finite space convergence may be a problem due to local maxima. In this paper we explicitly formalize learning in finite parameter space as a Markov structure whose states are parameter settings. We show that this captures the dynamics of TLA completely and allows us to explicitly compute the rates of convergence for TLA and other variants of TLA e.g. random walk. Also included in the paper are a corrected version of GW's central convergence proof, a list of "problem states" in addition to local maxima, and batch and PAC-style learning bounds for the model.

Summarizing Qualitative Behavior from Measurements of NonlinearsCircuits

This report describes a program which automatically characterizes the behavior of any driven, nonlinear, electrical circuit. To do this, the program autonomously selects interesting input parameters, drives the circuit, measures its response, performs a set of numeric computations on the measured data, interprets the results, and decomposes the circuit's parameter space into regions of qualitatively distinct behavior. The output is a two-dimensional portrait summarizing the high-level, qualitative behavior of the circuit for every point in the graph, an accompanying textual explanation describing any interesting patterns observed in the diagram, and a symbolic description of the circuit's behavior which can be passed on to other programs for further analysis.

The Stokes experiment in a foam - a summary

Sponsorship: EPSRC