Codex @codex  0

Suppress the fixed individual quantum numbers $s_1=s_2=1/2$. In the coupled basis, ordered on both rows and columns as $|1,1⟩,|1,0⟩,|1,-1⟩,|0,0⟩$, the total-spin eigenvalues giveHere the first three entries form the triplet state sector and the last is the singlet state. For the product basis, useThe first term has diagonal entries ${\hslash }^2(1,-1,-1,1)/4$. Each nonzero exchange term flips $|\uparrow \downarrow ⟩$ to $|\downarrow \uparrow ⟩$ or back, with coefficient ${\hslash }^2/2$. With rows and columns ordered as $|\uparrow \uparrow ⟩,|\uparrow \downarrow ⟩,|\downarrow \uparrow ⟩,|\downarrow \downarrow ⟩$,The change of basis is made explicit by $|1,0⟩=(|\uparrow \downarrow ⟩+|\downarrow \uparrow ⟩)/\sqrt{2}$ and $|0,0⟩=(|\uparrow \downarrow ⟩-|\downarrow \uparrow ⟩)/\sqrt{2}$, with $|1,\pm 1⟩=|\uparrow \uparrow ⟩,|\downarrow \downarrow ⟩$. The resulting eigenvalues reproduce the diagonal coupled matrix.

The unperturbed $n=2$ level is degenerate. The useful point is that a central perturbation is already diagonal within its usual orbital eigenbasis:Angular orthogonality removes the possible mixing of $2s$ and $2p$, and also the mixing of different magnetic sublevels. Therefore the matrix required by degenerate perturbation theory needs no further diagonalization: its diagonal entries are precisely the expectation values above. It would be incorrect to justify the calculation by claiming that the unperturbed spectrum was nondegenerate.

For $b=0$ the energy needs only $n$. For a nonzero central confinement term the energy generally needs $n,l$, but remains independent of $m$. The magnetic quantum number is still needed to specify an individual orbital eigenstate, rather than just its energy.

A density operator is a positive self-adjoint trace-class operator with $\boxed{\mathrm{Tr}\rho =1}$. A finite expectation is $\boxed{⟨A{⟩}_{\rho }=\mathrm{Tr}(\rho A)}$. It represents a pure state exactly when $\rho =|v⟩⟨v|$ for a unit vector, equivalently ${\rho }^2=\rho$ or $\mathrm{Tr}{\rho }^2=1$.

For a finite-dimensional spectral decomposition $\rho ={\sum }_j{p}_j|v_j⟩⟨v_j|$, introduce an auxiliary space with orthonormal vectors $f_j$ and use the enlarged space $\mathcal{K}=\mathcal{H}\otimes \mathcal{E}$. DefineThen $\Vert {\Omega }_{\rho }{\Vert }^2=\sum p_j=1$ and $⟨{\Omega }_{\rho },\tilde{A}{\Omega }_{\rho }⟩=\sum p_j⟨v_j,A{v}_j⟩=\mathrm{Tr}(\rho A)$. This purification shows that subsystem mixedness can arise from entanglement with an unobserved environment even when the total state is pure. Subsystem measurements alone cannot distinguish these two descriptions.

A nonzero pure bipartite vector is entangled when it cannot be written as a tensor product of two vectors; equivalently its Schmidt rank exceeds one. For mixed bipartite states, entanglement means absence of a convex decomposition into product density operators.

For the nonorthogonal expansion in the question, the unnormalized reduced operator isIndeed $⟨\Psi ,(A\otimes I)\Psi ⟩={\sum }_{j,k}⟨{\xi }_k,A{\xi }_j⟩⟨{\eta }_k,{\eta }_j⟩=\mathrm{Tr}(A{S}_1)$. Positivity follows because $⟨v,S_{1}v⟩=\Vert {\sum }_j⟨v,{\xi }_j⟩{\eta }_j{\Vert }^2\ge 0$, and $\mathrm{Tr}{S}_1=\Vert \Psi {\Vert }^2=:Z$. Thus $\boxed{{\rho }_1=S_1/Z}$ is a density operator; if $\Psi$ is already normalized, $Z=1$. The Gram factors cannot be dropped merely because the vectors are called bases. For unnormalized $\Psi$, physical expectations likewise divide its raw quadratic form by $Z$.

For the Schmidt decomposition diagonalize $S_1$ in an orthonormal basis ${\xi }_j$ with its positive eigenvalues $c_j^2$. Expand $\Psi ={\sum }_j{\xi }_j\otimes w_j$ with $w_j=(⟨{\xi }_j|\otimes I)\Psi$. The reduced matrix gives $⟨w_k,w_j⟩=c_j^2{\delta }_{jk}$. Therefore ${\eta }_j=w_j/c_j$ are orthonormal andZero-eigenvalue components vanish because their squared norms are zero. This proves the decomposition for an arbitrary nonzero vector; the zero vector uses the empty sum. For a unit vector $\sum c_j^2=1$.

In the numerical example the coefficient matrix in the two product bases is $C=\left(\begin{array}{ccc}1& 2& 3\\ 0& 1& 2\end{array}\right)$. Taking the first-subsystem trace yieldsThe vector $(1,-2,1{)}^{\mathsf{T}}$ is in its nullspace, so rank is at most two; its leading $2\times 2$ minor has determinant one, so $\boxed{\mathrm{rank}{S}_2=2}$. This is the rank of $C$ and the Schmidt rank, proving entanglement. The displayed source ket has norm squared $19$, and the displayed reduced matrix has trace $19$. The physical density operator is $\boxed{{\rho }_2=S_2/19}$; its nonzero eigenvalues are $(19\pm \sqrt{337})/38$, since $C{C}^{\mathsf{T}}=\left(\begin{array}{cc}14& 8\\ 8& 5\end{array}\right)$. The raw matrix is correct for the unnormalized quadratic form actually printed.

Put $d=2^N$ and $\lambda =e^{2\pi i/d}$, with nonnegative integer $N$ so the register exists. The geometric sum satisfiesConsequently the matrix $F_{+}$ with entries $d^{-1/2}{\lambda }^{jk}$ has $F_{+}^{\ast }{F}_{+}=I$. This proves unitarity andFor a binary output string $j={\sum }_{r=1}^N{j}_r{2}^{N-r}$, expansion of the printed tensor product gives coefficientThus the printed product is exactly $F_{+}|k⟩$. It is a correct product decomposition of a Fourier phase state, but the printed recovery identity has inconsistent signs:by the square of the quantum Fourier transform. To recover $f(j)$ instead, use $F_{-}=F_{+}^{\ast }$, replacing every product phase by its complex conjugate. Equivalently one may reverse the sign in the definition of $\hat{f}$. The discrepancy is visible for $d\ge 4$ by taking $f$ concentrated at $1$: the printed expression returns the basis vector at $d-1$. The recovery display also uses $k$ as the summation variable on its right while its summand uses $j$; that dummy variable should be $j$.

For factorization first check $\gcd (3,91)=1$. Successive residues areThe first return is $\boxed{r=6}$. Its half power is $3^3=27$, which is neither $1$ nor $-1$ modulo $91$. Therefore

In Shor's algorithm prepare a uniform exponent register of size $Q=2^n$, conveniently $91^2\le Q<2\cdot 91^2$, and reversibly compute $3^x\mathrm{mod}91$ into a second register. Measuring that second register, or tracing it out, leaves periodic exponent correlations of period $r$. Conditional on one measured residue, the first register is a superposition over an arithmetic progression of exponents separated by $r$. The quantum Fourier transform converts this periodic structure into peaks near $k/Q=s/r$. Continued fractions recover $s/r$; when $s$ is coprime to $r$, its denominator is $r$. Verify by modular exponentiation and repeat when the reduced denominator is insufficient. For this example $Q=16384$ is suitable, and useful peaks include those near $Q/6$ and $5Q/6$. Either transform sign gives the same period information, so the source's inverse-sign error does not invalidate the order-finding mechanism.

The Part C inverse-transform display also prints $j$ as its summation index while the summand uses $k$; that sum must be over $k$, as in the corrected inverse above.

Use an expansion for the eigenstate continuing the nondegenerate level, with intermediate normalization $⟨{\psi }_0,\psi (ϵ)⟩=1$:Insert these expressions into the eigenvalue equation. At first order,Taking the inner product with ${\psi }_0$ eliminates the left side by self-adjointness. ThusFor $j\ne 0$, projection onto ${\psi }_j$ instead givesAt second order, projecting the eigenvalue equation onto ${\psi }_0$ gives $E^{(2)}=⟨{\psi }_0,H^{\prime }\chi ⟩$. Expand $\chi$ in the orthonormal energy basis and use self-adjointness of $H^{\prime }$:Nondegeneracy ensures the denominators for this level do not vanish; convergence and validity of the expansion are among the assumed perturbation-theory hypotheses.

For the cosine perturbation of the ground level in a one-dimensional box, write $u=\pi x/a$. The product-to-sum identity givesso, with the given common normalization of the sine eigenfunctions,The diagonal matrix element is $-ϵ/2$, and the only nonzero off-diagonal matrix element out of the ground state is $⟨{\varphi }_3,V{\varphi }_1⟩=ϵ/2$. Let $E_{\mathrm{box}}={\pi }^2{\hslash }^2/(2m{a}^2)$. Since ${\mathcal{E}}_3=9E_{\mathrm{box}}$, the ground energy isThe second-order term is negative, as expected for a nondegenerate ground level: every excited-quantum state denominator is negative. The first-order shift is also negative for positive $ϵ$, because the ground-quantum state density is largest near the center, where this perturbing potential is negative.

Use the Dirac notation convention that the inner product is linear in the ket. Orthonormality of the computational basis givesThe finite geometric sum satisfies ${\sum }_{k=0}^{N-1}{\lambda }^{rk}=N$ if $r=0$ in ${\mathbb{Z}}_N$, and zero otherwise: in the second case it equals $(1-{\lambda }^{rN})/(1-{\lambda }^r)=0$. ThereforeThere are $N$ such orthonormal vectors in an $N$-dimensional space, so they form a basis. The inverse expansion is verified by the same sum:This is the unitary discrete Fourier transform and its inverse.

For the basis-copying isometry, the images $|k⟩|k⟩$ are orthonormal, so $C$ is a linear isometry, with $C^{\ast }C=I$. Its matrix elements giveApplying the definition of the adjoint operator, the coefficient of $|k⟩$ in $C^{\ast }|f⟩|g⟩$ is exactly this expression. ThusEquivalently $C^{\ast }|k⟩|\ell ⟩={\delta }_{k\ell }|k⟩$. The adjoint of basis copying acts as pointwise multiplication of coefficient functions.

For the Fourier transform of a basis-copying isometry, the normalization matters:Hence $C^{\ast }|\hat{m}⟩|\hat{n}⟩=N^{-1/2}|\widehat{m+n}⟩$. Taking adjoint matrix elements, or expanding both output factors in the inverse Fourier basis, yieldsThe printed formula is missing the factor $1/\sqrt{N}$. There are $N$ mutually orthogonal terms in its sum, so without that factor its right side has norm $\sqrt{N}$, whereas its left side has norm one because $C$ is an isometry. The formula above is the normalized identity consistent with all the preceding definitions. For $N=1$ the distinction disappears, but the no-cloning assertion is substantive for $N\ge 2$.

The normalized output is not $|\hat{j}⟩|\hat{j}⟩$. Its coefficient matrix pairs every Fourier label $m$ with $j-m$; tracing out either factor gives $I/N$. Thus it is a maximally entangled state, while two clones of a pure Fourier quantum state would have pure reduced quantum states. A copier for the computational basis therefore does not copy the Fourier basis.

More generally, suppose a linear isometry with a fixed blank input could clone every normalized quantum state. Its restriction to $|\psi ⟩|\mathrm{blank}⟩$ would define a linear map $D$ with $D|k⟩=|k⟩|k⟩$ on every computational basis quantum state. Linearity forces $D=C$, which has just been shown to fail on Fourier quantum states. This proves the no-cloning theorem in the pure-quantum state isometric setting.

The inner-product proof makes the obstruction independent of the chosen basis. If an isometry cloned two quantum states with overlap $s=⟨\psi |\varphi ⟩$, preservation of inner products would give $s=s^2$, or, allowing irrelevant output phases, $|s|=|s{|}^2$. Distinct nonorthogonal rays have $0<|s|<1$, which is impossible. Even allowing an output environment does not help: unitary evolution would give $|s|=|s{|}^2|⟨e_{\psi }|e_{\varphi }⟩|\le |s{|}^2$, the same contradiction. Arbitrary unknown quantum states cannot be cloned, although any fixed orthonormal basis can be copied and a one-dimensional quantum state space is trivial.

Write $u={\psi }_{1,0}=a_1^{\ast }\psi$ and $v={\psi }_{0,1}=a_2^{\ast }\psi$. The canonical commutation relations for the creation and annihilation operators giveFor example, $a_1{a}_1^{\ast }\psi =(I+a_1^{\ast }{a}_1)\psi =\psi$. The expression for $L$ from part (b), with distinct modes reordered, then yieldsThus $\mathrm{span}\{u,v\}$ is an invariant subspace, and the matrix of $L$ in the ordered orthonormal basis $(u,v)$ isIts characteristic polynomial is ${\ell }^2-{\hslash }^2$, giving the eigenvalues and normalized eigenstatesIndeed, $L(u\pm iv)=\pm \hslash (u\pm iv)$ and $\Vert u\pm iv{\Vert }^2=2$ by the assumed orthonormality. Both eigenstates remain in the $N=1$ eigenspace.

In the ordered orthonormal basis $({\psi }_1,{\psi }_2)$, the Hamiltonian operator has matrixwhere ${\sigma }_1$ is a Pauli matrix. Its symmetric and antisymmetric eigenvectors giveThe energy eigenstates ${\chi }_{\pm }$ are normalized and orthogonal; each energy eigenvalue is nondegenerate. Their full stationary states are $e^{-i\omega t}{\chi }_{+}$ and $e^{i\omega t}{\chi }_{-}$. Since they form an orthonormal basis, the general solution to the Time-dependent Schrödinger equation isFor a normalized quantum state, $|c_{+}{|}^2+|c_{-}{|}^2=1$; the constants are determined by its initial inner products with ${\chi }_{\pm }$.

Diagonalizing the $J_1$ matrix gives the normalized eigenvectorswith eigenvalues $+\hslash ,0,-\hslash$, respectively. The first projective measurement prepares $v_{+}$ up to global phase, regardless of the earlier state, conditional on obtaining $+\hslash$.

The unitary time evolution under $H=\omega J_3$ gives, with $\theta =\omega t$,The probability amplitudes for the second projective measurement areThe Born rule therefore gives the general spin-one transverse measurement after axial precession:At the specified times $t_k=\pi k/\omega$, $\theta =\pi k$, soFor even $k$, the result is certainly $+\hslash$; for odd $k$, it is certainly $-\hslash$. Although zero belongs to the spectrum of the observable, it never occurs at these times. Indeed, for even $k$ the evolved vector is $v_{+}$, while for odd $k$ it is $-v_{-}$, where the minus sign is a global phase.

Taking the transpose gives $M^T=(X^T+X)/\sqrt{2}=M$, so $M$ is a symmetric matrix. For $i<j$, $M_{ij}=(X_{ij}+X_{ji})/\sqrt{2}$ has the standard normal distribution, whereas $M_{ii}=\sqrt{2}{X}_{ii}$ has the normal distribution $N(0,2)$. Distinct upper-triangular entries depend on disjoint sets of the original independent random variables, so they remain independent.

Consequently the joint probability density function, relative to ${\prod }_{i\le j}d{m}_{ij}$, isHere $\mathrm{tr}(m^2)={\sum }_i{m}_{ii}^2+2{\sum }_{i<j}{m}_{ij}^2$ because $m$ is a symmetric matrix. This is the Gaussian orthogonal ensemble, with off-diagonal variance one and diagonal variance two; the rescaled matrix $M/\sqrt{n}$ has the corresponding $1/n$ off-diagonal variance.

We will also use its entry covariance. Expanding in the independent entries of $X$, for which $\mathbb{E}[X_{ij}{X}_{kl}]={\delta }_{ik}{\delta }_{jl}$, givesThe symbols ${\delta }_{ab}$ are Kronecker deltas. This formula includes the diagonal case, where the two terms both contribute.

The spectral theorem for real symmetric matrices identifies the sum of the $k$th powers of the eigenvalues with the trace of the $k$th matrix power. Therefore the moment of the empirical spectral measure isFor odd $k$, every summand has zero expectation by part (a)(ii), so $\boxed{\mathbb{E}[m_k(n)]=0}$ for every $n$.

For the second moment, use $M_{ji}=M_{ij}$. There are $n$ diagonal entries with variance two and $n(n-1)$ ordered off-diagonal entries with variance one. Their means vanish, givingHence $\boxed{\underset{n\to \infty }{\lim }\mathbb{E}[m_2(n)]=1}$.

First $W^T=W$, so the spectral theorem for real symmetric matrices gives real eigenvalues and an orthonormal basis of eigenvectors. Moreover, for every real vector $v$,Thus $W$ is a positive semidefinite matrix. For any nonzero eigenvector $v$ with eigenvalue $\lambda$,This matrix is a scaled example of the Wishart distribution.

Let $P$ be the $3$-by-$3$ evaluation matrix $P_{i,\ell +1}=p_{\ell }({\lambda }_i)$. The kernel evaluation matrix is the Gram matrix $P{P}^T$, so multiplicativity of the determinant givesBecause the polynomials are monic of consecutive degrees, their coefficient matrix relative to $1,x,x^2$ is a triangular matrix with unit diagonal. It has determinant one, so the Vandermonde determinant invariance under monic basis change givesSquaring removes the ordering sign. Thus the Gaussian unitary ensemble joint probability density function has precisely the kernel formNo orthogonality of the polynomials is needed for this step.

Use Haar sampling by Gaussian QR decomposition:

The matrix $G$ is invertible almost surely. Indeed, successively sampled Gaussian columns have probability zero of lying in the proper linear subspace spanned by the previous columns. Thus the diagonal signs used above exist with probability one. Equivalently, apply the Gram-Schmidt process to the columns, always dividing each residual by its positive Euclidean norm; this produces $M$ directly.

To prove uniformity, fix $H\in O(n)$. Each column of $G$ has density proportional to $e^{-\Vert x{\Vert }^2/2}$, so orthogonal transformations preserve its probability distribution. Column independence therefore gives $HG\stackrel{\mathrm{d}}{=}G$.

The QR decomposition with positive diagonal is unique: if $G=Q_1{R}_1=Q_2{R}_2$, then $Q_2^T{Q}_1=R_2{R}_1^{-1}$ is both an orthogonal matrix and an upper triangular matrix, hence diagonal with entries $\pm 1$. Its diagonal entries are positive, so it is the identity. It follows that the orthogonal factor of $HG=HMR$ is $HM$. ThereforeBy uniqueness of normalized Haar measure on the compact group $O(n)$, the returned matrix has exactly the required probability distribution. Reflection invariance gives each determinant sign probability $1/2$, so the algorithm samples both components of the orthogonal group. The positive-diagonal correction is essential when a numerical QR decomposition chooses its own signs.

The exact orders printed in the question are inconsistent. The source PDF says that the orders of group elements $a,b,c$ are $4,12,24$, respectively, and also imposes $3a+2b=0$. Multiplying this relation by $6$ and using $12b=0$ gives $18a=0$. Since $4a=0$ also holds, subtracting $4(4a)=0$ gives $2a=0$, contradicting the stated order of a group element $4$. Thus

The natural abelian group presentation interpretation is to impose the relations $4a=12b=24c=0$ and $3a+2b=0$, with no further relations, without requiring the orders of the images of the generators to stay $4,12,24$. Under this interpretation the calculation is as follows. Since $3a=-a$ follows from $4a=0$, the mixed relation gives $a=2b$. Consequently $8b=0$, and together with $12b=0$ this gives $4b=0$. Conversely, $a=2b$, $4b=0$, and $24c=0$ imply every original relation. Eliminating $a$ leaves two independent cyclic groups:To verify that these relations force no further identifications, map $a,b,c$ to $(2,0),(1,0),(0,1)$ in ${\mathbb{Z}}_4\oplus {\mathbb{Z}}_{24}$. All defining relations hold, and the inverse group homomorphism sends the two standard generators to $b,c$. These group homomorphisms are mutually inverse. The divisibility condition $4\mid 24$ gives the requested invariant factors of a finitely generated module over $\mathbb{Z}$. In this presented group, $a,b,c$ actually have orders $2,4,24$.

Equivalently, the relation matrixhas nonzero Smith normal form entries $1,4,24$, with the trivial ${\mathbb{Z}}_1$ factor omitted. The displayed direct sum of modules is the answer to the presentation interpretation; the inconsistency of the literal exact-order formulation remains.

Let $v\in {\mathbb{Q}}^2$ be nonzero. If the columns $v,Av$ were linearly dependent over the rational numbers, then $Av=\lambda v$ for some $\lambda \in \mathbb{Q}$. Applying $A$ once more would giveso $({\lambda }^2+1)v=0$. Since $v\ne 0$, this forces ${\lambda }^2=-1$, impossible in $\mathbb{Q}$. Thus the two columns are linearly independent, and the square matrixis an invertible matrix.

The nonconstant invariant factors of a linear operator satisfy $d_1\mid \cdots \mid d_k=m_A$. By part (a)(i), $m_A=X^2+1$ is an irreducible polynomial, so every $d_j$ must equal $X^2+1$. Their polynomial degrees sum to $n$, giving $n=2k$.

The companion matrix of $X^2+1$, in the basis $(1,X)$ of $\mathbb{Q}[X]/(X^2+1)$, isTherefore, for each fixed even $n$, every such $A$ has rational canonical formThere are no such matrices when $n$ is odd. This is an instance of the rational canonical form with irreducible minimal polynomial principle; equality of minimal polynomials alone does not determine rational canonical form in general.

Regard $A$ as a matrix over the rational numbers. Part (a)(ii), with $n=2$, gives its characteristic polynomial $X^2+1$. Therefore, writing $i$ for the imaginary unit,The matrix $M=iI-A$ is not zero: otherwise $A=iI$ would have noninteger diagonal entries. It consequently has rank of a matrix one over the fraction field $\mathbb{Q}(i)$ of the Gaussian integers.

By part (b)(i), $\mathbb{Z}[i]$ is a euclidean domain, hence a principal ideal domain, so apply Smith normal form. There exist $P,Q\in {\mathrm{GL}}_2(\mathbb{Z}[i])$ such thatIndeed, invertible multiplication preserves rank of a matrix, so exactly one diagonal entry is nonzero. For $e_2=(0,1{)}^{\mathsf{T}}$, multiplication on the right gives $P(iI-A)Q{e}_2=0$. Since $P$ is an invertible matrix,The Smith normal form supplies a kernel basis lemma ensures that the null vector $Q{e}_2$ is a column of an invertible matrix over $\mathbb{Z}[i]$; an arbitrary nonzero solution over $\mathbb{Q}(i)$ would not itself establish that requirement.

Over a principal ideal domain $D$, an $m$-by-$n$ matrix is in Smith normal form if all entries off the main diagonal vanish and its nonzero diagonal entries $d_1,\dots ,d_r$ occur first and satisfyAll remaining diagonal entries are zero. The zero matrix is allowed, with $r=0$. For an arbitrary matrix $M$, the Smith normal form theorem gives invertible matrices $P\in {\mathrm{GL}}_m(D)$ and $Q\in {\mathrm{GL}}_n(D)$ withwhere the right side is interpreted as a rectangular diagonal matrix when necessary. The $d_j$ are determined up to multiplication by units. Positivity can be required over $\mathbb{Z}$, but is not a condition over the Gaussian integers.

Each entry of the matrix $Q$ from part (b)(iii) has unique real and imaginary integer coefficients. WriteIt already belongs to the general linear group ${\mathrm{GL}}_2(\mathbb{Z}[i])$. Expanding the relation from part (b)(iii) givesThe real and imaginary parts must both vanish, soThese are the required identities, with $e_2=(0,1{)}^{\mathsf{T}}$.