The quadratic formula gives . For , write this as
Each root of a polynomial has modulus , since . As increases from to , the upper root traverses the upper unit semicircle from to , and the lower root traverses the lower semicircle in the opposite angular direction. At the roots coincide at ; at they coincide at . For intermediate they are complex conjugates.
Figure 1.
The conjugate roots travel along the two unit semicircles
. Arrows show increasing , with endpoints and and the midpoint pair . The figure is an original plot of the solution formula. Details.
The factor has modulus and principal argument . Multiplication by adds to its complex argument, while the real factor halves its modulus. Hence
Taking the principal logarithm for gives
This complex number lies in the second quadrant. Its modulus follows from the sum of squared components, and its principal argument follows by the second-quadrant correction:
All arguments of this fixed add integer multiples of to the displayed one.
If is interpreted as the multivalued complex logarithm, there is instead a family , where and . Then , and the principal argument is for , or for . This records the logarithm-branch choice that the paper leaves implicit.
The complex conjugate is . From the complex sine and cosine identities,
Divide these expressions and multiply by the complex conjugate of the denominator. Its squared modulus is . The numerator simplifies using and the double-angle formulas. Thus
These expressions apply except at the poles , with integer , where the original tangent function is undefined.
With , the expression inside the modulus is . Squaring its nonnegative modulus therefore gives
The locus is an ellipse in the complex plane, centred at , with horizontal semiaxis and vertical semiaxis . Its vertices are and .
Figure 1.
Ellipse obtained from the complex modulus condition
. Details.
Write and let its constant modulus be . If , then . For , differentiating with respect to and , then applying the Cauchy-Riemann equations, gives
Thus
The coefficient matrix has determinant , forcing . The Cauchy-Riemann equations then also give . Therefore are constant on a connected domain and
This proves the constant-modulus holomorphic function result. If a domain is disconnected, the conclusion applies separately to each connected component.
If , then , and the matrix exponential series gives . Also , since . Thus
The assumed diagonalization ensures that these account for all eigenvalues, including multiplicities. In particular, for ,
both of modulus one, in agreement with the unitary matrix property.
For the upper Apollonius circle, let its constant ratio be . Matching the center and radius from part (i) gives
Subtracting their squares yields . Also
The relevant inverse hyperbolic cosine is positive. Now is the real part of a local holomorphic logarithm, hence a harmonic function away from . The point is inside the excluded disk, because , while is below the real axis. Thus neither singularity lies in the physical domain.
On the numerator and denominator have equal modulus, so the logarithm is zero. On the circular boundary its value is . Therefore the harmonic potential between a grounded line and a circular boundary is
This verifies Laplace's equation and both prescribed Dirichlet boundary conditions. It also tends to zero at infinity. A uniqueness assertion on this unbounded domain would require an additional boundedness or infinity condition; the displayed function is the natural decaying solution requested.
Write and square the modulus equation. It becomes
For , completing the square gives the Apollonius circle
Hence
The printed statement needs the exception : equal distances from and give their perpendicular bisector, a straight line rather than a finite-radius circle. For the center lies in the upper half-plane; this is the case needed in part (ii). The original PDF confirms the numerator , correcting the duplicated, corrupted OCR formula.
The prescribed pole lets us choose and . Write . Fixing and gives
so and . Thus
Its determinant is , so it is a genuine Möbius transformation. Direct substitution gives , , and a pole at ; the numerator at that pole is .
To identify the image of the unit disc, put
Then , , and
In particular, its pole is outside the closed unit disc. The identity
shows . For any , solve to get
Its denominator is nonzero, and the same modulus identity for this inverse gives . Thus , and global injectivity of a Möbius transformation proves the required one-to-one correspondence.
The two boundary circles meet at : subtracting their equations gives , and either equation then gives . Send these points to zero and infinity by the Möbius transformation
The point is on the boundary, so is holomorphic throughout . Write . Substituting the inverse gives the exact equivalences
Thus
There is no missing point from the inverse here: lies outside this sector. The image is therefore exactly the stated sector, rather than merely a subset of it.
Figure 1.
The circular crescent becomes a right-angle sector, then a sector of angle pi over three
. Details.
Rotate by to obtain the first quadrant, and take the power using its holomorphic logarithm with . The required map is
Here the principal logarithm works because the expression inside it lies in the first quadrant. It multiplies its argument by , giving , and maps its modulus to , which ranges over all positive real numbers.
For explicit surjectivity and injectivity, given with , define
where uses . Then , so and , and these are inverse operations. Each stage is holomorphic with nonzero derivative, proving that is the requested conformal bijection.
The Jacobian matrix at the origin and its characteristic polynomial are
The fixed-point multipliers are
By the Jury stability criterion, both lie strictly inside the unit circle exactly when
The last condition always holds, since . Therefore the origin is asymptotically stable exactly for among the hyperbolic parameter values.
For completeness, its linear stability throughout the parameter line is as follows. When , the complex-conjugate Floquet multipliers have modulus , so the origin repels. When , their modulus is , so it attracts with oscillation. At the repeated Floquet multiplier is , still attracting. For , both Floquet multipliers are real and in . For , and , so one Floquet multiplier lies in and the other exceeds one: the origin is a saddle fixed point of a map. For , the complex-conjugate Floquet multipliers have modulus , so it repels.
Possible changes at the unit circle are exhausted by , , or a nonreal pair with product one. They give
At the Floquet multipliers are ; the nonzero fixed point branches in part (a) meet the origin, and part (c) identifies a subcritical pitchfork bifurcation. Indeed at the exact restricted map is , which makes the origin unstable despite its neutral linear direction.
At the Floquet multipliers are : this is a 1:4 strong resonance of a planar map, coinciding with the change from repelling to attracting as increases. A generic nonresonant Neimark–Sacker bifurcation conclusion would be unjustified. At there is a double Floquet multiplier ; the matrix has a nontrivial Jordan block, and the two further nonzero fixed point branches meet the origin. Linear stability alone is inconclusive at these unit-modulus boundary values. The repeated Floquet multiplier at is strictly inside the unit circle and gives no bifurcation.
As appropriate for a genuine -step linear multistep method, take . Suppose the characteristic polynomial has a polynomial root with modulus . From and the triangle inequality,
Dividing by the positive number gives , a contradiction. Hence
The nonzero leading coefficient matters: if , the assumed inequality forces all the to vanish, so is the zero polynomial rather than a characteristic polynomial with roots. This degenerate case does not define the method in the question.
The power iteration starts from , normally scaled to unit Euclidean norm, and repeats
Normalization prevents the rapidly growing or decaying magnitudes of from obscuring its direction. The iteration requires .
Because is a diagonalizable matrix, choose an eigenbasis , with corresponding eigenvalues , and expand . First suppose the dominant eigenvalue is simple and . Then
For
the sum in parentheses after is in Big O notation. Consequently the power iteration direction error is : the iterates approach the line spanned by . For a real matrix, a uniquely dominant eigenvalue is real, since otherwise its complex conjugate would have the same modulus. If , approaches one unit eigenvector; if , successive alternate signs. In either case,
This is the sign and phase in power iteration qualification. Also, the necessary assumption is a nonzero coefficient in the eigenbasis, rather than merely a nonzero ordinary inner product with : those tests can differ for a nonsymmetric matrix. If , the method cannot generate the missing component.
If the dominant value has a higher-dimensional eigenspace, the same argument replaces by the nonzero sum of all starting components in that eigenspace. The limiting direction is then that of this sum, which is still an eigenvector for . The usual algorithm assumes ; a strict dominance gap guarantees this whenever another eigenvalue is present.
An estimate retaining the sign of is the Rayleigh quotient
Continuity of this quotient at a nonzero eigenvector proves the limit even when is a nonsymmetric matrix. Alternatively, for a coordinate whose limiting eigenvector component is nonzero. The scaling factor alone tends only to .
Use sinusoidal phasors with physical values , , and . The complex impedance is : its modulus is the voltage-to-current amplitude ratio, and its complex argument is the voltage-minus-current phase difference. Here is the electrical resistance part and the signed electrical reactance; the source's notation differs from the common notation.
The component laws , and give the complex impedances of an ideal resistor, inductor and capacitor:
In a series circuit the same electric current flows through both components, so Kirchhoff's voltage law adds their voltage phasors. The resulting series RC impedance is
For positive , the negative phase difference means that electric current leads the applied electric potential difference.
The Born rule uses the squared modulus of the wavefunction, not its real part. Write . Then
Wavefunction normalization requires
Integrating the probability density over the required interval gives the dimensionless probability
Thus approximately of position measurements fall in this interval.
In sinusoidal steady state, complex impedance is the voltage phasor divided by the current phasor. Its modulus is the amplitude ratio and its argument is the voltage lead. With the convention, , so .
The PDF labels across , not across the complete driven network. If is the current through , then and , where
The resistor-referenced parallel-resonator response is
Writing ,
The maximum occurs at , where and the phase is zero. Below resonance the phase is positive, tending to ; above it the phase is negative, tending to . The amplitude tends to zero at both extremes: the inductor shorts the parallel network at low frequency, and the capacitor does so at high frequency. At resonance their reactive currents cancel, leaving only the resistor current.
Figure 1.
Resistor-referenced parallel-resonator amplitude and phase, normalized without choosing circuit component values
. Details.
If instead denotes the voltage across the complete circuit, its different transfer ratio is , with peak . This distinction follows directly from the printed voltage endpoints.
For electric potential difference and electric current varying as sinusoids, choose the phasor convention
The electrical impedance is
Its modulus gives the voltage-to-current amplitude ratio, and its complex argument gives the voltage-minus-current phase difference. For ideal components,
These follow from Ohm's law , the capacitor law , and the inductor law . In phasor notation a time derivative becomes multiplication by , turning the linear ordinary differential equations of the electrical circuit into algebraic equations. Electrical impedances add in a series circuit; their reciprocals, the electrical admittances, add in a parallel circuit. This computes both amplitudes and phase differences in the steady response to sinusoids.
Classically, a particle with is completely reflected and never enters the raised region. In quantum mechanics the wavefunction penetrates the classically forbidden side as an evanescent wave. Define the tunnelling decay constant
Boundedness at positive infinity excludes , so
Matching the wavefunction and its derivative gives
The reflected probability amplitude has unit modulus, hence quantum reflection probability one. The decaying tail carries zero probability current: this semi-infinite potential step has no transmitted beam at infinity, unlike a finite potential barrier. The sketch shows an exact matched example rather than attempting to assign a single real value to a complex wavefunction.
Figure 1.
Real and imaginary components of a matched wavefunction at a potential step, with exponential decay of the probability density
. Details.
Because the probability density decays as , the requested penetration distance is
For a proton, take , and . Then
The wavefunction amplitude decay length is twice this value. The nonrelativistic equation supplied in the paper is appropriate for this estimate since is small compared with the proton rest energy, approximately . If the original PDF's unbarred is instead read literally as the Planck constant, the displayed differential equation gives . This is a valid calculation for that literal coefficient, but is the physical estimate using the standard Schrödinger equation.
Use physical fields proportional to the real part of . The Drude frequency response is
Thus the sketch is an even Lorentzian line shape, with maximum at zero frequency, half that value at , and an tail. The positive-frequency half is shown in the preceding figure. Squaring the modulus differs from plotting either or the real part of .
For a single carrier species, the plasma frequency and electrical conductivity satisfy
Therefore
The corresponding relaxation time is . This is an angular frequency; the ordinary frequency is . The units of electrical conductivity used in the calculation are , consistently with the SI plasma-frequency relation.