1,292 research outputs found
Sparse Bayesian Nonlinear System Identification using Variational Inference
IEEE Bayesian nonlinear system identification for one of the major classes of dynamic model, the nonlinear autoregressive with exogenous input (NARX) model, has not been widely studied to date. Markov chain Monte Carlo (MCMC) methods have been developed, which tend to be accurate but can also be slow to converge. In this contribution, we present a novel, computationally efficient solution to sparse Bayesian identification of the NARX model using variational inference, which is orders of magnitude faster than MCMC methods. A sparsity-inducing hyper-prior is used to solve the structure detection problem. Key results include: 1. successful demonstration of the method on low signal-to-noise ratio signals (down to 2dB); 2. successful benchmarking in terms of speed and accuracy against a number of other algorithms: Bayesian LASSO, reversible jump MCMC, forward regression orthogonalisation, LASSO and simulation error minimisation with pruning; 3. accurate identification of a real world system, an electroactive polymer; and 4. demonstration for the first time of numerically propagating the estimated nonlinear time-domain model parameter uncertainty into the frequency-domain
Sparse Bayesian Nonlinear System Identification using Variational Inference
IEEE Bayesian nonlinear system identification for one of the major classes of dynamic model, the nonlinear autoregressive with exogenous input (NARX) model, has not been widely studied to date. Markov chain Monte Carlo (MCMC) methods have been developed, which tend to be accurate but can also be slow to converge. In this contribution, we present a novel, computationally efficient solution to sparse Bayesian identification of the NARX model using variational inference, which is orders of magnitude faster than MCMC methods. A sparsity-inducing hyper-prior is used to solve the structure detection problem. Key results include: 1. successful demonstration of the method on low signal-to-noise ratio signals (down to 2dB); 2. successful benchmarking in terms of speed and accuracy against a number of other algorithms: Bayesian LASSO, reversible jump MCMC, forward regression orthogonalisation, LASSO and simulation error minimisation with pruning; 3. accurate identification of a real world system, an electroactive polymer; and 4. demonstration for the first time of numerically propagating the estimated nonlinear time-domain model parameter uncertainty into the frequency-domain
Transient HDO rovibrational satellite peaks in solid parahydrogen: evidence of hydrogen atoms or vacancies?
We present FTIR studies of the 193 nm photolysis of fully deuterated formic acid (DCOOD) isolated in solid
parahydrogen at 1.9 K which show evidence of transient HDO rovibrational satellite peaks. The S1 and S2 satellite
peaks are readily detected for α-type (1₀₁ ← 0₀₀) rovibrational transitions of HDO either during or immediately
after photolysis. Intensity measurements show the HDO b-type (1₁₁ ← 0₀₀) rovibrational transitions have
satellite peaks as well, but due to the greater linewidth of these absorptions, the satellite peaks cannot be spectroscopically
resolved from the monomer transition and are therefore difficult to detect. These newly identified
HDO satellite peaks may result from the HDO photoproduct being formed next to an H atom or a vacancy in the
parahydrogen solid. The development of the infrared spectroscopy of these satellite peaks can provide a new
means to study radiation effects on low-temperature hydrogen solids doped with chemical species
Frequency-domain analysis for nonlinear systems with time-domain model parameter uncertainty
Frequency-domain analysis of dynamic systems is important across many areas of engineering. However, whilst there are many analysis methods for linear systems, the problem is much less widely studied for nonlinear systems. Frequency-domain analysis of nonlinear systems using frequency response functions (FRFs) is particularly important to reveal resonances, super/sub-harmonics and energy transfer across frequencies. In this paper the novel contribution is a time-domain model-based approach to describing the uncertainty of nonlinear systems in the frequency-domain. The method takes a nonlinear input-output model that has normally distributed parameters, and propagates that uncertainty into the frequency-domain using analytic expressions based on FRFs. We demonstrate the approach on both synthetic examples of nonlinear systems and a real-world nonlinear system identified from experimental data. We benchmark the proposed approach against a brute-force technique based on Monte Carlo sampling and show that there is good agreement between the methods
Trapping of Projectiles in Fixed Scatterer Calculations
We study multiple scattering off nuclei in the closure approximation. Instead
of reducing the dynamics to one particle potential scattering, the scattering
amplitude for fixed target configurations is averaged over the target
groundstate density via stochastic integration. At low energies a strong
coupling limit is found which can not be obtained in a first order optical
potential approximation. As its physical explanation, we propose it to be
caused by trapping of the projectile. We analyse this phenomenon in mean field
and random potential approximations.
(PACS: 24.10.-i)Comment: 15 page
Radiative charge transfer lifetime of the excited state of (NaCa)
New experiments were proposed recently to investigate the regime of cold
atomic and molecular ion-atom collision processes in a special hybrid
neutral-atom--ion trap under high vacuum conditions. The collisional cooling of
laser pre-cooled Ca ions by ultracold Na atoms is being studied. Modeling
this process requires knowledge of the radiative lifetime of the excited
singlet A state of the (NaCa) molecular system. We calculate
the rate coefficient for radiative charge transfer using a semiclassical
approach. The dipole radial matrix elements between the ground and the excited
states, and the potential curves were calculated using Complete Active Space
Self-Consistent field and M\"oller-Plesset second order perturbation theory
(CASSCF/MP2) with an extended Gaussian basis, 6-311+G(3df). The semiclassical
charge transfer rate coefficient was averaged over a thermal Maxwellian
distribution. In addition we also present elastic collision cross sections and
the spin-exchange cross section. The rate coefficient for charge transfer was
found to be cm/sec, while those for the elastic and
spin-exchange cross sections were found to be several orders of magnitude
higher ( cm/sec and cm/sec,
respectively). This confirms our assumption that the milli-Kelvin regime of
collisional cooling of calcium ions by sodium atoms is favorable with the
respect to low loss of calcium ions due to the charge transfer.Comment: 4 pages, 5 figures; v.2 - conceptual change
Bcc He as a Coherent Quantum Solid
In this work we investigate implications of the quantum nature of bcc %
He. We show that it is a unique solid phase with both a lattice structure and
an Off-Diagonal Long Range Order of coherently oscillating local electric
dipole moments. These dipoles arise from the local motion of the atoms in the
crystal potential well, and oscillate in synchrony to reduce the dipolar
interaction energy. The dipolar ground-state is therefore found to be a
coherent state with a well defined global phase and a three-component complex
order parameter. The condensation energy of the dipoles in the bcc phase
stabilizes it over the hcp phase at finite temperatures. We further show that
there can be fermionic excitations of this ground-state and predict that they
form an optical-like branch in the (110) direction. A comparison with
'super-solid' models is also discussed.Comment: 12 pages, 8 figure
Cardiovascular disease diagnoses among older women with endometrial cancer
Background: Endometrial cancer (EC) shares risk factors (e.g. obesity) with cardiovascular disease (CVD), yet little research has investigated CVD diagnoses among EC survivors. We aimed to describe the burden of CVD diagnoses among older women with EC compared to women without a cancer history. Methods: Women aged 66+ years with an EC diagnosis during 2004–2017 (N = 44,386) and matched women without cancer (N = 221,219) were identified in the SEER-Medicare linked data. An index date was defined as the cancer diagnosis date of the EC case in that matched set. ICD-9/10 diagnosis codes were used to define CVD outcomes in the Medicare claims. Prevalent CVD was identified using diagnosis codes in the year before the index date. Hazard ratios (HRs) for incident CVD diagnoses after the index date were estimated using multivariable Cox proportional hazards regression. Women with a prevalent CVD were excluded from incidence analyses for that outcome. Results: Compared to women without cancer, women with EC had a higher prevalence of CVD diagnoses at the index date. In analyses beginning follow-up at 1 year post-index date, EC survivors had an increased risk of incident CVD diagnoses including ischemic heart diseases (HR = 1.73; 95% CI: 1.69–1.78), pulmonary heart disease (HR = 1.95; 95% CI: 1.88–2.02), and diseases of the veins and lymphatics (HR = 2.71; 95% CI: 95% CI: 2.64–2.78). Risk of CVD diagnoses among women with EC was also elevated within the first year post-index date. Conclusions: Management of pre-existing CVD and monitoring for incident CVD may be critical during EC treatment and throughout long-term survivorship
Species abundance dynamics under neutral assumptions: a Bayesian approach to the controversy
1. Hubbell's 'Unified Neutral Theory of Biodiversity and Biogeography' (UNTB) has generated much controversy about both the realism of its assumptions and how well it describes the species abundance dynamics in real communities. 2. We fit a discrete-time version of Hubbell's neutral model to long-term macro-moth (Lepidoptera) community data from the Rothamsted Insect Survey (RIS) light-traps network in the United Kingdom. 3. We relax the assumption of constant community size and use a hierarchical Bayesian approach to show that the model does not fit the data well as it would need parameter values that are impossible. 4. This is because the ecological communities fluctuate more than expected under neutrality. 5. The model, as presented here, can be extended to include environmental stochasticity, density-dependence, or changes in population sizes that are correlated between different species
- …