Explicit local time-stepping methods for Maxwell’s equations

Explicit local time-stepping methods are derived for time dependent Maxwell equations in conducting and non-conducting media. By using smaller time steps precisely where smaller elements in the mesh are located, these methods overcome the bottleneck caused by local mesh refinement in explicit time integrators. When combined with a finite element discretisation in space with an essentially diagonal mass matrix, the resulting discrete time-marching schemes are fully explicit and thus inherently parallel. In a non-conducting source-free medium they also conserve a discrete energy, which provides a rigorous criterion for stability. Starting from the standard leap-frog scheme, local time-stepping methods of arbitrarily high accuracy are derived for non-conducting media. Numerical experiments with a discontinuous Galerkin discretisation in space validate the theory and illustrate the usefulness of the proposed time integration schemes.     

Functional Additive Mixed Models

We propose an extensive framework for additive regression models for correlated functional responses, allowing for multiple partially nested or crossed functional random effects with flexible correlation structures for, for example, spatial, temporal, or longitudinal functional data. Additionally, our framework includes linear and nonlinear effects of functional and scalar covariates that may vary smoothly over the index of the functional response. It accommodates densely or sparsely observed functional responses and predictors which may be observed with additional error and includes both spline-based and functional principal component-based terms. Estimation and inference in this framework is based on standard additive mixed models, allowing us to take advantage of established methods and robust, flexible algorithms. We provide easy-to-use open source software in the pffr() function for the R package refund. Simulations show that the proposed method recovers relevant effects reliably, handles small sample sizes well, and also scales to larger datasets. Applications with spatially and longitudinally observed functional data demonstrate the flexibility in modeling and interpretability of results of our approach.     

Functional Generalized Additive Models

We introduce the functional generalized additive model (FGAM), a novel regression model for association studies between a scalar response and a functional predictor. We model the link-transformed mean response as the integral with respect to t of F{X(t), t} where F( ·, ·) is an unknown regression function and X(t) is a functional covariate. Rather than having an additive model in a finite number of principal components as by Müller and Yao (2008 Müller, H.G., and Yao, F. (2008), “Functional Additive Models,” Journal of the American Statistical Association, 103, 1534–1544.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), our model incorporates the functional predictor directly and thus our model can be viewed as the natural functional extension of generalized additive models. We estimate F( ·, ·) using tensor-product B-splines with roughness penalties. A pointwise quantile transformation of the functional predictor is also considered to ensure each tensor-product B-spline has observed data on its support. The methods are evaluated using simulated data and their predictive performance is compared with other competing scalar-on-function regression alternatives. We illustrate the usefulness of our approach through an application to brain tractography, where X(t) is a signal from diffusion tensor imaging at position, t, along a tract in the brain. In one example, the response is disease-status (case or control) and in a second example, it is the score on a cognitive test. The FGAM is implemented in R in the refund package. There are additional supplementary materials available online.     

Combinatorics of reaction-network posets

Reaction networks are viewed as derived from ordinary molecular structures related in reactant-product pairs so as to manifest a chemical super-structure. Such super-structures then are candidates for applications in a general combinatoric chemistry. Notable additional characterization of a reaction super-structure occurs when such reaction graphs are directed, as for example when there is progressive substitution (or addition) on a fixed molecular skeleton. Such a set of partially ordered entities is in mathematics termed a poset, which further manifests a number of special properties, as then might be utilized in different applications. Focus on the overall "super-structural" poset goes beyond ordinary molecular structure in attending to how a structure fits into a (reaction) network, and thereby brings an extra "dimension" to conventional stereochemical theory. The possibility that different molecular properties vary smoothly along chains of interconnections in such a super-structure is a natural assumption for a novel approach to molecular property and bioactivity correlations. Different manners to interpolate/extrapolate on a poset network yield quantitative super-structure/activity relationships (QSSARs), with some numerical fits, e.g., for properties of polychlorinated biphenyls (PCBs) seemingly being quite reasonable. There seems to be promise for combinatoric posetic ideas.     

Multiplicity theorems for superlinear elliptic problems

In this paper we study second order elliptic equations driven by the Laplacian and p-Laplacian differential operators and a nonlinearity which is (p-)superlinear (it satisfies the Ambrosetti–Rabinowitz condition). For the p-Laplacian equations we prove the existence of five nontrivial smooth solutions, namely two positive, two negative and a nodal solution. Finally we indicate how in the semilinear case, Morse theory can be used to produce six nontrivial solutions. This paper was completed while the first author was visiting the University of Aveiro as an Invited Scientist. The hospitality and financial support of the host institution are gratefully acknowledged. The second and third authors acknowledge the partial financial support of the Portuguese Foundation for Science and Technology (FCT) under the research project POCI/MAT/55524/2004.     

The Creutz-Taube complex revisited: DFT study of the infrared frequencies

The structural parameters, electronic properties, and infrared frequencies of three binuclear ruthenium complexes, [(NH(3))(5)Ru(pyrazine)Ru(NH(3))(5)](n+), n = 4-6, have been investigated with density functional theory. Structural analysis demonstrates that the structure of the mixed-valence 5+, or [II,III], is not an intermediate of the reduced 4+, or [II,II], and the oxidized 6+, or [III,III], compounds. Electronic structure comparison shows that the Ru d(yz) antibonding orbital is empty when n = 5 and 6 and occupied, when n = 4. The infrared frequencies have been calculated for a sequence of models with increasingly detailed accounts of counterions, ranging from the free Creutz-Taube 5+ ion, over jellium embedded, COnductor-like Screening MOdel (COSMO), to the experimental structure of the triclinic [II,III](tos)(5).4H(2)O (tos = p-toluenesulfonate) crystal. Analysis of the Ru vibrations shows that the spectra for the two symmetry-inequivalent Ru atoms are essentially the same. We find that Ru-Ru modes exist near three well-defined frequencies in the solid: at 145, 285, and 345 cm(-1). Similar results are also obtained for the simplified jellium and COSMO models. The spectral properties of these vibrational correlations testify to the existence of two coupled Ru atoms in the same charge state.     

Molecular dynamics simulation of premelting and melting phase transitions in stoichiometric uranium dioxide

Results of molecular dynamics (MD) simulation of UO2 in a wide temperature range are presented and discussed. A new approach to the calibration of a partly ionic Busing-Ida-type model is proposed. A potential parameter set is obtained reproducing the experimental density of solid UO2 in a wide range of temperatures. A conventional simulation of the high-temperature stoichiometric UO2 on large MD cells, based on a novel fast method of computation of Coulomb forces, reveals characteristic features of a premelting lambda transition at a temperature near to that experimentally observed (T(lambda)=2670 K). A strong deviation from the Arrhenius behavior of the oxygen self-diffusion coefficient was found in the vicinity of the transition point. Predictions for liquid UO2, based on the same potential parameter set, are in good agreement with existing experimental data and theoretical calculations.     

Synthesis,characterization, and cytotoxicity evaluation of high-magnetization multifunctional nanoclusters

The paper presents the synthesis, characterization, and in vitro cytotoxicity tests of Fe₃O₄ magnetic nanoclusters coated with ethylenediaminetetraacetic acid disodium salt (EDTA). Electron microscopy analysis (SEM) evidences that magnetite nanoparticles are closely packed into the clusters stabilized with EDTA with well-defined near spherical shapes and sizes in the range 100–200 nm. From XRD measurements, we determined the mean size of the crystallites inside the magnetic cluster about 36 nm. The saturation magnetization determined for the magnetic clusters stabilized with EDTA has high value, about 81.7 emu/g at 300 K. X-ray photoelectron spectroscopy has been used to determine both the elemental and chemical structure of the magnetic cluster surface. In vitro studies have shown that the magnetic clusters at low doses did not induce toxicity on human umbilical vein endothelial cells or lesions of the cell membrane. In contrast, at high doses, the magnetic clusters increased the lipid peroxidation and reduced the leakage of a cytoplasmic enzyme, lactate dehydrogenase (LDH), in parallel with increasing the antioxidant defense.

Open image in new window

Graphical abstract SEM images of EDTA-coated magnetic clusters (MCs) and the HUVEC viability at different MC doses

     

Explicit dynamics equations of the constrained robotic systems

Recursive matrix relations concerning the kinematics and the dynamics of a constrained robotic system, schematized by several kinematical chains, are established in this paper. Introducing frames and bases, we first analyze the geometrical properties of the mechanism and derive a general set of relations. Kinematics of the vector system of velocities and accelerations for each element of robot are then obtained. Expressed for every independent loop of the robot, useful conditions of connectivity regarding the relative velocities and accelerations are determined for direct or inverse kinematics problem. Based on the general principle of virtual powers, final matrix relations written in a recursive compact form express just the explicit dynamics equations of a constrained robotic system. Establishing active forces or actuator torques in an inverse dynamic problem, these equations are useful in fact for real-time control of a robot.     

Foreword

Journal of Radioanalytical and Nuclear Chemistry -