An application of the finite element method to curve fitting

An application of the Finite Element Method (FEM) to the solution of a geometric problem is shown. The problem is related to curve fitting i.e. pass a curve trough a set of given points even if they are irregularly spaced. Situations where cur ves with cusps can be encountered in the practice and therefore smooth interpolatting curves may be unsuitable. In this paper the possibilities of the FEM to deal with this type of problems are shown. A particular example of application to road planning is discussed. In this case the funcional to be minimized should express the unpleasent effects of the road traveller. Some comparative numerical examples are also given.

Canonical Correlation Analysis for Interpreting Airborne Laser Scanning Metrics along the Lorenz Curve of Tree Size Inequality

Canonical Correlation Analysis for Interpreting Airborne Laser Scanning Metrics along the Lorenz Curve of Tree Size Inequality

Evolutionary computation of forests with Degree- and Role-Constrained Minimum Spanning Trees

Finding the degree-constrained minimum spanning tree (DCMST) of a graph is a widely studied NP-hard problem. One of its most important applications is network design. Here we deal with a new variant of the DCMST problem, which consists of finding not only the degree- but also the role-constrained minimum spanning tree (DRCMST), i.e., we add constraints to restrict the role of the nodes in the tree to root, intermediate or leaf node. Furthermore, we do not limit the number of root nodes to one, thereby, generally, building a forest of DRCMSTs. The modeling of network design problems can benefit from the possibility of generating more than one tree and determining the role of the nodes in the network. We propose a novel permutation-based representation to encode these forests. In this new representation, one permutation simultaneously encodes all the trees to be built. We simulate a wide variety of DRCMST problems which we optimize using eight different evolutionary computation algorithms encoding individuals of the population using the proposed representation. The algorithms we use are: estimation of distribution algorithm, generational genetic algorithm, steady-state genetic algorithm, covariance matrix adaptation evolution strategy, differential evolution, elitist evolution strategy, non-elitist evolution strategy and particle swarm optimization. The best results are for the estimation of distribution algorithms and both types of genetic algorithms, although the genetic algorithms are significantly faster.

Benthic and planktic foraminiferal stable isotope record and coccolith abundance spanning the last 500,000 years of ODP Hole 121-758C

The Indian Summer Monsoon (ISM) is an inter-hemispheric and highly variable ocean-atmosphere-land interaction that directly affects the densely populated Indian subcontinent. Here, we present new records of palaeoceanographic variability that span the last 500,000 years from the eastern equatorial Indian Ocean, a relatively under-sampled area of ISM influence. We have generated carbon and oxygen stable isotope records from three foraminiferal species from Ocean Drilling Program Site 758 (5°N, 90°E) to investigate the oceanographic history of this region. We interpret our resultant Dd18O (surface-thermocline) record of upper water-column stratification in the context of past ISM variability, and compare orbital phase relationships in our Site 758 data to other climate and monsoon proxies in the region. Results suggest that upper water-column stratification at Site 758, which is dominated by variance at precession and half-precession frequencies (23, 19 and 11 ka), is forced by both local (5°N) insolation and ISM winds. In the precession (23 ka) band, stratification minima at Site 758 lag northern hemisphere summer insolation maxima (precession minima) by 9 ka, which is consistent with Arabian Sea ISM phase estimates and suggests a common wind forcing in both regions. This phase implicates a strong sensitivity to both ice volume and southern hemisphere insolation forcing via latent heat export from the southern subtropical Indian Ocean. Additionally, we find evidence of possible overprinting of millennial-scale events during glacial terminations in our stratification record, which suggests an influence of remote abrupt climate events on ISM dynamics.

Automatic curve fitting for interactive display /

"AEC report no. C00-2383-002."

Analysis of algorithms for finding all spanning trees of a graph.

"Supported in part by the following grants: US NSF GJ 217, US NSF GJ 812, and project BUILD."

An O(/E/loglog/V/) algorithm for finding minimum spanning trees /

Includes bibliographical references.

Winter sea surface temperature reconstruction from planktic foraminiferal transfer functions in the northeastern Arabian Sea over the last two millennia

The Indian monsoon system is an important climate feature of the northern Indian Ocean. Small variations of the wind and precipitation patterns have fundamental influence on the societal, agricultural, and economic development of India and its neighboring countries. To understand current trends, sensitivity to forcing, or natural variation, records beyond the instrumental period are needed. However, high-resolution archives of past winter monsoon variability are scarce. One potential archive of such records are marine sediments deposited on the continental slope in the NE Arabian Sea, an area where present-day conditions are dominated by the winter monsoon. In this region, winter monsoon conditions lead to distinctive changes in surface water properties, affecting marine plankton communities that are deposited in the sediment. Using planktic foraminifera as a sensitive and well-preserved plankton group, we first characterize the response of their species distribution on environmental gradients from a dataset of surface sediment samples in the tropical and sub-tropical Indian Ocean. Transfer functions for quantitative paleoenvironmental reconstructions were applied to a decadal-scale record of assemblage counts from the Pakistan Margin spanning the last 2000?years. The reconstructed temperature record reveals an intensification of winter monsoon intensity near the year 100 CE. Prior to this transition, winter temperatures were >1.5°C warmer than today. Conditions similar to the present seem to have established after 450 CE, interrupted by a singular event near 950 CE with warmer temperatures and accordingly weak winter monsoon. Frequency analysis revealed significant 75-, 40-, and 37-year cycles, which are known from decadal- to centennial-scale resolution records of Indian summer monsoon variability and interpreted as solar irradiance forcing. Our first independent record of Indian winter monsoon activity confirms that winter and summer monsoons were modulated on the same frequency bands and thus indicates that both monsoon systems are likely controlled by the same driving force.

Catalyst screening: Refinement of the origin of the volcano curve and its implication in heterogeneous catalysis

Understanding the overall catalytic activity trend for rational catalyst design is one of the core goals in heterogeneous catalysis. In the past two decades, the development of density functional theory (DFT) and surface kinetics make it feasible to theoretically evaluate and predict the catalytic activity variation of catalysts within a descriptor-based framework. Thereinto, the concept of the volcano curve, which reveals the general activity trend, usually constitutes the basic foundation of catalyst screening. However, although it is a widely accepted concept in heterogeneous catalysis, its origin lacks a clear physical picture and definite interpretation. Herein, starting with a brief review of the development of the catalyst screening framework, we use a two-step kinetic model to refine and clarify the origin of the volcano curve with a full analytical analysis by integrating the surface kinetics and the results of first-principles calculations. It is mathematically demonstrated that the volcano curve is an essential property in catalysis, which results from the self-poisoning effect accompanying the catalytic adsorption process. Specifically, when adsorption is strong, it is the rapid decrease of surface free sites rather than the augmentation of energy barriers that inhibits the overall reaction rate and results in the volcano curve. Some interesting points and implications in assisting catalyst screening are also discussed based on the kinetic derivation. Moreover, recent applications of the volcano curve for catalyst design in two important photoelectrocatalytic processes (the hydrogen evolution reaction and dye-sensitized solar cells) are also briefly discussed.

Probabilistic Evaluation of Three-Dimensional Reconstructions from X-Ray Images Spanning a Limited Angle

An important part of computed tomography is the calculation of a three-dimensional reconstruction of an object from series of X-ray images. Unfortunately, some applications do not provide sufficient X-ray images. Then, the reconstructed objects no longer truly represent the original. Inside of the volumes, the accuracy seems to vary unpredictably. In this paper, we introduce a novel method to evaluate any reconstruction, voxel by voxel. The evaluation is based on a sophisticated probabilistic handling of the measured X-rays, as well as the inclusion of a priori knowledge about the materials that the object receiving the X-ray examination consists of. For each voxel, the proposed method outputs a numerical value that represents the probability of existence of a predefined material at the position of the voxel while doing X-ray. Such a probabilistic quality measure was lacking so far. In our experiment, false reconstructed areas get detected by their low probability. In exact reconstructed areas, a high probability predominates. Receiver Operating Characteristics not only confirm the reliability of our quality measure but also demonstrate that existing methods are less suitable for evaluating a reconstruction.

Design and Applications of a Decade-Spanning Terahertz Frequency Comb Spectrometer: Doppler-limited Rotational Spectroscopy of Methanol and Methanol-OD

This thesis details the design and applications of a terahertz (THz) frequency comb spectrometer. The spectrometer employs two offset locked Ti:Sapphire femtosecond oscillators with repetition rates of approximately 80 MHz, offset locked at 100 Hz to continuously sample a time delay of 12.5 ns at a maximum time delay resolution of 15.6 fs. These oscillators emit continuous pulse trains, allowing the generation of a THz pulse train by the master, or pump, oscillator and the sampling of this THz pulse train by the slave, or probe, oscillator via the electro-optic effect. Collecting a train of 16 consecutive THz pulses and taking the Fourier transform of this pulse train produces a decade-spanning frequency comb, from 0.25 to 2.5 THz, with a comb tooth width of 5 MHz and a comb tooth spacing of ~80 MHz. This frequency comb is suitable for Doppler-limited rotational spectroscopy of small molecules. Here, the data from 68 individual scans at slightly different pump oscillator repetition rates were combined, producing an interleaved THz frequency comb spectrum, with a maximum interval between comb teeth of 1.4 MHz, enabling THz frequency comb spectroscopy.

The accuracy of the THz frequency comb spectrometer was tested, achieving a root mean square error of 92 kHz measuring selected absorption center frequencies of water vapor at 10 mTorr, and a root mean square error of 150 kHz in measurements of a K-stack of acetonitrile. This accuracy is sufficient for fitting of measured transitions to a model Hamiltonian to generate a predicted spectrum for molecules of interest in the fields of astronomy and physical chemistry. As such, the rotational spectra of methanol and methanol-OD were acquired by the spectrometer. Absorptions from 1.3 THz to 2.0 THz were compared to JPL catalog data for methanol and the spectrometer achieved an RMS error of 402 kHz, improving to 303 kHz when excluding low signal-to-noise absorptions. This level of accuracy compares favorably with the ~100 kHz accuracy achieved by JPL frequency multiplier submillimeter spectrometers. Additionally, the relative intensity performance of the THz frequency comb spectrometer is linear across the entire decade-spanning bandwidth, making it the preferred instrument for recovering lineshapes and taking absolute intensity measurements in the THz region. The data acquired by the Terahertz Frequency Comb Spectrometer for methanol-OD is of comparable accuracy to the methanol data and may be used to refine the fit parameters for the predicted spectrum of methanol-OD.

Breaking the curve with candels: a Bayesian approach to reveal the non-universality of the dust-attenuation law at high redshift

Dust attenuation affects nearly all observational aspects of galaxy evolution, yet very little is known about the form of the dust-attenuation law in the distant universe. Here, we model the spectral energy distributions of galaxies at z ~ 1.5–3 from CANDELS with rest-frame UV to near-IR imaging under different assumptions about the dust law, and compare the amount of inferred attenuated light with the observed infrared (IR) luminosities. Some individual galaxies show strong Bayesian evidence in preference of one dust law over another, and this preference agrees with their observed location on the plane of infrared excess (IRX, L_TIR/L_UV) and UV slope (β). We generalize the shape of the dust law with an empirical model, A_ λ,σ =E(B-V)k_ λ (λ / λ v)^ σ where k_λ is the dust law of Calzetti et al., and show that there exists a correlation between the color excess E(B-V) and tilt δ with δ =(0.62±0.05)log(E(B-V))+(0.26±0.02). Galaxies with high color excess have a shallower, starburst-like law, and those with low color excess have a steeper, SMC-like law. Surprisingly, the galaxies in our sample show no correlation between the shape of the dust law and stellar mass, star formation rate, or β. The change in the dust law with color excess is consistent with a model where attenuation is caused by scattering, a mixed star–dust geometry, and/or trends with stellar population age, metallicity, and dust grain size. This rest-frame UV-to-near-IR method shows potential to constrain the dust law at even higher redshifts (z>3).