Machine learning for transient discovery in Pan-STARRS1 difference imaging

Efficient identification and follow-up of astronomical transients is hindered by the need for humans to manually select promising candidates from data streams that contain many false positives. These artefacts arise in the difference images that are produced by most major ground-based time-domain surveys with large format CCD cameras. This dependence on humans to reject bogus detections is unsustainable for next generation all-sky surveys and significant effort is now being invested to solve the problem computationally. In this paper, we explore a simple machine learning approach to real-bogus classification by constructing a training set from the image data of similar to 32 000 real astrophysical transients and bogus detections from the Pan-STARRS1 Medium Deep Survey. We derive our feature representation from the pixel intensity values of a 20 x 20 pixel stamp around the centre of the candidates. This differs from previous work in that it works directly on the pixels rather than catalogued domain knowledge for feature design or selection. Three machine learning algorithms are trained (artificial neural networks, support vector machines and random forests) and their performances are tested on a held-out subset of 25 per cent of the training data. We find the best results from the random forest classifier and demonstrate that by accepting a false positive rate of 1 per cent, the classifier initially suggests a missed detection rate of around 10 per cent. However, we also find that a combination of bright star variability, nuclear transients and uncertainty in human labelling means that our best estimate of the missed detection rate is approximately 6 per cent.

Dissociative ionization of molecules in intense laser fields

Accurate and efficient grid based techniques for the solution of the time-dependent Schrodinger equation for few-electron diatomic molecules irradiated by intense, ultrashort laser pulses are described. These are based on hybrid finite-difference, Lagrange mesh techniques. The methods are applied in three scenarios, namely H-2(+) with fixed internuclear separation, H-2(+) with vibrating nuclei and H-2 with fixed internuclear separation and illustrative results presented.

Efficient grid treatment of the ionization dynamics of laser-driven H-2(+)

We implement a parallel, time-dependent hybrid finite-difference Lagrange mesh code to model the electron dynamics of the fixed-nuclei hydrogen molecular ion subjected to intense ultrashort laser Pulses, Ionization rates are calculated and compared with results from a previous finite-difference approach and also with published Floquet results. The sensitivity of the results to the gauge describing the electron-field interaction is studied. Visualizations of the evolving wave packets are also presented in which the formation of a stable bound-state resonance is observed.

Reduced-complexity iterative equalization for severe time-dispersive MIMO channels

In this paper, a reduced-complexity soft-interference-cancellation minimum mean-square-error.(SIC-MMSE) iterative equalization method for severe time-dispersive multiple-input-multiple-output (MIMO) channels is proposed. To mitigate the severe time dispersiveness of the channel, a single carrier with cyclic prefix is employed, and the equalization is per-formed in the frequency domain. This simplifies the challenging problem of equalization in MIMO channels due to both the intersymbol interference (ISI) and the coantenna interference (CAI). The proposed iterative algorithm works in two stages. The first stage estimates the transmitted frequency-domain symbols using a low-complexity SIC-MMSE equalizer. The second stage converts the estimated frequency-domain symbols in the time domain and finds their means and variances to incorporate in the SIC-MMSE equalizer in the next iteration. Simulation results show the bit-/symbol-error-rate performance of the SIC-MMSE equalizer, with and without coding, for various modulation schemes.

Quality factor assessment of subwavelength resonant cavities at FIR frequencies

Subwavelength resonators at FIR are presented and studied. The structures consist of 1D cavities formed between a metallized (silver) surface and a metamaterial surface comprising a periodic array of silver patches on a silver-backed silicon substrate. The concept derives from recent discoveries of artificial magnetic conductors (AMC). By studying the currents excited on the metamaterial surface by a normally incident plane wave, the nature of the emerging resonant phenomena and the physical mechanism underlying the AMC operation are investigated. Full wave simulations, based on finite element method and time-domain transmission line modelling technique, have been carried out to demonstrate the effective AMC boundary condition and prove the possibilities for subwavelength cavities. The quality factor of the resonant cavities is assessed as a function of the cavity profile. It is demonstrated that the quality factor drops to about 1/8 of the half-wavelength value for lambda/8 resonant cavity.

Two-dimensional hygrothermal transfer in porous building materials

A two-dimensional mathematical model for evaluating the simultaneous heat and moisture migration in porous building materials was proposed. Vapor content and temperature were chosen as the principal driving potentials. The numerical solution was based on the control volume finite difference technique with fully implicit scheme in time. Two validation experiments were developed in this study. The evolution of transient moisture distributions in both one-dimensional and two-dimensional cases was measured. A comparison between experimental results and those obtained by the numerical model proves that they are fully consistent with each other. The model can be easily integrated into a whole building heat, air and moisture transfer model. Another main advantage of the present numerical method lies in the fact that the required moisture transport properties are comparatively simple and easy to determine.

Modelling the mechanical response of the reed-mouthpiece-lip system of a clarinet. Part II: A lumped model approximation

A non-linear lumped model of the reed-mouthpiece-lip system of a clarinet is formulated, in which the lumped parameters are derived from numerical experiments with a finite-difference simulation based on a distributed reed model. The effective stiffness per unit area is formulated as a function of the pressure signal driving the reed, in order to simulate the effects of the reed bending against the lay, and mass and damping terms are added as a first approximation to the dynamic behaviour of the reed. A discrete-time formulation is presented, and its response is compared to that of the distributed model. In addition, the lumped model is applied in the simulation of clarinet tones, enabling the analysis of the effects of using a pressure-dependent stiffness per unit area on sustained oscillations. The analysed effects and features are in qualitative agreement with players' experiences and experimental results obtained in prior studies.

Time-resolved studies of ultrafast wavepacket dynamics in hydrogen molecules

Recent advances in the study of quantum vibrations and rotations in the fundamental hydrogen molecules are reported. Using the deuterium molecules (D-2(+) and D-2) as exemplars, the application of ultrafast femtosecond pump-probe experiments to study the creation and time-resolved imaging of coherent nuclear wavepackets is discussed. The ability to study the motion of these fundamental molecules in the time-domain is a notable milestone, made possible through the advent of ultrashort intense laser pulses with durations on sub-vibrational (and sub-rotational) timescales. Quantum wavepacket revivals are characterised for both vibrational and rotational degrees of freedom and quantum models are used to provide a detailed discussion of the underlying ultrafast physical dynamics for the specialist and non-specialist alike. (C) 2009 Elsevier B.V. All rights reserved.

The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

Time-resolved resonance Raman spectroscopy

Vibrational Raman spectroscopy is now widely recognized as a useful technique for chemical analysis. It has become increasingly popular for the characterization of stable species since the technology which underpins Raman measurements has matured. Time-resolved Raman spectroscopy has also become established as an excellent method for the characterization of transient chemical species but it is not so widely applied. However, the technical advances which have reduced the cost and increased the reliability of conventional: Raman systems can also be exploited in studies of transient species. In some cases it is just as straightforward to record the Raman-spectra of a short-lived transient species as it is to monitor a more stable sample. This raises the possibility of routinely adding time-domain Raman measurements to more conventional Raman techniques, increasing the selectivity of the analysis while retaining its ability to provide spectral information which is characteristic of the species under investigation.

Thermal fatigue analysis of solar panel structure for micro-satellite applications

Thermal fatigue analysis based on 2D finite difference and 3D finite element methods is carried out to study the performance of solar panel structure during micro-satellite life time. Solar panel primary structure consists of honeycomb structure and composite laminates. The 2D finite difference (I-DEAS) model yields predictions of the temperature profile during one orbit. Then, 3D finite element analysis (ANSYS) is applied to predict thermal fatigue damage of solar panel structure. Meshing the whole structure with 2D multi-layer shell elements with sandwich option is not efficient, as it misses thermal response of the honeycomb structure. So we applied a mixed approach between 3D solid and 2D shell elements to model the solar panel structure without the sandwich option.

Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.