Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Calibration of BAS-TR image plate response to high energy (3-300 MeV) carbon ions

The paper presents the calibration of Fuji BAS-TR image plate (IP) response to high energy carbon ions of different charge states by employing an intense laser-driven ion source, which allowed access to carbon energies up to 270 MeV. The calibration method consists of employing a Thomson parabola spectrometer to separate and spectrally resolve different ion species, and a slotted CR-39 solid state detector overlayed onto an image plate for an absolute calibration of the IP signal. An empirical response function was obtained which can be reasonably extrapolated to higher ion energies. The experimental data also show that the IP response is independent of ion charge states.

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach

Complex spatiotemporal haemodynamic response following sensory stimulation in the awake rat

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

The hemodynamic impulse response to a single neural event

This article investigates the relation between stimulus-evoked neural activity and cerebral hemodynamics. Specifically, the hypothesis is tested that hemodynamic responses can be modeled as a linear convolution of experimentally obtained measures of neural activity with a suitable hemodynamic impulse response function. To obtain a range of neural and hemodynamic responses, rat whisker pad was stimulated using brief (less than or equal to2 seconds) electrical stimuli consisting of single pulses (0.3 millisecond, 1.2 mA) combined both at different frequencies and in a paired-pulse design. Hemodynamic responses were measured using concurrent optical imaging spectroscopy and laser Doppler flowmetry, whereas neural responses were assessed through current source density analysis of multielectrode recordings from a single barrel. General linear modeling was used to deconvolve the hemodynamic impulse response to a single "neural event" from the hemodynamic and neural responses to stimulation. The model provided an excellent fit to the empirical data. The implications of these results for modeling schemes and for physiologic systems coupling neural and hemodynamic activity are discussed.

Use of intrinsic modes in biology: Examples of indicial response of pulmonary blood pressure to ± step hypoxia

Recently, a new method to analyze biological nonstationary stochastic variables has been presented. The method is especially suitable to analyze the variation of one biological variable with respect to changes of another variable. Here, it is illustrated by the change of the pulmonary blood pressure in response to a step change of oxygen concentration in the gas that an animal breathes. The pressure signal is resolved into the sum of a set of oscillatory intrinsic mode functions, which have zero “local mean,” and a final nonoscillatory mode. With this device, we obtain a set of “mean trends,” each of which represents a “mean” in a definitive sense, and together they represent the mean trend systematically with different degrees of oscillatory content. Correspondingly, the oscillatory content of the signal about any mean trend can be represented by a set of partial sums of intrinsic mode functions. When the concept of “indicial response function” is used to describe the change of one variable in response to a step change of another variable, we now have a set of indicial response functions of the mean trends and another set of indicial response functions to describe the energy or intensity of oscillations about each mean trend. Each of these can be represented by an analytic function whose coefficients can be determined by a least-squares curve-fitting procedure. In this way, experimental results are stated sharply by analytic functions.

Dynamic, convex, and robust optimization with Bayesian learning for response-guided dosing

Thesis (Ph.D.)--University of Washington, 2016-08

Distinct contrast response functions in striate and extra-striate regions of visual cortex revealed with magnetoencephalography (MEG)

Objective: To spatially and temporally characterise the cortical contrast response function to pattern onset stimuli in humans. Methods: Magnetoencephalography (MEG) was used to investigate the human cortical contrast response function to pattern onset stimuli with high temporal and spatial resolution. A beamformer source reconstruction approach was used to spatially localise and identify the time courses of activity at various visual cortical loci. Results: Consistent with the findings of previous studies, MEG beamformer analysis revealed two simultaneous generators of the pattern onset evoked response. These generators arose from anatomically discrete locations in striate and extra-striate visual cortex. Furthermore, these loci demonstrated notably distinct contrast response functions, with striate cortex increasing approximately linearly with contrast, whilst extra-striate visual cortex followed a saturating function. Conclusions: The generators that underlie the pattern onset visual evoked response arise from two distinct regions in striate and extra-striate visual cortex. Significance: The spatially, temporally and functionally distinct mechanisms of contrast processing within the visual cortex may account for the disparate results observed across earlier studies and assist in elucidating causal mechanisms of aberrant contrast processing in neurological disorders. © 2005 International Federation of Clinical Neurophysiology. Published by Elsevier Ireland Ltd. All rights reserved.

CAUSAL INFERENCE WITH A CONTINUOUS TREATMENT AND OUTCOME: ALTERNATIVE ESTIMATORS FOR PARAMETRIC DOSE-RESPONSE FUNCTIONS WITH APPLICATIONS.

Causal inference with a continuous treatment is a relatively under-explored problem. In this dissertation, we adopt the potential outcomes framework. Potential outcomes are responses that would be seen for a unit under all possible treatments. In an observational study where the treatment is continuous, the potential outcomes are an uncountably infinite set indexed by treatment dose. We parameterize this unobservable set as a linear combination of a finite number of basis functions whose coefficients vary across units. This leads to new techniques for estimating the population average dose-response function (ADRF). Some techniques require a model for the treatment assignment given covariates, some require a model for predicting the potential outcomes from covariates, and some require both. We develop these techniques using a framework of estimating functions, compare them to existing methods for continuous treatments, and simulate their performance in a population where the ADRF is linear and the models for the treatment and/or outcomes may be misspecified. We also extend the comparisons to a data set of lottery winners in Massachusetts. Next, we describe the methods and functions in the R package causaldrf using data from the National Medical Expenditure Survey (NMES) and Infant Health and Development Program (IHDP) as examples. Additionally, we analyze the National Growth and Health Study (NGHS) data set and deal with the issue of missing data. Lastly, we discuss future research goals and possible extensions.

Modelling and studies of the spectral response of some optoelectronic components

Solar radiation takes in today's world, an increasing importance. Different devices are used to carry out spectral and integrated measurements of solar radiation. Thus the sensors can be divided into the fallow types: Calorimetric, Thermomechanical, Thermoelectric and Photoelectric. The first three categories are based on components converting the radiation to temperature (or heat) and then into electrical quantity. On the other hand, the photoelectric sensors are based on semiconductor or optoelectronic elements that when irradiated change their impedance or generate a measurable electric signal. The response function of the sensor element depends not only on the intensity of the radiation but also on its wavelengths. The radiation sensors most widely used fit in the first categories, but thanks to the reduction in manufacturing costs and to the increased integration of electronic systems, the use of the photoelectric-type sensors became more interesting. In this work we present a study of the behavior of different optoelectronic sensor elements. It is intended to verify the static response of the elements to the incident radiation. We study the optoelectronic elements using mathematical models that best fit their response as a function of wavelength. As an input to the model, the solar radiation values are generated with a radiative transfer model. We present a modeling of the spectral response sensors of other types in order to compare the behavior of optoelectronic elements with other sensors currently in use.

The feasibility of vibration analysis as a technique to detect osseointegration of transfemoral implants

One of the main causes of above knee or transfemoral amputation (TFA) in the developed world is trauma to the limb. The number of people undergoing TFA due to limb trauma, particularly due to war injuries, has been increasing. Typically the trauma amputee population, including war-related amputees, are otherwise healthy, active and desire to return to employment and their usual lifestyle. Consequently there is a growing need to restore long-term mobility and limb function to this population. Traditionally transfemoral amputees are provided with an artificial or prosthetic leg that consists of a fabricated socket, knee joint mechanism and a prosthetic foot. Amputees have reported several problems related to the socket of their prosthetic limb. These include pain in the residual limb, poor socket fit, discomfort and poor mobility. Removing the socket from the prosthetic limb could eliminate or reduce these problems. A solution to this is the direct attachment of the prosthesis to the residual bone (femur) inside the residual limb. This technique has been used on a small population of transfemoral amputees since 1990. A threaded titanium implant is screwed in to the shaft of the femur and a second component connects between the implant and the prosthesis. A period of time is required to allow the implant to become fully attached to the bone, called osseointegration (OI), and be able to withstand applied load; then the prosthesis can be attached. The advantages of transfemoral osseointegration (TFOI) over conventional prosthetic sockets include better hip mobility, sitting comfort and prosthetic retention and fewer skin problems on the residual limb. However, due to the length of time required for OI to progress and to complete the rehabilitation exercises, it can take up to twelve months after implant insertion for an amputee to be able to load bear and to walk unaided. The long rehabilitation time is a significant disadvantage of TFOI and may be impeding the wider adoption of the technique. There is a need for a non-invasive method of assessing the degree of osseointegration between the bone and the implant. If such a method was capable of determining the progression of TFOI and assessing when the implant was able to withstand physiological load it could reduce the overall rehabilitation time. Vibration analysis has been suggested as a potential technique: it is a non destructive method of assessing the dynamic properties of a structure. Changes in the physical properties of a structure can be identified from changes in its dynamic properties. Consequently vibration analysis, both experimental and computational, has been used to assess bone fracture healing, prosthetic hip loosening and dental implant OI with varying degrees of success. More recently experimental vibration analysis has been used in TFOI. However further work is needed to assess the potential of the technique and fully characterise the femur-implant system. The overall aim of this study was to develop physical and computational models of the TFOI femur-implant system and use these models to investigate the feasibility of vibration analysis to detect the process of OI. Femur-implant physical models were developed and manufactured using synthetic materials to represent four key stages of OI development (identified from a physiological model), simulated using different interface conditions between the implant and femur. Experimental vibration analysis (modal analysis) was then conducted using the physical models. The femur-implant models, representing stage one to stage four of OI development, were excited and the modal parameters obtained over the range 0-5kHz. The results indicated the technique had limited capability in distinguishing between different interface conditions. The fundamental bending mode did not alter with interfacial changes. However higher modes were able to track chronological changes in interface condition by the change in natural frequency, although no one modal parameter could uniquely distinguish between each interface condition. The importance of the model boundary condition (how the model is constrained) was the key finding; variations in the boundary condition altered the modal parameters obtained. Therefore the boundary conditions need to be held constant between tests in order for the detected modal parameter changes to be attributed to interface condition changes. A three dimensional Finite Element (FE) model of the femur-implant model was then developed and used to explore the sensitivity of the modal parameters to more subtle interfacial and boundary condition changes. The FE model was created using the synthetic femur geometry and an approximation of the implant geometry. The natural frequencies of the FE model were found to match the experimental frequencies within 20% and the FE and experimental mode shapes were similar. Therefore the FE model was shown to successfully capture the dynamic response of the physical system. As was found with the experimental modal analysis, the fundamental bending mode of the FE model did not alter due to changes in interface elastic modulus. Axial and torsional modes were identified by the FE model that were not detected experimentally; the torsional mode exhibited the largest frequency change due to interfacial changes (103% between the lower and upper limits of the interface modulus range). Therefore the FE model provided additional information on the dynamic response of the system and was complementary to the experimental model. The small changes in natural frequency over a large range of interface region elastic moduli indicated the method may only be able to distinguish between early and late OI progression. The boundary conditions applied to the FE model influenced the modal parameters to a far greater extent than the interface condition variations. Therefore the FE model, as well as the experimental modal analysis, indicated that the boundary conditions need to be held constant between tests in order for the detected changes in modal parameters to be attributed to interface condition changes alone. The results of this study suggest that in a clinical setting it is unlikely that the in vivo boundary conditions of the amputated femur could be adequately controlled or replicated over time and consequently it is unlikely that any longitudinal change in frequency detected by the modal analysis technique could be attributed exclusively to changes at the femur-implant interface. Therefore further development of the modal analysis technique would require significant consideration of the clinical boundary conditions and investigation of modes other than the bending modes.

Effect of Oil Price on Emissions Trading Market: VAR Analysis

In this study, we investigated the relationship of European Union carbon dioxide CO2 allowances EUAs prices and oil prices by employing a VAR analysis, Granger causality test and impulse response function. If oil price continues increasing, companies will decrease dependency on fossil fuels because of an increase in energy costs. Therefore, the price of EUAs may be affected by variations in oil prices if the greenhouse gases discharged by the consumption of alternative energy are less than that of fossil fuels. There are no previous studies that investigated these relationships. In this study, we analyzed eight types of EUAs EUA05 to EUA12 with a time series daily data set during 2005-2007 collected from a European Climate Exchange time series data set. Differentiations in these eight types were redemption period. We used the New York Mercantile Exchange light sweet crude price as an oil price. From our examination, we found that only the EUA06 and EUA07 types of EUAs Granger-cause oil prices and vice versa and other six types of EUAs do not Granger-cause oil price. These results imply that the earlier redemption period types of EUAs are more sensitive to oil price. In employing the impulse response function, the results showed that a shock to oil price has a slightly positive effect on all types of EUAs for a very short period. On the other hand, we found that a shock to price of EUA has a slightly negative effect on oil price following a positive effect in only EUA06 and EUA07 types. Therefore, these results imply that fluctuations in EUAs prices and oil prices have little effect on each other. Lastly, we did not consider the substitute energy prices in this study, so we plan to include the prices of coal and natural gas in future analyses.

Modeling of the hemodynamic responses in block design fMRI studies

The hemodynamic response function (HRF) describes the local response of brain vasculature to functional activation. Accurate HRF modeling enables the investigation of cerebral blood flow regulation and improves our ability to interpret fMRI results. Block designs have been used extensively as fMRI paradigms because detection power is maximized; however, block designs are not optimal for HRF parameter estimation. Here we assessed the utility of block design fMRI data for HRF modeling. The trueness (relative deviation), precision (relative uncertainty), and identifiability (goodness-of-fit) of different HRF models were examined and test-retest reproducibility of HRF parameter estimates was assessed using computer simulations and fMRI data from 82 healthy young adult twins acquired on two occasions 3 to 4 months apart. The effects of systematically varying attributes of the block design paradigm were also examined. In our comparison of five HRF models, the model comprising the sum of two gamma functions with six free parameters had greatest parameter accuracy and identifiability. Hemodynamic response function height and time to peak were highly reproducible between studies and width was moderately reproducible but the reproducibility of onset time was low. This study established the feasibility and test-retest reliability of estimating HRF parameters using data from block design fMRI studies.

A Kalman filter based strategy for linear structural system identification based on multiple static and dynamic test data

The problem of identification of stiffness, mass and damping properties of linear structural systems, based on multiple sets of measurement data originating from static and dynamic tests is considered. A strategy, within the framework of Kalman filter based dynamic state estimation, is proposed to tackle this problem. The static tests consists of measurement of response of the structure to slowly moving loads, and to static loads whose magnitude are varied incrementally; the dynamic tests involve measurement of a few elements of the frequency response function (FRF) matrix. These measurements are taken to be contaminated by additive Gaussian noise. An artificial independent variable τ, that simultaneously parameterizes the point of application of the moving load, the magnitude of the incrementally varied static load and the driving frequency in the FRFs, is introduced. The state vector is taken to consist of system parameters to be identified. The fact that these parameters are independent of the variable τ is taken to constitute the set of ‘process’ equations. The measurement equations are derived based on the mechanics of the problem and, quantities, such as displacements and/or strains, are taken to be measured. A recursive algorithm that employs a linearization strategy based on Neumann’s expansion of structural static and dynamic stiffness matrices, and, which provides posterior estimates of the mean and covariance of the unknown system parameters, is developed. The satisfactory performance of the proposed approach is illustrated by considering the problem of the identification of the dynamic properties of an inhomogeneous beam and the axial rigidities of members of a truss structure.

Modelling of solute transport in a mild heterogeneous porous medium using stochastic finite element method: Effects of random source conditions

Randomness in the source condition other than the heterogeneity in the system parameters can also be a major source of uncertainty in the concentration field. Hence, a more general form of the problem formulation is necessary to consider randomness in both source condition and system parameters. When the source varies with time, the unsteady problem, can be solved using the unit response function. In the case of random system parameters, the response function becomes a random function and depends on the randomness in the system parameters. In the present study, the source is modelled as a random discrete process with either a fixed interval or a random interval (the Poisson process). In this study, an attempt is made to assess the relative effects of various types of source uncertainties on the probabilistic behaviour of the concentration in a porous medium while the system parameters are also modelled as random fields. Analytical expressions of mean and covariance of concentration due to random discrete source are derived in terms of mean and covariance of unit response function. The probabilistic behaviour of the random response function is obtained by using a perturbation-based stochastic finite element method (SFEM), which performs well for mild heterogeneity. The proposed method is applied for analysing both the 1-D as well as the 3-D solute transport problems. The results obtained with SFEM are compared with the Monte Carlo simulation for 1-D problems.