Submit an analog program using QuEra’s Aquila - Amazon Braket

# Submit an analog program using QuEra’s Aquila

This page provides a comprehensive documentation about the capabilities of the Aquila machine from QuEra. Details covered here are the following: 1) The parameterized Hamiltonian simulated by Aquila, 2) AHS program parameters, 3) AHS result content, 4) Aquila capabilities parameter. We suggest using Ctrl+F text search to find parameters relevant to your questions.

## Hamiltonian

The Aquila machine from QuEra simulates the following (time-dependent) Hamiltonian natively:

###### Note

Access to local detuning is an Experimental capability and is available by request through Braket Direct.

where

• Hdrive,k​(t)=( 1/2 ​Ω(t)eiϕ(t)S−,k​ + 1/2 ​Ω(t)e−iϕ(t) S+,k​) + (−Δglobal​(t)nk​),

• Ω(t) is the time-dependent, global driving amplitude (aka Rabi frequency), in units of (rad / s)

• ϕ(t) is the time-dependent, global phase, measured in radians

• S−,k​ and S+,k​ are the spin lowering and raising operators of atom k (in the basis |↓⟩=|g⟩, |↑⟩=|r⟩, they are S​=|g⟩⟨r|, S+​=(S​)=|r⟩⟨g|)

• Δglobal​(t) is the time-dependent, global detuning

• nk ​is the projection operator on the Rydberg state of atom k (i.e. n=|r⟩⟨r|)

• Hlocal detuning,k(t)=-Δlocal(t)hknk

• Δlocal(t) is the time-dependent factor of the local frequency shift, in units of (rad / s)

• hk is the site-dependent factor, a dimensionless number between 0.0 and 1.0

• Vvdw,k,l​=C6​/(dk,l​)6nk​nl​,

• C6​ is the van der Waals coefficient, in units of (rad / s) * (m)^6

• dk,l ​is the Euclidean distance between atom k and l, measured in meters.

Users have control over the following parameters through the Braket AHS program schema.

• 2-d atom arrangement (xk​ and yk​ coordinates of each atom k, in units of um), which controls the pairwise atomic distances dk,l​ with k,l=1,2,…N

• Ω(t), the time-dependent, global Rabi frequency, in units of (rad / s)

• ϕ(t), the time-dependent, global phase, in units of (rad)

• Δglobal​(t), the time-dependent, global detuning, in units of (rad / s)

• Δlocal(t), the time-dependent (global) factor of the magnitude of local detuning, in units of (rad / s)

• hk, the (static) site-dependent factor of the magnitude of local detuning, a dimensionless number between 0.0 and 1.0

###### Note

The user cannot control which levels are involved (i.e. S​,S+​, n operators are fixed) nor the strength of the Rydberg-Rydberg interaction coefficient (C6​).

## Braket AHS program schema

braket.ir.ahs.program_v1.Program object (example)

``````Program(
name='braket.ir.ahs.program',
version='1'
),
setup=Setup(
ahs_register=AtomArrangement(
sites=[
[Decimal('0'), Decimal('0')],
[Decimal('0'), Decimal('4e-6')],
[Decimal('4e-6'), Decimal('0')],
],
filling=[1, 1, 1]
)
),
hamiltonian=Hamiltonian(
drivingFields=[
DrivingField(
amplitude=PhysicalField(
time_series=TimeSeries(
values=[Decimal('0'), Decimal('15700000.0'), Decimal('15700000.0'), Decimal('0')],
times=[Decimal('0'), Decimal('0.000001'), Decimal('0.000002'), Decimal('0.000003')]
),
pattern='uniform'
),
phase=PhysicalField(
time_series=TimeSeries(
values=[Decimal('0'), Decimal('0')],
times=[Decimal('0'), Decimal('0.000001')]
),
pattern='uniform'
),
detuning=PhysicalField(
time_series=TimeSeries(
values=[Decimal('-54000000.0'), Decimal('54000000.0')],
times=[Decimal('0'), Decimal('0.000001')]
),
pattern='uniform'
)
)
],
localDetuning=[
LocalDetuning(
magnitude=PhysicalField(
times_series=TimeSeries(
values=[Decimal('0'), Decimal('25000000.0'), Decimal('25000000.0'), Decimal('0')],
times=[Decimal('0'), Decimal('0.000001'), Decimal('0.000002'), Decimal('0.000003')]
),
pattern=Pattern([Decimal('0.8'), Decimal('1.0'), Decimal('0.9')])
)
)
]
)
)``````

JSON (example)

``````{
"name": "braket.ir.ahs.program",
"version": "1"
},
"setup": {
"ahs_register": {
"sites": [
[0E-7, 0E-7],
[0E-7, 4E-6],
[4E-6, 0E-7],
],
"filling": [1, 1, 1]
}
},
"hamiltonian": {
"drivingFields": [
{
"amplitude": {
"time_series": {
"values": [0.0, 15700000.0, 15700000.0, 0.0],
"times": [0E-9, 0.000001000, 0.000002000, 0.000003000]
},
"pattern": "uniform"
},
"phase": {
"time_series": {
"values": [0E-7, 0E-7],
"times": [0E-9, 0.000001000]
},
"pattern": "uniform"
},
"detuning": {
"time_series": {
"values": [-54000000.0, 54000000.0],
"times": [0E-9, 0.000001000]
},
"pattern": "uniform"
}
}
],
"localDetuning": [
{
"magnitude": {
"time_series": {
"values": [0.0, 25000000.0, 25000000.0, 0.0],
"times": [0E-9, 0.000001000, 0.000002000, 0.000003000]
},
"pattern": [0.8, 1.0, 0.9]
}
}
]
}
}``````
Main fields
Program field type description

setup.ahs_register.sites

List[List[Decimal]]

List of 2-d coordinates where the tweezers trap atoms

setup.ahs_register.filling

List[int]

Marks atoms that occupy the trap sites with 1, and empty sites with 0

hamiltonian.drivingFields[].amplitude.time_series.times

List[Decimal]

time points of driving amplitude, Omega(t)

hamiltonian.drivingFields[].amplitude.time_series.values

List[Decimal]

values of driving amplitude, Omega(t)

hamiltonian.drivingFields[].amplitude.pattern

str

spatial pattern of driving amplitude, Omega(t); must be 'uniform'

hamiltonian.drivingFields[].phase.time_series.times

List[Decimal]

time points of driving phase, phi(t)

hamiltonian.drivingFields[].phase.time_series.values

List[Decimal]

values of driving phase, phi(t)

hamiltonian.drivingFields[].phase.pattern

str

spatial pattern of driving phase, phi(t); must be 'uniform'

hamiltonian.drivingFields[].detuning.time_series.times

List[Decimal]

time points of driving detuning, Delta_global(t)

hamiltonian.drivingFields[].detuning.time_series.values

List[Decimal]

values of driving detuning, Delta_global(t)

hamiltonian.drivingFields[].detuning.pattern

str

spatial pattern of driving detuning, Delta_global(t); must be 'uniform'

hamiltonian.localDetuning[].magnitude.time_series.times

List[Decimal]

time points of the time-dependent factor of the local detuning magnitude, Delta_local(t)

hamiltonian.localDetuning[].magnitude.time_series.values

List[Decimal]

values of the time-dependent factor of the local detuning magnitude, Delta_local(t)

hamiltonian.localDetuning[].magnitude.pattern

List[Decimal]

site-dependent factor of the local detuning magnitude, h_k (values corresponds to sites in setup.ahs_register.sites)

Program field type description

str

name of the schema; must be 'braket.ir.ahs.program'

str

version of the schema

## Braket AHS task result schema

``````AnalogHamiltonianSimulationQuantumTaskResult(
version='1'
),
shots=2,
deviceId='arn:aws:braket:us-east-1::device/qpu/quera/Aquila',
deviceParameters=None,
createdAt='2022-10-25T20:59:10.788Z',
endedAt='2022-10-25T21:00:58.218Z',
status='COMPLETED',
failureReason=None
),
measurements=[
ShotResult(
status=<AnalogHamiltonianSimulationShotStatus.SUCCESS: 'Success'>,

pre_sequence=array([1, 1, 1, 1]),
post_sequence=array([0, 1, 1, 1])
),

ShotResult(
status=<AnalogHamiltonianSimulationShotStatus.SUCCESS: 'Success'>,

pre_sequence=array([1, 1, 0, 1]),
post_sequence=array([1, 0, 0, 0])
)
]
)``````

JSON (example)

``````{
"version": "1"
},
"version": "1"
},
"shots": 2,
"deviceId": "arn:aws:braket:us-east-1::device/qpu/quera/Aquila",

"createdAt": "2022-10-25T20:59:10.788Z",
"endedAt": "2022-10-25T21:00:58.218Z",
"status": "COMPLETED"

},
"measurements": [
{
"shotResult": {
"preSequence": [1, 1, 1, 1],
"postSequence": [0, 1, 1, 1]
}
},
{
"shotResult": {
"preSequence": [1, 1, 0, 1],
"postSequence": [1, 0, 0, 0]
}
}
],
"action": {...}
"version": "1"
},
"numSuccessfulShots": 100
}
}
}``````
Main fields

measurements[].shotResult.preSequence

List[int]

Pre-sequence measurement bits (one for each atomic site) for each shot: 0 if site is empty, 1 if site is filled, measured before the sequences of pulses that run the quantum evolution

measurements[].shotResult.postSequence

List[int]

Post-sequence measurement bits for each shot: 0 if atom is in Rydberg state or site is empty, 1 if atom is in ground state, measured at the end of the sequences of pulses that run the quantum evolution

str

str

version of the schema

str

str

version of the schema

str

int

The number of shots for the quantum task

str

The ID of the device on which the quantum task ran. For AWS devices, this is the device ARN.

str

The timestamp of creation; the format must be in ISO-8601/RFC3339 string format YYYY-MM-DDTHH:mm:ss.sssZ. Default is None.

str

The timestamp of when the quantum task ended; the format must be in ISO-8601/RFC3339 string format YYYY-MM-DDTHH:mm:ss.sssZ. Default is None.

str

The status of the quantum task (CREATED, QUEUED, RUNNING, COMPLETED, FAILED). Default is None.

str

The failure reason of the quantum task. Default is None.

braket.ir.ahs.program_v1.Program

(See the Braket AHS program schema section)

str

str

version of the schema

int

number of completely successful shots; must be equal to the requested number of shots

int

The status of the shot, (Success, Partial success, Failure); must be "Success"

## QuEra device properties schema

``````QueraDeviceCapabilities(
service=DeviceServiceProperties(
name='braket.device_schema.device_service_properties',
version='1'
),
executionWindows=[
DeviceExecutionWindow(
executionDay=<ExecutionDay.MONDAY: 'Monday'>,
windowStartHour=datetime.time(1, 0),
windowEndHour=datetime.time(23, 59, 59)
),
DeviceExecutionWindow(
executionDay=<ExecutionDay.TUESDAY: 'Tuesday'>,
windowStartHour=datetime.time(0, 0),
windowEndHour=datetime.time(12, 0)
),
DeviceExecutionWindow(
executionDay=<ExecutionDay.WEDNESDAY: 'Wednesday'>,
windowStartHour=datetime.time(0, 0),
windowEndHour=datetime.time(12, 0)
),
DeviceExecutionWindow(
executionDay=<ExecutionDay.FRIDAY: 'Friday'>,
windowStartHour=datetime.time(0, 0),
windowEndHour=datetime.time(23, 59, 59)
),
DeviceExecutionWindow(
executionDay=<ExecutionDay.SATURDAY: 'Saturday'>,
windowStartHour=datetime.time(0, 0),
windowEndHour=datetime.time(23, 59, 59)
),
DeviceExecutionWindow(
executionDay=<ExecutionDay.SUNDAY: 'Sunday'>,
windowStartHour=datetime.time(0, 0),
windowEndHour=datetime.time(12, 0)
)
],
shotsRange=(1, 1000),
deviceCost=DeviceCost(
price=0.01,
unit='shot'
),
deviceDocumentation=
DeviceDocumentation(
summary='Analog quantum processor based on neutral atom arrays',
externalDocumentationUrl='https://www.quera.com/aquila'
),
deviceLocation='Boston, USA',
updatedAt=datetime.datetime(2024, 1, 22, 12, 0, tzinfo=datetime.timezone.utc),
),
action={
<DeviceActionType.AHS: 'braket.ir.ahs.program'>: DeviceActionProperties(
version=['1'],
actionType=<DeviceActionType.AHS: 'braket.ir.ahs.program'>
)
},
deviceParameters={},
name='braket.device_schema.quera.quera_device_capabilities',
version='1'
),
...
...
)
)``````

JSON (example)

``````{
"service": {
"name": "braket.device_schema.device_service_properties",
"version": "1"
},
"executionWindows": [
{
"executionDay": "Monday",
"windowStartHour": "01:00:00",
"windowEndHour": "23:59:59"
},
{
"executionDay": "Tuesday",
"windowStartHour": "00:00:00",
"windowEndHour": "12:00:00"
},
{
"executionDay": "Wednesday",
"windowStartHour": "00:00:00",
"windowEndHour": "12:00:00"
},
{
"executionDay": "Friday",
"windowStartHour": "00:00:00",
"windowEndHour": "23:59:59"
},
{
"executionDay": "Saturday",
"windowStartHour": "00:00:00",
"windowEndHour": "23:59:59"
},
{
"executionDay": "Sunday",
"windowStartHour": "00:00:00",
"windowEndHour": "12:00:00"
}
],
"shotsRange": [
1,
1000
],
"deviceCost": {
"price": 0.01,
"unit": "shot"
},
"deviceDocumentation": {
"summary": "Analog quantum processor based on neutral atom arrays",
"externalDocumentationUrl": "https://www.quera.com/aquila"
},
"deviceLocation": "Boston, USA",
"updatedAt": "2024-01-22T12:00:00+00:00"
},
"action": {
"braket.ir.ahs.program": {
"version": [
"1"
],
"actionType": "braket.ir.ahs.program"
}
},
"deviceParameters": {},
"name": "braket.device_schema.quera.quera_device_capabilities",
"version": "1"
},
...
# See Aquila device page > "Calibration" tab > "JSON" page
...
}
}``````
Service properties fields
Service properties field type description

service.executionWindows[].executionDay

ExecutionDay

Days of the execution window; must be 'Everyday', 'Weekdays', 'Weekend', 'Monday', 'Tuesday', 'Wednesday', Thursday', 'Friday', 'Saturday' or 'Sunday'

service.executionWindows[].windowStartHour

datetime.time

UTC 24-hour format of the time when the execution window starts

service.executionWindows[].windowEndHour

datetime.time

UTC 24-hour format of the time when the execution window ends

service.qpu_capabilities.service.shotsRange

Tuple[int, int]

Minimum and maximum number of shots for the device

service.qpu_capabilities.service.deviceCost.price

float

Price of the device in terms of US dollars

service.qpu_capabilities.service.deviceCost.unit

str

unit for charging the price, e.g: 'minute', 'hour', 'shot', 'task'

action[].version

str

version of the AHS program schema

action[].actionType

ActionType

AHS program schema name; must be 'braket.ir.ahs.program'

str

name of the schema; must be 'braket.device_schema.device_service_properties'

str

version of the schema

service.deviceDocumentation.imageUrl

str

URL for the image of the device

service.deviceDocumentation.summary

str

brief description on the device

service.deviceDocumentation.externalDocumentationUrl

str

external documentation URL

service.deviceLocation

str

geographic location fo the device

service.updatedAt

datetime

time when the device properties were last updated