Quantile-Based Optimization of Noisy Computer Experiments with Tunable Precision

This article addresses the issue of kriging-based optimization of stochastic simulators. Many of these simulators depend on factors that tune the level of precision of the response, the gain in accuracy being at a price of computational time. The contribution of this work is two-fold: first, we propose a quantile-based criterion for the sequential design of experiments, in the fashion of the classical expected improvement criterion, which allows an elegant treatment of heterogeneous response precisions. Second, we present a procedure for the allocation of the computational time given to each measurement, allowing a better distribution of the computational effort and increased efficiency. Finally, the optimization method is applied to an original application in nuclear criticality safety. This article has supplementary material available online. The proposed criterion is available in the R package DiceOptim.

Associating optical measurements and estimating orbits of geocentric objects through population-based meta-heuristic methods

Currently several thousands of objects are being tracked in the MEO and GEO regions through optical means. The problem faced in this framework is that of Multiple Target Tracking (MTT). In this context both, the correct associations among the observations and the orbits of the objects have to be determined. The complexity of the MTT problem is defined by its dimension S. The number S corresponds to the number of fences involved in the problem. Each fence consists of a set of observations where each observation belongs to a different object. The S ≥ 3 MTT problem is an NP-hard combinatorial optimization problem. There are two general ways to solve this. One way is to seek the optimum solution, this can be achieved by applying a branch-and- bound algorithm. When using these algorithms the problem has to be greatly simplified to keep the computational cost at a reasonable level. Another option is to approximate the solution by using meta-heuristic methods. These methods aim to efficiently explore the different possible combinations so that a reasonable result can be obtained with a reasonable computational effort. To this end several population-based meta-heuristic methods are implemented and tested on simulated optical measurements. With the advent of improved sensors and a heightened interest in the problem of space debris, it is expected that the number of tracked objects will grow by an order of magnitude in the near future. This research aims to provide a method that can treat the correlation and orbit determination problems simultaneously, and is able to efficiently process large data sets with minimal manual intervention.

Streak detection algorithm for space debris detection on optical images.

Any image processing object detection algorithm somehow tries to integrate the object light (Recognition Step) and applies statistical criteria to distinguish objects of interest from other objects or from pure background (Decision Step). There are various possibilities how these two basic steps can be realized, as can be seen in the different proposed detection methods in the literature. An ideal detection algorithm should provide high recognition sensitiv ity with high decision accuracy and require a reasonable computation effort . In reality, a gain in sensitivity is usually only possible with a loss in decision accuracy and with a higher computational effort. So, automatic detection of faint streaks is still a challenge. This paper presents a detection algorithm using spatial filters simulating the geometrical form of possible streaks on a CCD image. This is realized by image convolution. The goal of this method is to generate a more or less perfect match between a streak and a filter by varying the length and orientation of the filters. The convolution answers are accepted or rejected according to an overall threshold given by the ackground statistics. This approach yields as a first result a huge amount of accepted answers due to filters partially covering streaks or remaining stars. To avoid this, a set of additional acceptance criteria has been included in the detection method. All criteria parameters are justified by background and streak statistics and they affect the detection sensitivity only marginally. Tests on images containing simulated streaks and on real images containing satellite streaks show a very promising sensitivity, reliability and running speed for this detection method. Since all method parameters are based on statistics, the true alarm, as well as the false alarm probability, are well controllable. Moreover, the proposed method does not pose any extraordinary demands on the computer hardware and on the image acquisition process.

Rotational constants and structure of para-difluorobenzene determined by femtosecond Raman coherence spectroscopy: A new transient type

Femtosecond Raman rotational coherence spectroscopy (RCS) detected by degenerate four-wave mixing is a background-free method that allows to determine accurate gas-phase rotational constants of non-polar molecules. Raman RCS has so far mostly been applied to the regular coherence patterns of symmetric-top molecules, while its application to nonpolar asymmetric tops has been hampered by the large number of RCS transient types, the resulting variability of the RCS patterns, and the 10³–10⁴ times larger computational effort to simulate and fit rotational Raman RCS transients. We present the rotational Raman RCS spectra of the nonpolar asymmetric top 1,4-difluorobenzene (para-difluorobenzene, p-DFB) measured in a pulsed Ar supersonic jet and in a gas cell over delay times up to ~2.5 ns. p-DFB exhibits rotational Raman transitions with ΔJ = 0, 1, 2 and ΔK = 0, 2, leading to the observation of J −, K −, A −, and C–type transients, as well as a novel transient (S–type) that has not been characterized so far. The jet and gas cell RCS measurements were fully analyzed and yield the ground-state (v = 0) rotational constants Aₒ = 5637.68(20) MHz, Bₒ = 1428.23(37) MHz, and Cₒ = 1138.90(48) MHz (1σ uncertainties). Combining the Aₒ, Bₒ, and Cₒ constants with coupled-cluster with single-, double- and perturbatively corrected triple-excitation calculations using large basis sets allows to determine the semi-experimental equilibrium bond lengths rₑ(C₁–C₂) = 1.3849(4) Å, rₑ(C₂–C³) = 1.3917(4) Å, rₑ(C–F) = 1.3422(3) Å, and rₑ(C₂–H₂) = 1.0791(5) Å.

Computational biomechanics to simulate the femoropopliteal intersection during knee flexion: a preliminary study

To assess if finite element (FE) models can be used to predict deformation of the femoropopliteal segment during knee flexion.

Paternity of subordinates raises cooperative effort in cichlids

Background In cooperative breeders, subordinates generally help a dominant breeding pair to raise offspring. Parentage studies have shown that in several species subordinates can participate in reproduction. This suggests an important role of direct fitness benefits for cooperation, particularly where groups contain unrelated subordinates. In this situation parentage should influence levels of cooperation. Here we combine parentage analyses and detailed behavioural observations in the field to study whether in the highly social cichlid Neolamprologus pulcher subordinates participate in reproduction and if so, whether and how this affects their cooperative care, controlling for the effect of kinship. Methodology/Principal Findings We show that: (i) male subordinates gained paternity in 27.8% of all clutches and (ii) if they participated in reproduction, they sired on average 11.8% of young. Subordinate males sharing in reproduction showed more defence against experimentally presented egg predators compared to subordinates not participating in reproduction, and they tended to stay closer to the breeding shelter. No effects of relatedness between subordinates and dominants (to mid-parent, dominant female or dominant male) were detected on parentage and on helping behaviour. Conclusions/Significance This is the first evidence in a cooperatively breeding fish species that the helping effort of male subordinates may depend on obtained paternity, which stresses the need to consider direct fitness benefits in evolutionary studies of helping behaviour.

Au@Hg nanoalloy formation through direct amalgamation: Structural, spectroscopic, and computational evidence for slow nanoscale diffusion

Dynamic core-shell nanoparticles have received increasing attention in recent years. This paper presents a detailed study of Au-Hg nanoalloys, whose composing elements show a large difference in cohesive energy. A simple method to prepare Au@Hg particles with precise control over the composition up to 15 atom% mercury is introduced, based on reacting a citrate stabilized gold sol with elemental mercury. Transmission electron microscopy shows an increase of particle size with increasing mercury content and, together with X-ray powder diffraction, points towards the presence of a core-shell structure with a gold core surrounded by an Au-Hg solid solution layer. The amalgamation process is described by pseudo-zero-order reaction kinetics, which indicates slow dissolution of mercury in water as the rate determining step, followed by fast scavenging by nanoparticles in solution. Once adsorbed at the surface, slow diffusion of Hg into the particle lattice occurs, to a depth of ca. 3 nm, independent of Hg concentration. Discrete dipole approximation calculations relate the UV-vis spectra to the microscopic details of the nanoalloy structure. Segregation energies and metal distribution in the nanoalloys were modeled by density functional theory calculations. The results indicate slow metal interdiffusion at the nanoscale, which has important implications for synthetic methods aimed at core-shell particles.

Can wild ungulate carcasses provide enough biomass to maintain avian scavenger populations? An empirical assessment using a bio-inspired computational model

Background The reduction in the amount of food available for European avian scavengers as a consequence of restrictive public health policies is a concern for managers and conservationists. Since 2002, the application of several sanitary regulations has limited the availability of feeding resources provided by domestic carcasses, but theoretical studies assessing whether the availability of food resources provided by wild ungulates are enough to cover energetic requirements are lacking. Methodology/Findings We assessed food provided by a wild ungulate population in two areas of NE Spain inhabited by three vulture species and developed a P System computational model to assess the effects of the carrion resources provided on their population dynamics. We compared the real population trend with to a hypothetical scenario in which only food provided by wild ungulates was available. Simulation testing of the model suggests that wild ungulates constitute an important food resource in the Pyrenees and the vulture population inhabiting this area could grow if only the food provided by wild ungulates would be available. On the contrary, in the Pre-Pyrenees there is insufficient food to cover the energy requirements of avian scavenger guilds, declining sharply if biomass from domestic animals would not be available. Conclusions/Significance Our results suggest that public health legislation can modify scavenger population trends if a large number of domestic ungulate carcasses disappear from the mountains. In this case, food provided by wild ungulates could be not enough and supplementary feeding could be necessary if other alternative food resources are not available (i.e. the reintroduction of wild ungulates), preferably in European Mediterranean scenarios sharing similar and socio-economic conditions where there are low densities of wild ungulates. Managers should anticipate the conservation actions required by assessing food availability and the possible scenarios in order to make the most suitable decisions.