quantas seismic

The SEISMIC frontend solves the Christoffel eigenproblem over a spherical grid. The selected calculation level determines whether the archive contains only phase fields, first-derivative group properties, or second-derivative enhancement and caustic diagnostics.

Key controls

  • --hemisphere changes the sampled angular interval. The same ntheta gives a coarser spacing on a full sphere than on an upper hemisphere.

  • --level controls scientific content and archive size.

  • --ntheta and --nphi control angular resolution and therefore sampled extrema and candidate locations.

  • --batch-size controls vectorized throughput and temporary memory only; its default is 512 directions and it is not an accuracy parameter.

  • eigenvalue and degeneracy tolerances classify near-zero numerical cases. They should not be enlarged to hide an unstable stiffness tensor.

  • polarization tracking changes branch continuity and visualization, not phase velocities.

See Seismic-wave analysis: implementation and workflow for the detailed numerical workflow, Seismic-wave analysis for the hydroxylapatite example, and Elasticity and SEISMIC text input for the shared tensor input.

Generated command reference

quantas seismic

Calculate, persist, export, and plot anisotropic acoustic-wave properties.

The workflow solves the Christoffel problem on a spherical grid and can progressively add group velocity, power-flow, enhancement, caustic, and polarization-tracking data. The HDF5 archive is the common source for reports, CSV exports, and figures.

Usage

quantas seismic [OPTIONS] COMMAND [ARGS]...

export

Export sampled acoustic fields from a Quantas SEISMIC HDF5 result.

The CSV contains one row per sampled direction and wave mode, including every phase, group, enhancement, degeneracy, and tracking field available in FILENAME. No new Christoffel calculation is performed.

Usage

quantas seismic export [OPTIONS] FILENAME

Options

-o, --output <outfile>

Output CSV file. Default: input base name + ‘.csv’. The command reads an existing result and does not rerun the scientific workflow.

Arguments

FILENAME

Required argument

plot

Render spherical maps, acoustic surfaces, or a six-panel summary from HDF5.

FILENAME must contain the property level required by each requested plot. Map, surface, polarization, projection, and typography options affect visualization only; they do not change the stored wave solution.

Usage

quantas seismic plot [OPTIONS] FILENAME

Options

-o, --output <outfile>

Base name for output figures. Default: HDF5 base name. This affects the generated figure only and does not modify the stored scientific result.

--2d

Render spherical 2D maps. This affects the generated figure only and does not modify the stored scientific result.

--3d

Render 3D acoustic surfaces. This affects the generated figure only and does not modify the stored scientific result.

--summary

Render the six-panel seismic summary. This affects the generated figure only and does not modify the stored scientific result.

--property <properties>

Scalar property for 2D or 3D output. Repeat to select multiple properties. This affects the generated figure only and does not modify the stored scientific result.

Options:

phase_v_p | phase_v_s1 | phase_v_s2 | shear_anisotropy | shear_splitting | phase_v_p_over_v_s1 | phase_v_p_over_v_s2 | group_v_p | group_v_s1 | group_v_s2 | power_flow_v_p | power_flow_v_s1 | power_flow_v_s2 | log10_enhancement_v_p | log10_enhancement_v_s1 | log10_enhancement_v_s2

--surface <surface_types>

3D surface geometry. Repeat to select multiple types. This affects the generated figure only and does not modify the stored scientific result.

Options:

phase | slowness | group

--mode <modes>

Acoustic mode used for 3D surfaces. Repeat to select multiple modes. This affects the generated figure only and does not modify the stored scientific result.

Options:

v_p | v_s1 | v_s2

--projection <projection>

Spherical projection used for two-dimensional maps. This affects the generated figure only and does not modify the stored scientific result.

Default:

'equal_area'

Options:

equal_area | stereographic

--cmap <cmap>

Override the property-dependent default with a Matplotlib colormap. This affects the generated figure only and does not modify the stored scientific result.

--levels <levels>

Number of filled-contour levels used for two-dimensional maps. This affects the generated figure only and does not modify the stored scientific result.

Default:

16

--extrema, --no-extrema

Show or hide minimum and maximum markers. This affects the generated figure only and does not modify the stored scientific result.

Default:

True

--polarization, --no-polarization

Also render separate polarization-overlay variants whenever a tracked polarization field is available. This affects the generated figure only and does not modify the stored scientific result.

Default:

True

--polarization-stride <polarization_stride>

Sampling stride between plotted polarization axes. This affects the generated figure only and does not modify the stored scientific result.

Default:

8

--polarization-color <polarization_color>

Color of plotted polarization axes. This affects the generated figure only and does not modify the stored scientific result.

Default:

'black'

--polarization-linewidth <polarization_linewidth>

Line width of plotted polarization axes. This affects the generated figure only and does not modify the stored scientific result.

Default:

1.2

--polarization-scale <polarization_scale>

Length scale of plotted polarization axes. This affects the generated figure only and does not modify the stored scientific result.

Default:

0.07

--geometry <geometry>

Use a unit sphere or the natural physical acoustic geometry in 3D. This affects the generated figure only and does not modify the stored scientific result.

Default:

'unit-sphere'

Options:

unit-sphere | physical

--axis-labels <axis_labels>

Use x/y/z labels or [100]/[010]/[001] tensor-frame labels. This affects the generated figure only and does not modify the stored scientific result.

Default:

'cartesian'

Options:

cartesian | crystal

--complete-surface, --sampled-domain

Complete hemispherical 3D results by antipodal symmetry. This affects the generated figure only and does not modify the stored scientific result.

Default:

True

--preset <figure_preset>

Choose a named figure style. ‘screen’ favors interactive viewing, ‘publication’ uses print-oriented geometry and resolution, and ‘monochrome’ adds grayscale-safe rendering; the preset affects presentation only.

Default:

'screen'

Options:

screen | publication | monochrome

--format <image_format>

Choose the static figure file format supported by the Matplotlib renderer. The format changes serialization only; vector formats are preferable when downstream editing is required.

Default:

'png'

Options:

png | pdf | svg

--dpi <dpi>

Override the raster resolution supplied by the selected figure preset. This affects PNG or other raster output only and does not change sampled scientific data.

--show

Display figures interactively after files are generated. Leave this disabled on headless, batch, or continuous-integration systems.

Arguments

FILENAME

Required argument

run

Run the sampled SEISMIC workflow and save a native HDF5 result.

FILENAME contains a stiffness tensor and density. The angular grid determines spatial resolution, ‘–level’ determines the highest derivative-dependent property computed, and ‘–batch-size’ controls memory and throughput only. A mechanically unstable tensor is rejected before wave propagation is evaluated.

Usage

quantas seismic run [OPTIONS] FILENAME

Options

--hemisphere <hemisphere>

Sample the selected spherical hemisphere. This changes the physical model or interpretation and should be compared through the relevant scientific diagnostics.

Default:

'upper'

Options:

upper | lower | full

-L, --level <level>

Highest acoustic-property level calculated and persisted. This changes the physical model or interpretation and should be compared through the relevant scientific diagnostics.

Default:

'enhancement'

Options:

phase | group | enhancement

--track-polarizations, --no-track-polarizations

Track shear-wave polarization axes across the sampled sphere. This changes the physical model or interpretation and should be compared through the relevant scientific diagnostics.

Default:

True

--ntheta <ntheta>

Set the polar sampling count for the selected hemisphere. The default 91 gives about one-degree spacing on the upper hemisphere; increase it to converge narrow extrema or acoustic-axis candidates.

Default:

91

--nphi <nphi>

Set the azimuthal sampling count. The default 181 avoids a duplicate periodic endpoint while matching the angular scale of the default upper-hemisphere polar grid.

Default:

181

--batch-size <batch_size>

Set the number of directions solved in each vectorized numerical batch. The default 512 balances temporary memory and overhead; changing it affects performance and progress granularity, not the angular grid or numerical accuracy.

Default:

512

-o, --output <output>

HDF5 result file. Default: input base name + ‘_seismic.hdf5’. This controls persistence or presentation and does not change the scientific model.

-f, --force

Replace existing generated files without an interactive confirmation. Use this deliberately in reproducible scripts because overwritten results cannot be recovered by Quantas.

-r, --report <report>

Write a deterministic plain-text scientific report. When omitted, run commands use the primary input name with a ‘.log’ suffix; reports contain no ANSI styling or live progress.

-v, --verbosity <verbosity>

Select report detail: ‘standard’ gives the normal scientific summary, ‘extended’ adds complete tables, and ‘debug’ adds numerical diagnostics. The choice never changes stored numerical results.

Default:

'standard'

Options:

standard | extended | debug

-q, --quiet

Suppress terminal presentation while preserving calculation, HDF5 output, and the report. This is useful for scripts, redirected jobs, and automated validation.

--progress, --no-progress

Enable or disable transient Rich progress on an interactive terminal. Progress events are operational only and are never written to reports or HDF5 histories.

Default:

True

--rotate-xyz <X Y Z>

Rotate stiffness components into a new analysis frame using right-handed fixed-axis angles about source X, then Y, then Z. Wave directions are reported in the rotated frame while invariant speeds are preserved.

--rotation-matrix <rotation_matrix>

Read a proper 3 x 3 source-to-analysis component-transformation matrix. This changes the coordinate frame used for directions and polarizations, not the physical crystal.

Arguments

FILENAME

Required argument