Rehabilitation of the Gauss-Jordan algorithm

Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.     

Images of a Bose-Einstein condensate in position and momentum space

In the Bogoliubov theory, a condensate initially prepared in its ground state described by a stationary Bogoliubov vacuum and later perturbed by a time-dependent potential or interaction strength evolves into a time-dependent excited state which is a dynamical Bogoliubov vacuum. The dynamical vacuum has a simple diagonal form in a time-dependent orthonormal basis of single-particle modes. This diagonal representation leads to a Gaussian probability distribution for possible density-measurement outcomes in position and momentum space.     

Double-layer structural-acoustic coupling for cylindrical shell by using a combination of WDA and BIE

This work formulates the double-layer structural-acoustic coupling problem for cylindrical shell by using a combination of the wave-number domain approach (WDA) and the boundary integral equation (BIE). Expressions for the spectral radial velocity of the outer surface of a finite fluid-filled/submerged (FFS) cylindrical thin shell are formulated by means of the transfer matrix equation in wave-number domain. It is shown that the spectral variables on the inner surface of the shell are related to those on the outer surface of the shell. The far field sound radiation from this kind of shell is numerically evaluated for various fluid cases. An experimental verification is performed, and a good correlation between the theoretical results and the experimental results shows that the theoretical study work in this paper is correct.     

Non-linear Stability of Modulated Fronts¶for the Swift–Hohenberg Equation

We consider front solutions of the Swift–Hohenberg equation ∂ _t u= -(1+ ∂ _x ²)² u + ɛ² u -u ³. These are traveling waves which leave in their wake a periodic pattern in the laboratory frame. Using renormalization techniques and a decomposition into Bloch waves, we show the non-linear stability of these solutions. It turns out that this problem is closely related to the question of stability of the trivial solution for the model problem ∂ _t u(x,t) = ∂ _x ² u (x,t)+(1+tanh(x-ct))u(x,t)+u(x,t) ^p with p>3. In particular, we show that the instability of the perturbation ahead of the front is entirely compensated by a diffusive stabilization which sets in once the perturbation has hit the bulk behind the front. Received: 23 February 2001 / Accepted: 27 August 2001     

Characterization and Computation of Canonical Tight Windows for Gabor Frames

Let (\gnm)_{n,m ? \Zst}(\gnm)_{n,m\in\Zst} be a Gabor frame for \LtR\LtR for given window gg. We show that the window \ho = \SQI g\ho=\SQI g that generates the canonically associated tight Gabor frame minimizes ||g-h||\|g-h\| among all windows hh generating a normalized tight Gabor frame. We present and prove versions of this result in the time domain, the frequency domain, the time-frequency domain, and the Zak transform domain, where in each domain the canonical \ho\ho is expressed using functional calculus for Gabor frame operators. Furthermore, we derive a Wiener--Levy type theorem for rationally oversampled Gabor frames. Finally, a Newton-type method for a fast numerical calculation of \ho\ho is presented. We analyze the convergence behavior of this method and demonstrate the efficiency of the proposed algorithm by some numerical examples.     

Cosmic dust and our origins

The small solid particles in the space between the stars provide the surfaces for the production of many simple and complex molecules. Processes involving the effects of ultraviolet irradiation of the thin (hundredth micron) mantles are shown to produce a wide range of molecules and ions also seen in comets. Some of the more complex ones inferred from laboratory experiments are expected to play an important role in the origin of life. An outline of the chemical evolution of interstellar dust as observed and as studied in the laboratory is presented. Observations of comets are shown to provide substantial evidence for their being fluffy aggregates of interstellar dust as it was in the protosolar nebula, i.e. the interstellar cloud which collapsed to form the solar system. The theory that comets may have brought the progenitors of life to the earth is summarized.     

On the critical exponents of random k‐SAT

There has been much recent interest in the satisfiability of random Boolean formulas. A random k‐SAT formula is the conjunction of m random clauses, each of which is the disjunction of k literals (a variable or its negation). It is known that when the number of variables n is large, there is a sharp transition from satisfiability to unsatisfiability; in the case of 2‐SAT this happens when m/n → 1, for 3‐SAT the critical ratio is thought to be m/n ≈ 4.2. The sharpness of this transition is characterized by a critical exponent, sometimes called ν = ν_k (the smaller the value of ν the sharper the transition). Experiments have suggested that ν₃ = 1.5 ± 0.1. ν₄ = 1.25 ± 0.05, ν₅ = 1.1 ± 0.05, ν₆ = 1.05 ± 0.05, and heuristics have suggested that ν_k → 1 as k → ∞. We give here a simple proof that each of these exponents is at least 2 (provided the exponent is well defined). This result holds for each of the three standard ensembles of random k‐SAT formulas: m clauses selected uniformly at random without replacement, m clauses selected uniformly at random with replacement, and each clause selected with probability p independent of the other clauses. We also obtain similar results for q‐colorability and the appearance of a q‐core in a random graph. © 2002 Wiley Periodicals, Inc. Random Struct. Alg., 21: 182–195, 2002     

Miscibility of bisphenol-A polycarbonate with poly(methyl methacrylate)

The miscibility of bisphenol-A polycarbonate (PC) with poly(methyl methacrylate) (PMMA) has been reexamined using differential scanning calorimetry (DSC) and optical indications for phase separation on heating, i.e., lower critical solution temperature (LCST) behavior. Various methods have been used to prepare the blends including methylene chloride (CH₂Cl₂) and tetrahydrofuran (THF) solution casting, melt mixing, and precipitation of PC and PMMA simultaneously from THF solution by using the nonsolvents methanol and heptane. It is shown that the resulting phase behavior for PC/PMMA blends is strongly affected by the blend preparation method. However, these blends are miscible over the whole blend composition range (unambiguous single composition-dependent T_g's and LCST behavior) when prepared by precipitation from solution using heptane as the nonsolvent. To the contrary, solution-cast and melt-mixed PC/PMMA blends were all phase separated, which may be attributed to the “solvent” effect and LCST behavior, respectively, not discovered in previous reports. Methanol precipitation does not lead to fully mixed blends, which demonstrates the importance of the choice of nonsolvent when using the precipitation method.