Advances in three-dimensional geoelectric forward solver techniques

Modern geoelectrical data acquisition systems allow large amounts of data to be collected in a short time. Inversions of such data sets require powerful forward solvers for predicting the electrical potentials. State-of-the-art solvers are typically based on finite elements. Recent developments in numerical mathematics led to direct matrix solvers that allow the equation systems arising from such finite element problems to be solved very efficiently. They are particularly useful for 3-D geoelectrical problems, where many electrodes are involved. Although modern direct matrix solvers include optimized memory saving strategies, their application to realistic, large-scale 3-D problems is still somewhat limited. Therefore, we present two novel techniques that allow the number of gridpoints to be reduced considerably, while maintaining a high solution accuracy. In the areas surrounding an electrode array we attach infinite elements that continue the electrical potentials to infinity. This does not only reduce the number of gridpoints, but also avoids the artificial Dirichlet or mixed boundary conditions that are well known to be the cause of numerical inaccuracies. Our second development concerns the singularity removal in the presence of significant surface topography. We employ a fast multipole boundary element method for computing the singular potentials. This renders unnecessary mesh refinements near the electrodes, which results in substantial savings of gridpoints of up to more than 50 per cent. By means of extensive numerical tests we demonstrate that combined application of infinite elements and singularity removal allows the number of gridpoints to be reduced by a factor of ≈6-10 compared with traditional finite element methods. This will be key for applying finite elements and direct matrix solver techniques to realistic 3-D inversion problem

Tapered N-helical metamaterials with three-fold rotational symmetry as improved circular polarizers

Chiral helix-based metamaterials can potentially serve as compact and broadband circular polarizers. We have recently shown that the physics of structures composed of multiple intertwined helices, so called N-helices with N being an integer multiple of 4, is distinct from that of structures made of single circular helices (N = 1). In particular, undesired circular polarization conversion is strictly eliminated for N = 4 helices arranged on a square lattice. However, the fabrication of such structures for infrared/visible operation wavelengths still poses very significant challenges. Thus, we here revisit the possibility of reducing N from 4 to 3, which would ease micro-fabrication considerably. We show analytically that N = 3 helices arranged on a hexagonal lattice exhibit strictly vanishing circular polarization conversion. N = 3 is the smallest option as N = 2 obviously leads to linear birefringence. To additionally improve the circular-polarizer operation bandwidth and the extinction ratio while maintaining high transmission for the wanted polarization and zero conversion, we also investigate by numerical calculations N = 3 helices with tapered diameter along the helix axis. We find operation bandwidths as large as 2.4 octaves

3-D electrical resistivity tomography using adaptive wavelet parameter grids

We present a novel adaptive model parametrization strategy for the 3-D electrical resistivity tomography problem and demonstrate its capabilities with a series of numerical examples. In contrast to traditional parametrization schemes, which are based on fixed disjoint blocks, we discretize the subsurface in terms of Haar wavelets and adaptively adjust the parametrization as the iterative inversion proceeds. This results in a favourable balance of cell sizes and parameter reliability, that is, in regions where the data constrain the subsurface properties well, our parametrization strategy leads to a fine grid, whereas poorly resolved areas are represented only by a few large blocks. This is documented with eigenvalue analyses and by computing model resolution matrices. During the initial iteration steps, only a few model parameters are involved, which reduces the risk that the regularization dominates the inversion. The algorithm also automatically accounts for non-linear effects caused by pronounced conductivity contrasts. Inside conductive features a finer grid is generated than inside more resistive structures. The automated parameter adaptation is computationally efficient, because the coarsening and refinement subroutines have a nearly linear numerical complexity with respect to the number of model parameters. Because our approach is not tightly coupled to electrical resistivity tomography, it should be straightforward to adapt it to other data type

Advanced finite-element methods for design and analysis of nanooptical structures: Applications

An overview on recent applications of the finite-element method Maxwell-solver JCMsuite to simulation tasks in nanooptics is given. Numerical achievements in the fields of optical metamaterials, plasmonics, photonic crystal fibers, light emitting devices, solar cells, optical lithography, optical metrology, integrated optics, and photonic crystals are summarized

Fate of the Universe, Age of the Universe, Dark Matter, and the Decaying Vacuum Energy

It is shown that in the cosmological models based on a vacuum energy decaying as a^{-2}, where a is the scale factor of the universe, the fate of the universe in regard to whether it will collapse in future or expand forever is determined not by the curvature constant k but by an effective curvature constant k_{eff}. It is argued that a closed universe with k=1 may expand forever, in other words simulate the expansion dynamics of a flat or an open universe because of the possibility that k_{eff}=0 or -1, respectively. Two such models, in one of which the vacuum does not interact with matter and in another of which it does, are studied. It is shown that the vacuum equation of state p_{vac}= -\rho_{vac} may be realized in a decaying vacuum cosmology provided the vacuum interacts wuth matter. The optical depths for gravitational lensing as a function of the matter density and other parameters in the models are calculated at a source redshift of 2. The age of the universe is discussed and shown to be compatible with the new Hipparcos lower limit of 11Gyr. The possibility that a time-varying vacuum energy may serve as dark matter is suggested.Comment: AAS LaTex, 29 pages, published in the Astrophysical Journal, 520, 45, 199

Cosmic String Network Evolution in arbitrary Friedmann-Lemaitre models

We use the velocity-dependent one-scale model by Martins & Shellard to investigate the evolution of a GUT long cosmic string network in arbitrary Friedmann-Lemaitre models. Four representative models are used to show that in general there is no scaling solution. The implications for structure formation are briefly discussed.Comment: 8 pages, 4 postscript figures included, submitted to Phys. Rev.

Belief, Credence and Statistical Evidence

According to the Rational Threshold View, a rational agent believes p if and only if her credence in p is equal to or greater than a certain threshold. One of the most serious challenges for this view is the problem of statistical evidence: statistical evidence is often not sufficient to make an outright belief rational, no matter how probable the target proposition is given such evidence. This indicates that rational belief is not as sensitive to statistical evidence as rational credence. The aim of this paper is twofold. First, we argue that, in addition to playing a decisive role in rationalizing outright belief, non-statistical evidence also plays a preponderant role in rationalizing credence. More precisely, when both types of evidence are present in a context, non-statistical evidence should receive a heavier weight than statistical evidence in determining rational credence. Second, based on this result, we argue that a modified version of the Rational Threshold View can avoid the problem of statistical evidence. We conclude by suggesting a possible explanation of the varying sensitivity to different types of evidence for belief and credence based on the respective aims of these attitudes