Time domain simulations of dynamic river networks

The problem of simulating a river network is considered. A river network is considered to comprise of rivers, dams/lakes as well as weirs. We suggest a numerical approach with specific features that enable the correct representation of these assets. For each river the flow of water is described by the shallow water equations which is a system of hyperbolic partial differential equations and at the junctions of the rivers, suitable coupling conditions, viewed as interior boundary conditions are used to couple the dynamics. A different model for the dams is also presented. Numerical test cases are presented which show that the model is able to reproduce the expected dynamics of the system. Other aspects of the modelling such as rainfall, run-off, overflow/flooding, evaporation, absorption/seepage, bed-slopes, bed friction have not been incorporated in the model due to their specific nature

Measuring Planck beams with planets

Aims. Accurate measurement of the cosmic microwave background (CMB) anisotropy requires precise knowledge of the instrument beam. We explore how well the Planck beams will be determined from observations of planets, developing techniques that are also appropriate for other experiments. Methods. We simulate planet observations with a Planck-like scanning strategy, telescope beams, noise, and detector properties. Then we employ both parametric and non-parametric techniques, reconstructing beams directly from the time-ordered data. With a faithful parameterization of the beam shape, we can constrain certain detector properties, such as the time constants of the detectors, to high precision. Alternatively, we decompose the beam using an orthogonal basis. For both techniques, we characterize the errors in the beam reconstruction with Monte Carlo realizations. For a simplified scanning strategy, we study the impact on estimation of the CMB power spectrum. Finally, we explore the consequences for measuring cosmological parameters, focusing on the spectral index of primordial scalar perturbations, n_s. Results. The quality of the power spectrum measurement will be significantly influenced by the optical modeling of the telescope. In our most conservative case, using no information about the optics except the measurement of planets, we find that a single transit of Jupiter across the focal plane will measure the beam window functions to better than 0.3% for the channels at 100–217 GHz that are the most sensitive to the CMB. Constraining the beam with optical modeling can lead to much higher quality reconstruction. Conclusions. Depending on the optical modeling, the beam errors may be a significant contribution to the measurement systematics for n_s

Emissivity measurements of reflective surfaces at near-millimeter wavelengths

We have developed an instrument for directly measuring the emissivity of reflective surfaces at near-millimeter wavelengths. The thermal emission of a test sample is compared with that of a reference surface, allowing the emissivity of the sample to be determined without heating. The emissivity of the reference surface is determined by one’s heating the reference surface and measuring the increase in emission. The instrument has an absolute accuracy of Δe = 5 x 10^-4 and can reproducibly measure a difference in emissivity as small as Δe = 10^-4 between flat reflective samples. We have used the instrument to measure the emissivity of metal films evaporated on glass and carbon fiber-reinforced plastic composite surfaces. We measure an emissivity of (2.15 ± 0.4) x 10^-3 for gold evaporated on glass and (2.65 ± 0.5) x 10^-3 for aluminum evaporated on carbon fiber-reinforced plastic composite

Deterministic cavity quantum electrodynamics with trapped ions

We have employed radio-frequency trapping to localize a single 40Ca+-ion in a high-finesse optical cavity. By means of laser Doppler cooling, the position spread of the ion's wavefunction along the cavity axis was reduced to 42 nm, a fraction of the resonance wavelength of ionized calcium (λ = 397 nm). By controlling the position of the ion in the optical field, continuous and completely deterministic coupling of ion and field was realized. The precise three-dimensional location of the ion in the cavity was measured by observing the fluorescent light emitted upon excitation in the cavity field. The single-ion system is ideally suited to implement cavity quantum electrodynamics under cw conditions. To this end we operate the cavity on the D3/2–P1/2 transition of 40Ca+ (λ = 866 nm). Applications include the controlled generation of single-photon pulses with high efficiency and two-ion quantum gates

Markov Chain Beam Randomization: a study of the impact of PLANCK beam measurement errors on cosmological parameter estimation

We introduce a new method to propagate uncertainties in the beam shapes used to measure the cosmic microwave background to cosmological parameters determined from those measurements. The method, which we call Markov Chain Beam Randomization, MCBR, randomly samples from a set of templates or functions that describe the beam uncertainties. The method is much faster than direct numerical integration over systematic `nuisance' parameters, and is not restricted to simple, idealized cases as is analytic marginalization. It does not assume the data are normally distributed, and does not require Gaussian priors on the specific systematic uncertainties. We show that MCBR properly accounts for and provides the marginalized errors of the parameters. The method can be generalized and used to propagate any systematic uncertainties for which a set of templates is available. We apply the method to the Planck satellite, and consider future experiments. Beam measurement errors should have a small effect on cosmological parameters as long as the beam fitting is performed after removal of 1/f noise.Comment: 17 pages, 23 figures, revised version with improved explanation of the MCBR and overall wording. Accepted for publication in Astronomy and Astrophysics (to appear in the Planck pre-launch special issue

Hydrothermal synthesis of perovskite and pyrochlore powders of potassium tantalate

Potassium tantalate powders were hydrothermally synthesized at 100 to 200 °C in 4 to 15 M aqueous KOH solutions. A defect pyrochlore, Kta_(2)O_(5)(OH). nH2O (n ≈ 1.4), was obtained at 4 M KOH, but at 7–12 M KOH, this pyrochlore was gradually replaced by a defect perovskite as the stable phase. At 15 M KOH, there was no intermediate pyrochlore, only a defect perovskite, K_(0.85)Ta_(0.92)O_(2.43)(OH)_(0.57) 0.15H_(2)O. Synthesis at higher KOH concentrations led to greater incorporation of protons in the perovskite structures. The potassium vacancies required for charge compensation of incorporated protons could accommodate water molecules in the perovskite structure

Hydrothermal synthesis of KNbO_3 and NaNbO_3 powders

Orthorhombic KNbO_3 and NaNbO_3 powders were hydrothermally synthesized in KOH and NaOH solutions (6.7–15 M) at 150 and 200 °C. An intermediate hexaniobate species formed first before eventually converting to the perovskite phase. For synthesis in KOH solutions, the stability of the intermediate hexaniobate ion increased with decreasing KOH concentrations and temperatures. This led to significant variations in the induction periods and accounted for the large disparity in the mass of recovered powder for different processing parameters. It is also believed that protons were incorporated in the lattice of the as-synthesized KNbO_3 powders as water molecules and hydroxyl ions

MM Algorithms for Geometric and Signomial Programming

This paper derives new algorithms for signomial programming, a generalization of geometric programming. The algorithms are based on a generic principle for optimization called the MM algorithm. In this setting, one can apply the geometric-arithmetic mean inequality and a supporting hyperplane inequality to create a surrogate function with parameters separated. Thus, unconstrained signomial programming reduces to a sequence of one-dimensional minimization problems. Simple examples demonstrate that the MM algorithm derived can converge to a boundary point or to one point of a continuum of minimum points. Conditions under which the minimum point is unique or occurs in the interior of parameter space are proved for geometric programming. Convergence to an interior point occurs at a linear rate. Finally, the MM framework easily accommodates equality and inequality constraints of signomial type. For the most important special case, constrained quadratic programming, the MM algorithm involves very simple updates.Comment: 16 pages, 1 figur