Use equal intervals, , and the interior piecewise-linear hat functions , which equal one at node , vanish at all other nodes and vanish outside the two adjacent intervals. Write . Testing with and integrating the second derivative by parts gives
The integrals are local: , , and . Nonadjacent functions do not overlap. Thus the consistent mass matrix and stiffness matrix are
Equivalently the explicit nodal equations are
For merely L2 initial data a natural choice is the Galerkin projection ; nodal interpolation is also possible for suitably smooth data. The relation identifies the natural finite-element energy. Differentiating gives .
In mean-field electrodynamics, use Reynolds averaging to write and , with zero-mean fluctuations. Assume the averaging operation commutes with the derivatives and that resolved quantities vary little across an averaging region. Averaging the induction equation then gives
This is the mean-field electromotive force, a vector correlation with units of velocity times magnetic field. The closure problem is to express this correlation in terms of the resolved field.
In the kinematic regime the flow statistics are prescribed, so the response is linear in the mean magnetic field, after independently excited fluctuation transients have decayed. If the mean-field scale is much greater than the turbulent scale , and its evolution is slow compared with the response memory, expand that response in the local mean magnetic field and its gradients:
Without scale separation one generally has a spatially nonlocal and temporally retarded response rather than constant local coefficients.
For homogeneous turbulence the coefficients are independent of position. Invariance under proper rotations makes the rank-two response proportional to , and the isotropic rank-three gradient response proportional to the alternating tensor. Choosing the diffusivity sign convention gives
The alpha effect can regenerate a mean magnetic field, while supplies turbulent magnetic diffusivity. For constant positive , because the magnetic field is solenoidal. Thus the total diffusivity is .
The electromotive vector is a polar vector, whereas the magnetic field is an axial vector. Therefore is a pseudoscalar and changes sign under reflection; is an ordinary scalar. A nonzero isotropic alpha effect requires broken reflection symmetry, such as nonzero kinetic helicity. If “isotropic” includes reflection symmetry as well as rotations, must vanish. The dimensions are and .
The magnetic vector potential gives the upper poloidal magnetic field , and is its toroidal magnetic field; play the same roles below. This representation automatically satisfies . Diffusion acts in both regions. Above the interface the alpha effect converts toroidal magnetic field into poloidal magnetic field, while below it a mean shear converts the poloidal component into toroidal magnetic field. Separating these source regions still allows an alpha-Omega dynamo through diffusion across the interface.
For electromagnetic matching, keep the closure actually represented by these equations. Its retained turbulent mean-field electromotive force is above and zero below: only toroidal-to-poloidal alpha generation is retained. Let the mean velocity be , continuous at the interface, and write the mean electric field as
This yields the specified scalar equations with above and below. In particular, the upper tangential electric component is
where is the scalar toroidal-field amplitude. The other tangential component is .
Assume equal permeability, no imposed singular surface current and ordinary electromagnetic matching. Continuity of the normal magnetic field gives , and tangential magnetic continuity gives and . For the nonzero horizontal Fourier wavenumber used below, normal-field continuity is equivalent to ; the spatially constant potential difference is a gauge freedom. Continuity of , with continuous and , gives the final condition. Thus electromagnetic matching for an alpha-Omega interface is
Continuity of follows from and their time evolution, rather than supplying a fifth independent amplitude condition. A full isotropic mean-field electromotive force would add toroidal alpha induction and an alpha contribution to ; that would be a different bulk and interface problem. Retaining it only in a boundary condition would be inconsistent with the equations being solved here.
Set , without loss of generality for a real wave and its conjugate, and write . Define the two complex vertical wavenumbers and so the upper fields decay as and the lower fields as . The original PDF has a lowercase in the lower phase. The homogeneous upper and lower equations give
Both and have positive real part. Since their squares agree, they must be equal, rather than opposite. Put
so and follow from the equations rather than being assumed in advance.
Suppress the common factor and write
The identity gives the upper forcing equation, while the same operator on gives . Hence
The first and third interface dynamo matching give , . The derivative conditions then give
Thus and . Eliminating the nonzero amplitudes produces the dispersion relation of an equal-diffusivity interface dynamo:
In terms of ,
Each branch has a square root with positive real part and is therefore compatible with vertical localization. The branch with the plus sign is the potentially growing one:
The minus branch has . At the first branch is neutrally oscillatory. The frequency sign follows from the stated convention: positive give propagation in the negative direction. Reversing the coupling sign reverses propagation, while the threshold for the conjugate sign depends on the magnitude of the dynamo number.
Figure 1.
Separated interface-dynamo source layers and the growth threshold at dynamo number 32
. Details.
The equal-diffusivity interface dynamo idealizes a solar configuration in which helical convection regenerates poloidal magnetic field above the tachocline and differential rotation generates strong toroidal magnetic field in or below it. Diffusion communicates the fields between the two regions, closing the feedback loop. This separation is the mechanism emphasized by Parker's interface-wave model. An oscillatory traveling mode gives migrating magnetic patterns and polarity reversals; buoyant toroidal flux can appear as sunspot belts. If is a local latitude coordinate, the direction of their migration depends on the alpha–shear sign and the coordinate orientation.
The magnetic period of this branch is ; polarity reverses after half a period, and activity that depends on field strength rather than sign can repeat on that half-period. This is the qualitative relation between the magnetic solar cycle and its activity cycle. The model does not itself predict the observed period or guarantee equatorward migration: those require actual geometry, transport coefficients and shear signs. Real stars also have diffusivity contrasts, meridional transport, spatially varying regeneration and nonlinear feedback. The kinematic exponential mode describes onset, while saturation and a maintained finite-amplitude cycle require further physics.
In the outer solution, derivative terms are small and , so . To resolve its mismatch with the zero left boundary, set , . The original equation becomes exactly . The left inner variable is therefore , with and as .
Multiplying by gives the first integral of a square-root reaction layer. Matching , fixes its constant:
The positive derivative branch below begins with and approaches the plateau exponentially. Thus the left layer has thickness , not despite the coefficient on the original second derivative. Its small amplitude is essential in the dominant balance.
First assume the usual diffusion regime , with real; the complex-diffusion extension is described below. Substitution of a global normal mode gives
Set to remove the first derivative, then put . The resulting oscillator equation is
Writing gives the Hermite differential equation, . For a power series , its recurrence is
Termination at degree requires . If it does not terminate in the appropriate parity, the series has a growing Gaussian component at one infinity; the solution decaying at one end then fails at the other. The Hermite oscillator quantization can also be seen from the self-adjoint operator , whose square-integrable eigenfunctions are the Hermite functions. Their spectrum is with . The ground state shows that the paper's phrase “positive integer” must include zero.
The global modes of a quadratically confined Ginzburg-Landau equation are therefore
Their spatial functions are proportional to ; the quadratic Gaussian decay dominates the linear drift factor at both infinities. For real coefficients, the most dangerous mode is . The system has global hydrodynamic instability precisely when
Equality is marginal and smaller gives decay of all these modes. The finite-width confinement correction explains why a locally absolutely unstable region need not immediately support a growing global mode.
If complex diffusion is intended, impose the parabolic condition and choose with positive real part. The same formulas continue along the rotated ray, with the corresponding branch of ; instability is determined by the real part of the square bracket. Arbitrary constant , especially negative diffusion, does not justify the asserted whole-line Hermite decay or a well-posed evolution.
Assume the coefficients are sufficiently differentiable in the region considered, with no material discontinuity or caustic at the front. Let describe arrival time, the slowness vector, and the unit elastic wave propagation normal. For any field, its jump is the value immediately behind the front minus that immediately ahead. If the first nonzero jump occurs in the kth time derivative, , write
A smooth background may be added without changing these leading jumps. Since , differentiation expresses spatial jumps in terms of and time-derivative jumps. The restriction lets displacement field and its first time derivative remain continuous and avoids an imposed impulsive source at the front.
Insert this local expansion in . The coefficient of is
For isotropic elasticity the slowness acoustic matrix is . Its longitudinal eigenvalue is . The P branch therefore gives the eikonal equation and polarization
The next coefficient, , is
The matrix on the left is symmetric and annihilates . Project onto to obtain its solvability condition. Put . The two projected terms can be evaluated without dropping gradients of mass density or moduli:
These identities use , hence , and . Their curvature terms cancel, leaving
Multiplying by proves the requested P-wavefront discontinuity transport law,
The wave amplitude here is real for a real jump; a complex wave amplitude uses the corresponding modulus-squared law. An elastic ray follows , with arclength and . If is the infinitesimal cross-sectional area of a neighboring ray tube, . Thus is constant along that tube, the ray-tube conservation for a P-wave jump. This defines the geometrical spreading and shows why the wave amplitude changes with material properties even before attenuation is considered. At a caustic the simple ray-tube wave amplitude is singular and a uniform local description is needed.
Refer all layer wave amplitudes to the lower interface. Let a downward incident SH displacement field have unit wave amplitude there, an upward reflected wave have wave amplitude , and the decaying half-space wave wave amplitude . Interface continuity gives
Solving these equations proves the Love-wave interface reflection phase and transmission wave amplitude:
Thus : the evanescent substrate field has nonzero wave amplitude but no net downward mean flux. At the traction-free top, an upward incident wave and its downward reflection must have canceling derivatives, so their displacement field reflection coefficient is .
A layer wave acquires wave phase during a complete round trip. Constructive interference requires
Taking the tangent gives again . The reflection phases, rather than only geometric travel distance, are essential to the mode condition.
Use the same real elastic stiffness tensor for both fields, including its minor and major symmetries, and let , with outward unit boundary normal . No complex conjugation is needed for the initial bilinear identity. Define
The divergence has two force-balance terms and two cross-gradient terms. The latter cancel because
this follows by exchanging the two strain-gradient factors using major symmetry and the minor pair symmetries. The frequency-domain equations then give
Integrating and applying the divergence theorem proves
This is the Betti identity for elastodynamic fields. Spatially varying stiffness is allowed because its derivatives remain inside the divergence. The original PDF's constitutive contraction is ; the TeX's repeated displacement field index is erroneous.
If two harmonic solutions have the same body force and prescribed boundary traction, their difference has no body force and is traction-free. A nonzero difference would therefore be a free eigenfunction at that frequency. For any frequency outside the countable exceptional set found above, the difference must vanish, so the forced boundary-value solution, when it exists, is unique. At an eigenfrequency, a free eigenfunction can be added to any existing solution, and a resonant forcing may fail to admit a steady harmonic response. This distinguishes uniqueness from existence.
For the full time-dependent problem, let be the difference between two fields with identical forces, tractions and initial displacement field and velocity. It obeys , zero boundary traction, and zero initial data. Its energy uniqueness for traction-driven elasticity follows by direct integration:
Multiply the equation by and integrate by parts. Major symmetry turns the strain term into the derivative of elastic energy, giving . Stable elasticity makes , and gives . Thus equals its initial value zero. Even rigid zero-strain motions cannot survive these initial conditions.
The PDF correctly specifies zero initial displacement field and zero initial velocity; the TeX duplicates displacement field in the second condition. The complete causal displacement field is unique, including when its Fourier representation contains resonant modal contributions fixed by the initial data.
Let be holomorphic on a neighborhood of the closure of a bounded region with piecewise smooth positively oriented boundary, with no boundary zeros and with not identically zero. The argument principle states
where is the zero multiplicity. More generally a closed contour avoiding the zeros gives , and the integral is the winding number of about zero.
To prove it, factor near a zero as with . Then . Remove small disjoint circles about the finitely many zeros in the region. On the remainder is holomorphic, so the Cauchy integral theorem reduces the boundary integral to the sum of the small-circle integrals. Each is because is holomorphic on its small disc. This proves the formula. Along a contour, the imaginary part of is the change in a continuous argument, giving the winding interpretation and the indexed version.
Now choose a closed disc on which is the only zero of the nonconstant limit . Its boundary has . Local uniform convergence gives there for all sufficiently large . The contour functions , , have no boundary zero. Their winding numbers are constant under this homotopy, so the argument principle shows that has exactly the same positive number of zeros, counted with multiplicity, inside the disc. This is the Rouche theorem argument, here derived from the contour count.
Choose a zero in that disc. For every , the compact annulus contains no zero of , so uniform convergence makes nonzero on it eventually. Therefore all of the chosen zeros eventually lie within of . Consequently
The maps are eventually not identically zero, since they converge at any fixed point where is nonzero; thus their zero counts above are legitimate.
A double zero need not persist as a double zero. On the unit disc take , . Its zeros converge to zero but are simple, while vanishes only at zero and . Hence the additional simultaneous derivative condition cannot always be imposed. Local uniform convergence preserves total nearby multiplicity, not its concentration at a single point.
Use the curvature-minus-one normalization of the hyperbolic metric:
The associated distance is . A connected Riemann surface has its complete canonical Poincare metric on a Riemann surface precisely when its universal cover is conformally the unit disc. By the uniformization theorem, the other simply connected covering types are the complex plane and Riemann sphere. Disc automorphisms preserve , as direct differentiation of shows; therefore the metric descends through the deck group. This is the uniformization notion of hyperbolicity, not the different potential-theoretic Green-function classification.
First prove the needed Schwarz lemma. If is holomorphic and , then extends holomorphically at zero. On its modulus is at most . The maximum modulus principle and give and .
Write . Conjugating by the source and target disc automorphisms, fixes zero. Applying the derivative bound gives
Integrate along any piecewise smooth path and then take the infimum of lengths: . This proves the Schwarz-Pick theorem. Contraction means nonexpansion, not a uniform Lipschitz constant strictly smaller than one; automorphisms preserve all distances.
If both and have disc universal covers, lift to through . The disc is simply connected, so the lift exists, and it is holomorphic because the covering maps are local biholomorphisms. Both covering maps are local hyperbolic isometries. The disc contraction therefore descends to and, by the path-length argument, . If only is known to have a hyperbolic metric, a nonconstant forces to have one too: otherwise lift from the plane or sphere universal cover of to the disc cover of . A bounded entire function is constant by the Liouville theorem, and a holomorphic map from the compact sphere into the disc is constant by the maximum modulus principle. Thus a nonconstant map into a hyperbolic surface is distance-nonincreasing between the canonical hyperbolic metrics; a nonhyperbolic source admits only constant such maps.
For the prescribed point, Schwarz-Pick theorem gives . Therefore
is holomorphic after filling the removable singularity at , and . Inverting the target automorphism proves
Conversely every holomorphic with this bound gives an admissible disc map: , its inverse target automorphism remains in the disc, and the denominator cannot vanish.
Fix and put . The possible intermediate values fill , since constant choices , , attain them all. Thus the one-point value region of a holomorphic disc map is . Completing the square in gives its Euclidean centre and radius:
This closed disc lies strictly inside the unit disc because . For it degenerates to the single point . Its boundary is attainable by unimodular constant ; requiring to map into the open disc does not remove that boundary.
If is identically zero there are no exceptions. Otherwise choose with and put . Monotonicity implies thereafter. On this set, , whence
Adding the finite initial interval proves
This derivative growth outside a finite-measure set controls exceptional radii in Nevanlinna estimates. For example let be the smooth spherical characteristic and its increasing area-counting function. Apply the lemma to and : outside the union of two finite-measure sets,
Their logarithms are consequently . Such derivative/area bounds in estimates from the Poisson-Jensen formula give the Nevanlinna logarithmic derivative lemma, outside a finite-length set. The second main theorem uses that estimate as its error term. The elementary growth lemma controls that error and its exceptional set; it is not itself a substitute for the analytic estimates.
We use the truncated second main theorem stated in the preceding solution, together with the fact that a transcendental meromorphic function has . Here is a justification of the latter fact. If there are infinitely many poles, fixing any arbitrarily large finite pole count gives . If there are only finitely many poles, multiply by a polynomial removing them to obtain a transcendental entire function . Its Taylor series has nonzero coefficients of arbitrarily high degrees . Cauchy's coefficient estimate gives ; the subharmonic function Poisson estimate gives . Since , letting be arbitrarily large proves the claimed ratio limit. Thus the second-theorem error is along large nonexceptional radii.
Let count just the simple -points. A multiple -point of degree contributes one to the truncated count and at least two to the full count; a simple point contributes one to each. Hence, with the same positive logarithmic weights for large radii,
If each of the five values had only finitely many simple preimages, each simple counting function would be . Summing and using the first main theorem would give . But the second theorem requires , a contradiction. This proves the stronger result five values force infinitely many simple preimages:
In particular more than one simple zero is necessary, but more than one target with simple zeros is not necessary.
For the last distinction, take a nonsingular Weierstrass elliptic function scaled so that ; it has branch values . These are the four totally ramified Weierstrass values: finite branch preimages have local degree two, and all poles are double. Set
The first four targets are the images of those branch values under the Möbius transformation. Their preimages remain double; at a lattice pole, for example, . The last target corresponds to , where , so every such preimage is simple. The degree-two torus map attains that regular value, and periodicity produces infinitely many preimages in the plane. This nonconstant elliptic function is transcendental because it has infinitely many poles. Exactly one of the five selected finite targets therefore has simple preimages.
There are two necessary qualifications to the printed formula. Critical points satisfy , not . At a point with , the local degree is one and the printed summand is zero. In addition, is finite only when ; constant maps and origin-critical maps need separate treatment.
These are genuine issues even for otherwise ordinary disc maps. The map maps the unit disc into itself and has . Its derivative never equals one in the disc, so the printed count is zero, but it has a ramification point at of multiplicity one. For the correct identity has the additional term , proving that the displayed identity with the literal printed count is false. Also shows the undefined basepoint logarithm in an origin-critical example.
Prove the intended identity first for a nonconstant map with . Let be its hyperbolic derivative density and . Away from derivative zeros, is harmonic, and direct differentiation gives
At a zero of of multiplicity , factor with . Its logarithm contains , whose distributional Laplacian is . Thus
The local map degree is , so this is the required excess-degree weighting.
For a smooth function, differentiating its circular mean and applying the divergence theorem gives . Integrating in gives the logarithmic Green-Jensen formula. Factoring the isolated logarithmic singularities extends it to the present :
Substitution yields the critical-point Jensen identity for the hyperbolic derivative
where . This is the requested formula with the critical-point typo corrected. Initially choose a circle containing no ramification point on its boundary; continuity and integrability of logarithmic boundary singularities extend the result to the remaining radii.
By Schwarz-Pick theorem, . Therefore
The normalized hyperbolic distance is , and
Consequently
One can take and . Since the corrected count is nonnegative in the noncritical-basepoint case, this also proves the growth bound for with the literal printed, identically zero count, though not its false identity.
For completeness, if vanishes to order at zero, put and
Apply the same formula to , which is regular at zero. It gives
The origin contribution is a regularized Jensen term, not the undefined . If one instead omits the origin term, the upper bound gains ; this is bounded for , while the characteristic and off-origin count are bounded on any smaller compact disc. Thus a bound still holds for that nonnegative-count convention after enlarging . For a constant map and , so its characteristic growth bound is trivial; the logarithmic identity and a finite critical-divisor count are not defined. These qualifications resolve all maps allowed by the original opening sentence without asserting an undefined equality.
Regard the space of real matrices as , and use the smooth determinant map. For invertible , factor . Expanding the determinant to first order gives
since the linear coefficient of is the sum of the diagonal entries of . At every with determinant one, taking makes this derivative . Thus one is a regular value, and part (ii) makes a submanifold of dimension .
Its tangent space is the kernel of this differential. At the identity,
Indeed tangent derivatives of curves with constant determinant lie in that kernel, and the local regular-value coordinates show every kernel vector is tangent to the level set.
Let , , and . In coordinate charts centered at and , write the map as . Surjectivity of its derivative gives an invertible minor; reorder the domain coordinates so it is the first one. The map
has invertible derivative at the point. By the inverse function theorem, it is a local coordinate change. In these coordinates the fiber is simply the set whose first coordinates are zero, so it is locally a coordinate plane of dimension . These charts are inherited from , proving
If the fiber is empty, it is the empty submanifold under the usual convention. The no-boundary assumptions avoid any extra boundary conditions in this local model.
A point is a critical point of a smooth map if is not surjective. A critical value is the image of at least one critical point. A regular value is a point for which every point of has surjective differential. A point outside the image is therefore a regular value vacuously. These definitions use rank relative to , not merely whether a derivative vanishes in chosen coordinates.