Mechanisms of optical angular momentum transfer to nematic liquid crystalline droplets

A detailed study is presented that evaluates the relative importance of wave plate behavior, scattering processes and absorption phenomena in transferring optical torque from circularly polarized light to optically trapped nematic droplets. A wide range of parameters is considered: droplet diameters between 1 and 15 µm, birefringence values from 0.15 to 0.26 and trapping beam powers from 50 mW to 400 mW. Wave plate behavior is verified through the dependence of torque on droplet diameter and material birefringence. The dependence of the magnitude of the torque on material birefringence confirms the additional importance of the scattering mechanism. Absorption processes are found to be negligible

Some topics in the analysis of spherical data.

This thesis is concerned with the statistical analysis of directions in 3 dimensions. An important reference is the book by Mardia (1972). At the time of publication of this book, the repertoire of spherical distributions used for modelling purposes was rather limited, and there was clearly a need to investigate other possibilities. In the last few years there has been some interest in the 8 parameter family of distributions mentioned by Mardia (1975), which is known as the Fisher-Bingham family. In Chapter 1 an outline of the thesis is given. The Fisher-Bingham family is discussed in Chapter 2, and an effective method for calculating the normalising constant is presented. Attention is then focussed on an interesting 6 parameter subfamily, and a simple rule is given for classifying the distributions in this subfamily according to type (unimodal, bimodal, ’closed curve'). Estimation and inference are then discussed, and the Chapter is concluded with a numerical example. In Chapter 3, the family of bimodal distributions presented in Wood (1982) is described. Other bimodal models are also mentioned briefly. The problem of simulating Fisher-Bingham distributions is considered in Chapter 4. Some inequalities are derived and then used to construct suitable envelopes so that an acceptance-rejection procedure can be used. In Chapter 5, the robust estimation of concentration for a Fisher distribution is considered, and L-estimators of the type suggested by Fisher (1982) are investigated. It is shown that the best of these estimators have desirable all-round properties. Indications are also given as to how these ideas can be adapted to other contexts. Possibilities for further research are mentioned in Chapter 6

Using an Ellipsoid Model to Track and Predict the Evolution and Propagation of Coronal Mass Ejections

We present a method for tracking and predicting the propagation and evolution of coronal mass ejections (CMEs) using the imagers on the STEREO and SOHO satellites. By empirically modeling the material between the inner core and leading edge of a CME as an expanding, outward propagating ellipsoid, we track its evolution in three-dimensional space. Though more complex empirical CME models have been developed, we examine the accuracy of this relatively simple geometric model, which incorporates relatively few physical assumptions, including i) a constant propagation angle and ii) an azimuthally symmetric structure. Testing our ellipsoid model developed herein on three separate CMEs, we find that it is an effective tool for predicting the arrival of density enhancements and the duration of each event near 1 AU. For each CME studied, the trends in the trajectory, as well as the radial and transverse expansion are studied from 0 to ~.3 AU to create predictions at 1 AU with an average accuracy of 2.9 hours.Comment: 18 pages, 11 figure

Operating at the extreme: Estimating the upper yield boundary of winter wheat production in commercial practice

© 2020 The Authors. Wheat farming provides 28.5% of global cereal production. After steady growth in average crop yield from 1950 to 1990, wheat yields have generally stagnated, which prompts the question of whether further improvements are possible. Statistical studies of agronomic parameters such as crop yield have so far exclusively focused on estimating parameters describing the whole of the data, rather than the highest yields specifically. These indicators include the mean or median yield of a crop, or finding the combinations of agronomic traits that are correlated with increasing average yields. In this paper, we take an alternative approach and consider high yields only. We carry out an extreme value analysis of winter wheat yield data collected in England and Wales between 2006 and 2015. This analysis suggests that, under current climate and growing conditions, there is indeed a finite upper bound for winter wheat yield, whose value we estimate to be 17.60 tonnes per hectare. We then refine the analysis for strata defined by either location or level of use of agricultural inputs. We find that there is no statistical evidence for variation of maximal yield depending on location, and neither is there statistical evidence that maximum yield levels are improved by high levels of crop protection and fertilizer use

Nonparametric hypothesis testing for equality of means on the simplex

In the context of data that lie on the simplex, we investigate use of empirical and exponential empirical likelihood, and Hotelling and James statistics, to test the null hypothesis of equal population means based on two independent samples. We perform an extensive numerical study using data simulated from various distributions on the simplex. The results, taken together with practical considerations regarding implementation, support the use of bootstrap-calibrated James statistic

Accuracy and Limitations of Fitting and Stereoscopic Methods to Determine the Direction of Coronal Mass Ejections from Heliospheric Imagers Observations

Using data from the Heliospheric Imagers (HIs) onboard STEREO, it is possible to derive the direction of propagation of coronal mass ejections (CMEs) in addition to their speed with a variety of methods. For CMEs observed by both STEREO spacecraft, it is possible to derive their direction using simultaneous observations from the twin spacecraft and also, using observations from only one spacecraft with fitting methods. This makes it possible to test and compare different analyses techniques. In this article, we propose a new fitting method based on observations from one spacecraft, which we compare to the commonly used fitting method of Sheeley et al. (1999). We also compare the results from these two fitting methods with those from two stereoscopic methods, focusing on 12 CMEs observed simultaneously by the two STEREO spacecraft in 2008 and 2009. We find evidence that the fitting method of Sheeley et al. (1999) can result in significant errors in the determination of the CME direction when the CME propagates outside of 60deg \pm 20 deg from the Sun-spacecraft line. We expect our new fitting method to be better adapted to the analysis of halo or limb CMEs with respect to the observing spacecraft. We also find some evidence that direct triangulation in the HI fields-of-view should only be applied to CMEs propagating approximatively towards Earth (\pm 20deg from the Sun-Earth line). Last, we address one of the possible sources of errors of fitting methods: the assumption of radial propagation. Using stereoscopic methods, we find that at least seven of the 12 studied CMEs had an heliospheric deflection of less than 20deg as they propagated in the HI fields-of-view, which, we believe, validates this approximation.Comment: 17 pages, 6 figures, 2 tables, accepted to Solar Physic

Effect of Solar Wind Drag on the Determination of the Properties of Coronal Mass Ejections from Heliospheric Images

The Fixed-\Phi (F\Phi) and Harmonic Mean (HM) fitting methods are two methods to determine the average direction and velocity of coronal mass ejections (CMEs) from time-elongation tracks produced by Heliospheric Imagers (HIs), such as the HIs onboard the STEREO spacecraft. Both methods assume a constant velocity in their descriptions of the time-elongation profiles of CMEs, which are used to fit the observed time-elongation data. Here, we analyze the effect of aerodynamic drag on CMEs propagating through interplanetary space, and how this drag affects the result of the F\Phi and HM fitting methods. A simple drag model is used to analytically construct time-elongation profiles which are then fitted with the two methods. It is found that higher angles and velocities give rise to greater error in both methods, reaching errors in the direction of propagation of up to 15 deg and 30 deg for the F\Phi and HM fitting methods, respectively. This is due to the physical accelerations of the CMEs being interpreted as geometrical accelerations by the fitting methods. Because of the geometrical definition of the HM fitting method, it is affected by the acceleration more greatly than the F\Phi fitting method. Overall, we find that both techniques overestimate the initial (and final) velocity and direction for fast CMEs propagating beyond 90 deg from the Sun-spacecraft line, meaning that arrival times at 1 AU would be predicted early (by up to 12 hours). We also find that the direction and arrival time of a wide and decelerating CME can be better reproduced by the F\Phi due to the cancellation of two errors: neglecting the CME width and neglecting the CME deceleration. Overall, the inaccuracies of the two fitting methods are expected to play an important role in the prediction of CME hit and arrival times as we head towards solar maximum and the STEREO spacecraft further move behind the Sun.Comment: Solar Physics, Online First, 17 page

Physical and electrical characteristics of EDM debris

AbstractIn EDM, debris plays a key role in the electrical conditions of the discharge gap prior to each spark. Despite this, analysis of debris at all length-scales has not yet been performed, and therefore the nature of debris produced by electrical discharge processes is not fully understood. In this study debris created by the machining of two electrode materials set as negative polarity, silicon and titanium carbide, was centrifuged and imaged using SEM and TEM. From this analysis it was shown that electrode debris is 1nm or lower and up to 10μm in size. Population analysis of the particle size distribution was used to inform an electric field model based on a lattice Boltzmann method framework, simulating the effect of the presence of such debris on the electric field strength. This method is shown to be able to capture the local variation of the electric field and predict qualitatively the correct trend of the electric field strength increasing against the debris concentration. Such data is important for prediction and control of discharge gap size, as well as understanding the impact of a build-up of debris on uncontrolled sparking

Digital Single-Cell Analysis of Plant Organ Development Using 3DCellAtlas

Diverse molecular networks underlying plant growth and development are rapidly being uncovered. Integrating these data into the spatial and temporal context of dynamic organ growth remains a technical challenge. We developed 3DCellAtlas, an integrative computational pipeline that semiautomatically identifies cell types and quantifies both 3D cellular anisotropy and reporter abundance at single-cell resolution across whole plant organs. Cell identification is no less than 97.8% accurate and does not require transgenic lineage markers or reference atlases. Cell positions within organs are defined using an internal indexing system generating cellular level organ atlases where data from multiple samples can be integrated. Using this approach, we quantified the organ-wide cell-type-specific 3D cellular anisotropy driving Arabidopsis thaliana hypocotyl elongation. The impact ethylene has on hypocotyl 3D cell anisotropy identified the preferential growth of endodermis in response to this hormone. The spatiotemporal dynamics of the endogenous DELLA protein RGA, expansin gene EXPA3, and cell expansion was quantified within distinct cell types of Arabidopsis roots. A significant regulatory relationship between RGA, EXPA3, and growth was present in the epidermis and endodermis. The use of single-cell analyses of plant development enables the dynamics of diverse regulatory networks to be integrated with 3D organ growth.</p

Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations

We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyse the factors that limit the size of the integration time step in approaches based on plane-wave expansions. The maximum allowed time step is dictated by the highest frequency components of the fictitious electronic dynamics. These can result either from the large wavevector components of the kinetic energy or from the small wavevector components of the Coulomb potential giving rise to the so called {\it charge sloshing} problem. We show how to eliminate large wavevector instabilities by adopting a preconditioning scheme that is implemented here for the first-time in the context of Car-Parrinello ab-initio molecular dynamics simulations of the ionic motion. We also show how to solve the charge-sloshing problem when this is present. We substantiate our theoretical analysis with numerical tests on a number of different silicon and carbon systems having both insulating and metallic character.Comment: RevTex, 9 figures available upon request, to appear in Phys. Rev.